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interactions,  transverse  effects  tend  to  be  dominant,  and  it  appears  that  such 
ion-electron  interactions  are  (i)  probably  unsuitable  as  a wave  growth/ion  load- 
r ing  mechanism,  but  (ii)  unlikely  to  disrupt  the  process  of  collective  accelera- 

tion. A stability  analysis  of  a two-dimensional  sheet  beam  of  electrons  was 
completed.  We  have  also  determined  accurately  the  limiting  currents  for  magnet- 
ized injection  of  an  unneutralized  relativistic  electron  beam  into  a cylindrical 
drift  tube.  The  previously  cited  interpolation  formula  was  found  to  underesti- 
mate the  true  limiting  current  by  over  20$  in  some  cases,  and  the  corresponding 
expression  for  the  axial  electron  energy  sometimes  overestimates  the  ture  limit- 
ing energy  by  more  than  an  order  of  magnitude.  The  influence  of  a spat ia I . var i a1 
tion  in  y on  the  structure  of  the  relativistic  electron  beam  cyclotron  mode  was 
examined,  and  found  to  increase  the  peakedness  of  the  radial  electric  fields  nea 
the  beam  edge  when  y(r)  increases.  The  influence  of  this  radial  variation  in  y 
on  the  longitudinal  plasma  mode  was  also  examined  in  the  strong  magnetic  field 
limit.  Previous  results,  that  the  plasma  wave  phase  velocity  could  only  be 
made  small  (for  ion  loading,  prior  to  collective  acceleration)  by  increasing  the 
beam  current  nearly  to  the  limiting  current,  were  extended  from  the  pencil  beam 
regime  to  arbitrary  values  of  beam  and  waveguide  radii  and  arbitrary  density  pro- 
files. The  plasma  wave  eigenfunctions  were  found  to  have  increased  peakedness 
near  the  beam  axis  as  the  limiting  current  conditions  are  approached.  An  itera- 
tive calculation  was  performed  of  the  slowly  precess-ng  fluid  beam  equilibrium 
which  is  achieved  by  a 9-symmetric,  steady  state,  initially  solid  electron  beam 
injected  into  a cylindrical  waveguide  along  a finite  magnetic  guide  field.  This 
calculation  incorporates  the  space  charge  slowing  associated  with  the  waveguide 
injection,  as  well  as  the  acquisition  of  precessional  motion  and  associated  dia- 
magnetism, centrifugal  and  Coriolis  forces.  A generalization  of  the  limiting 
current  formula  was  obtained  to  incorporate  the  finite  axial  magnetic  field 
effects — reducing  Bz  increases  1^.  Several  small  parameters  were  isolated  durin' 
the  equilibrium  calculation,  which  may  be  used  to  quantify  the  validity  of  the 
common  assumptions  of  (i)  infinite  magnetic  fields,  or  (ii)  u 1 trare lat i v i st ic 
electron  beams.  An  examination  of  linearized  collective  beam  waves  was  com- 
pleted for  the  important  case  in  which  the  parameter  character izing  centrifugal 
and  Coriolis  forces,  and  plasma-cyclotron  mode  coupling,  was  retained  in  the 
analysis.  The  resulting  fourth-order  radial  eigenvalue  problem  was  solved  in 
terms  of  Bessel  functions,  and  the  mode  coupling  effects  on  the  dispersion  rela- 
tion were  examined.  This  general  theoretical  framework  allows  the  inclusion  of 
longitudinal  beam  compression  effects  upon  cyclotron  waves,  and  the  inclusion  of 
transverse  gyroeffects  upon  plasma  waves.  An  analysis  was  completed  of  the 
effect  of  a finite  thermal  spread  in  the  axial  beam  velocity  on  the  electron 
cyclotron  wave.  A fairly  nonrestr icti ve  sufficient  condition  was  obtained  such 
that  the  cyclotron  modes  may  propagate  undamped  by  kinetic  effects.  Finally,  we 
reviewed  some  previous  work  on  relativistic  electron  beam-plasma  heating,  whereii 
we  found  that  very  high  density  plasmas  may  be  heated  with  reasonably  efficiency 
by  high  quality,  ul trarel at i v i st ic  electron  beams  via  the  primary  beam-plasma 
streaming  interaction. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY  OF  RESULTS 
OF  RESEARCH  ACCOMPLISHED 


This  research  program  concerns  very  high  energy, 
high  power  beams  of  charged  particles  and  the  interac- 
tions which  they  may  have  with  one  another,  with  other 
species  of  charged  particles,  or  with  surrounding 
waveguide-diode  structures  via  external  or  self-consistent 
electromagnetic  fields.  The  potential  of  such  beams  for 
rapidly  transporting  energy  may  prove  to  have  a number  of 
important  applications  for  the  Air  Force.  One  ultimate 
possibility  is  the  creation  of  advanced  new  charged 

particle  weapons  for  the  Air  Force.  A more  immediate  goal, 

* 

which  has  received  significant  attention  in  recent  years, 
is  the  possibility  of  generating  ion  beam  pulses  at 
relativistic  (i.e.,  GeV)  ion  energy  by  utilizing  the 
powerful  collective  electromagnetic  fields  which  may  be 
sustained  by  relativistic  electron  beams.  This  research 
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program  is  aimed  at  supporting,  among  other  things,  the 
development  of  such  collective  ion  accelerators. 

The  research  problems  of  central  importance  to 
this  objective  focus  generally  on  obtaining  a good  under- 
standing of  the  stable  behavior  of  unneutralized  or 
partially  neutralized  electron  and  ion  beams,  the  inter- 
action of  such  electron  beams  with  ions,  and  the 
interaction  of  such  beams  with  background  plasmas.  The 
study  of  these  beam  interactions  is  naturally 
prerequisite  to  achieving  the  beam  control  required  by 
the  various  schemes  for  ion  acceleration.  The  inter- 
actions are  predominately  collisionless,  within  an 
environment  without  gross  charge  neutrality,  and  occur  via 
powerful  electromagnetic  fields  which  are  frequently 
influenced  appreciably  by  the  presence  of  both  passive 
external  waveguide  structures  and  active  charge  or  current 
sources  (e.g.,  guide  magnetic  field  coils). 

We  initiated  our  investigation  of  these  problems 
of  relativistic  electron  and  ion  beam  interactions  under 
our  previous  AFOSR  Contract  F44620-72-C-0071 , which 


■«••• 
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terminated  on  June  30,  1976 


A technical  description  of 


this  initial  research  is  given  in  reference  1.  Work  in 
this  area  then  continued  under  AFOSR  Contract  No. 
F49620-76-C-0002 , which  began  on  July  1,  1976  and  termi- 
nated on  September  30,  1978.  This  research  is  summarized 
in  this  final  report.  Some  additional  research  is  also 
included  for  completeness,  which  was  not  accomplished 
under  this  contract,  as  is  noted  in  the  remainder  of  this 
section.  We  believe  that  substantial  progress  was 
achieved  in  this  research  on  the  interactions  of  high 
energy,  charged  particle  beams,  as  described  more  fully 
below. 

A problem  of  considerable  relevance  to  schemes 
for  collective  ion  acceleration  concerns  the  interactions 
which  may  occur  between  a relativistic  electron  beam  and 


1.  W.  E.  Drummond,  J.  R.  Thompson,  et  al.,  "Turbulent 
Heating  of  Plasmas  by  High  Energy  Particle  Beams,"  Final 
Scientific  Report  for  AFOSR  Contract  No.  F44620-72-C-0071, 
August  1976.  I-ARA-76-U-80 . 


a less  dense  species  of  ions.  The  motivation  for  con- 
sidering this  problem  is  severalfold.  Of  primary  impor- 
tance to  ion  acceleration  schemes  is  the  fact  that  such 
interactions  between  the  subject  ions  and  the  background 
electron  beam  will  occur  during  the  process  of  accelera- 
tion, when  the  interactions  have  no  desired  purpose,  but 
must  be  assessed  to  determine  how  disruptive  they  are 
on  the  electron  beam  and  on  the  process  of  acceleration. 
On  the  other  hand,  interaction  between  the  electron  beam 
and  the  subject  ions  is  also  of  interest  as  a possible 
mechanism  for  both  wave  growth  and  ion  loading  into  the 
wave,  provided  that  the  growth  saturates  via  ion 
trapping.  This  possibility  has  been  considered  by 
Indykul,  et  al.  (reference  2),  whose  analysis  was 
supported  by  a one-dimensional,  single  wave  computer 
calculation  of  the  nonlinear  ion  dynamics.  Finally, 


2.  V.  P.  Indykul,  I.  P.  Panchenko,  V.  D.  Shapiro,  and 
V.  I.  Shevchenko,  JETP  Letters  20,  65  (1974). 


such  electron-ion  interactions  stimulate  the  growth  of 
unstable  eigenmodes  on  the  electron  beam,  which  may 
eventually  achieve  large  amplitudes  before  nonlinear 
effects  cause  saturation  of  the  growth.  Consequently, 
the  electron-ion  interaction  problem  is  a handy  vehicle 
for  the  study  of  large  wave,  nonlinear  effects,  which 
in  turn  are  crucial  to  the  determination  of  the  maximum 
accelerating  electric  field  amplitudes  which  may  be 
sustained  by  the  electron  beam. 

Reference  1 contains  a discussion  and  references 
to  our  early  work  on  this  problem,  including  both  work 
sponsored  under  AFOSR  Contract  F44620-72-C-0071  and  also 
work  sponsored  under  Contract  F29601-75-C-0029  on  the 
Auto-Resonant  acceleration  scheme.  One  finds  that  the 
ions  may  interact  with  either  of  two  "negative  energy" 
electron  modes:  longitudinal  electron  plasma  oscillations 
or  electron  cyclotron  interactions.  The  former  inter- 
actions, in  general,  are  much  better  understood  from  the 
vast  amount  of  prior  research  conducted  across  the  country 
on  beam-plasma  interaction  via  plasma  waves.  We  have 
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concentrated  our  research  on  the  latter  cyclotron 

interactions,  which  will  tend  to  be  the  faster  growing 

when  the  electron  beam  is  sufficiently  energetic  and 

intense  that  (ft  /y  oj  )1//j<  1.  [Here  ft  = eB  /y  me  is 
e e pe  e z e 

the  relativistic  cyclotron  frequency,  w = ( 4ir  Ng  e2/yem)  1//z 

is  the  relativistic  plasma  frequency,  and  (y  -l)mc2  is  the 

e 

electron  kinetic  energy  per  particle.] 

The  work  on  these  ion-electron  cyclotron  inter- 
actions which  we  did  prior  to  the  July  1,  1976  beginning 
of  this  contract  was  primarily  concerned  with  the  linear 
theory  of  the  interactions  in  waveguide  geometry.  One 
important  change  from  the  usual  charge-neutral  linear 
theory  is  that  the  ion  response  is  strongly  influenced  by 
the  radial  ion  bouncing  which  occurs  in  the  self  electric 

field  of  the  electron  beam,  at  frequency  = (2y  m/M)  l/2u) 

Bi  e pe 

In  fact,  when  the  ion  density  is  low  enough  that  the 

growth  rate  of  the  interaction  falls  below  w , then  this 

bouncing  will  be  crucial  to  the  interaction,  and  the  wave 

phase  velocity  will  be  proportional  to  u)  and  independent 

B J- 

of  the  ion  density. 
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Research  on  the  nonlinear  phase  of  these  ion- 


electron  cyclotron  interactions  was  begun  under  Contract 
F29601-75-C-0046  on  the  Auto-Resonant  Accelerator 
experiment,  and  was  continued  in  the  early  months  of  this 
AFOSR  Contract  F49620-76-C-0002 . An  internal  ARA  report 
on  these  nonlinear  effects  was  compiled  in  September,  1976 
and  is  included  as  Appendix  A,  herein,  in  its  entirety, 
for  the  benefit  of  the  reader.  However,  it  should  be 
noted  that  the  material  in  Sections  III,  V,  and  a portion 
of  VII  of  Appendix  A was  sponsored  under  Contract  F29601- 
75-C-0046,  and  has  been  previously  published  in  reference  3. 
The  rest  of  Appendix  A describes  the  AFOSR-sponsored 
nonlinear  work.  This  material  in  Appendix  A was  also 
presented  in  Paper  2T3  of  the  18th  Annual  Meeting  of  the 
Division  of  Plasma  Physics  of  the  American  Physical  Society 
in  November,  1976. 


3.  W.  E.  Drummond,  et.  al.  "Conceptual  Design  of  an 
Auto-Resonant  Accelerator  Experiment,"  Final  Report  for 
Contract  F29601-75-C-0046  , May  1976.  I-ARA-76-U-44 . 
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The  chief  results  of  this  research  on  the  nonlinear 
phase  of  ion-electron  cyclotron  interactions  were  the 
identification  of  a number  of  mechanisms  by  which  the  wave 
growth  could  be  saturated.  Some  were  electron  effects, 
such  as  a nonlinear  shift  in  the  frequency  of  the  cyclotron 
wave,  and  some  were  ion  effects,  such  as  an  evolution 
of  the  ion  trajectories  in  perpendicular  (xx,  vx)  phase 
space  to  a point  where  the  ions  ceased  to  gain  energy. 

As  shown  in  Eq.  (75)  of  Appendix  A,  the  negative  energy 
cyclotron  mode  can  only  grow  so  long  as  the  ions  also  gain 
energy.  As  summarized  in  Eq.  (78)  , the  nonlinear  effects 
are  important  for  short  wavelength  modes  at  lower  wave 
amplitudes  than  for  long  wavelength  modes,  and  hence  the 
longer  wavelength  (e.g.,  small  kx)  modes  should  achieve 
the  largest  amplitudes  if  their  growth  is  otherwise  allowed. 
One  of  the  most  important  results  of  this  analysis  (which 
is  also  confirmed  by  the  two-dimensional  computer  simula- 
tions which  were  performed)  is  that  for  a wide  range  of 
parameters,  nonlinear  effects  associated  with  the  transverse 
electron  or  ion  dynamics  in  the  cyclotron  wave  enter  to 
saturate  the  wave  amplitude  before  it  becomes  large  enough 
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for  longitudinal  ion  trapping.  Hence,  the  intuition 


gained  from  the  study  of  plasma  wave  interactions  does 
not  extend  to  cyclotron  wave  interactions,  where  trans- 
verse effects  tend  to  dominate.  The  neglect  of  such 
transverse  effects  in  the  Indykul  analysis  (reference  2) 
therefore  appears  to  be  a serious  deficiency.  Our  own 
calculations  support  the  general  conclusion  that  ion- 
electron  cyclotron  interactions  tend  to  saturate  at 
wave  amplitudes  too  low  for  longitudinal  ion  trapping, 
and  before  there  is  catastrophic  disruption  of  the 
electron  beam.  Therefore,  it  appears  that  such  inter- 
actions are  (i)  probably  unsuitable  as  a wave  growth/ion 
loading  mechanism,  but  (ii)  unlikely  to  disrupt  the 
acceleration  process  itself.  We  believe  that  this 
research  on  the  interactions  between  a sparse  ion  species 
and  the  cyclotron  modes  of  an  electron  beam,  which  was 
sponsored  about  50%  under  these  AFOSR  Contracts  F44620- 
72-C-0071  and  F49620-76-C-0002 , represents  a significant 
research  contribution  in  this  field. 


A second  problem  attacked  under  this  AFOSR  contract 


was  a stability  analysis  of  a two-dimensional  sheet  beam 
of  electrons,  in  a magnetic  field  in  a conducting  wave- 
guide. This  calculation  is  included  as  Appendix  B herein. 
It  is  common  practice  in  the  theoretical  analysis  of 
problems  concerning  unneutralized  electron  beams  within 
waveguides,  to  employ  both  cylindrical  and  cartesian  (slab) 
geometry.  For  certain  analyses,  the  trigonometric 
eigenfunctions  which  apply  for  slab  geometry  afford  a 
measure  of  convenience  compared  to  cylindrical  Bessel 
eigenfunctions,  and  yet  there  is  usually  little  difference 
in  the  underlying  physical  interactions.  Nevertheless, 
there  is  a constant  need  for  vigilance  and  comparison,  to 
insure  that  the  slab  analyses  overlook  no  important  effects 
which  might  be  rooted  in  cylindrical  geometry.  To  this 
end,  the  slab  geometry  stability  calculation  is  important, 
and  it  may  be  compared  with  the  previous  calculation  we 
have  published  for  cylindrical  geometry  (reference  4) . 


4.  H.  V.  Wong,  M.  L.  Sloan,  J.  R.  Thompson,  and 
A.  T.  Drobot,  Physics  of  Fluids  1^6,  902  (1973)  . 
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Naturally,  it  is  crucial  to  understand  the  conditions 
under  which  the  electron  beam  itself  will  be  stable, 
before  addressing  more  complex  questions  of  beam 
interaction. 

For  this  analysis,  we  also  imposed  a restriction 
to  two  dimensional  perturbations,  assuming  that  one 
coordinate  is  ignorable  for  both  equilibrium  and  pertur- 
bed quantities.  Although  this  limits  the  generality  of 
our  result,  such  a restriction  is  appropriate  to  two- 
dimensional  slab  geometry  computer  codes  which  have  been 
frequently  employed  to  examine  relativistic  electron  beam 
phenomena. 

The  approach  of  the  calculation  is  to  compute  the 
second  order  quantity  / d3x  Jj *E! , which  drives  the 
electromagnetic  field  energy  in  time  according  to  Maxwell's 
equations.  Under  our  assumptions,  it  is  possible  to 
employ  the  linearized  Vlasov  equation  to  express  this 
quantity  as  a total  time  derivative,  thereby  obtaining  an 
energy  constant  of  the  motion.  This  in  turn  leads  to  an 
energy  principle  for  the  stability  of  a common  class  of 
slab  equilibria.  One  finds  that  N b/y  B must  be  bounded 

6 6 Z 
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above  to  insure  stability,  where  N is  the  beam  density 

e 

and  b is  the  half-width  of  the  conducting  waveguide. 
This  scaling  agrees  with  the  previous  cylindrical 
result. 

Another  task  we  have  undertaken  for  this  AFOSR 
contract  is  to  determine  accurately  the  limiting  currents 
for  magnetized  injection  of  an  unneutralized  relativistic 
electron  beam  into  a cylindrical  drift  tube.  This 
problem  has  been  extant  for  years,  but  so  far  only  the 
nonrelativistic  solution  (reference  5)  and  an  approxi- 
mate interpolation  formula  (reference  6)  for  the  rela- 
tivistic solution  have  been  published.  Our  work  is 
included  as  Appendix  C herein,  in  the  form  of  an  article 
appearing  in  the  November,  1978  issue  of  The  Physics 
of  Fluids.  A portion  of  this  work  was  also  supported 


5.  L.  P.  Smith  and  P.  L.  Hartman,  Journal  of  Applied 
Physics  11,  220  (1940) . 

6.  L.  S.  Bogdankevich  and  A.  A.  Rukhadze,  Soviet  Physics 
Uspekhi,  14,  163  (1971). 
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under  Contract  DASG60-76-C-0045  by  the  Ballistic  Missile 
Defense  Advanced  Technology  Center. 

This  problem  is  to  calculate  the  steady  state, 
propagating  beam  equilibrium  for  a laminar,  unneutral- 
ized, cylindrical  relativistic  electron  beam,  which  is 
initially  solid  and  monoenergetic  at  the  anode  plane. 

The  beam,  of  radius  a and  initial  kinetic  energy 

(y  - l)mc2  is  injected  along  a very  strong  axial 
o 

magnetic  guide  field  into  a coaxial  drift  tube  of 

radius  b > a,  and  length  L >>  b.  The  tube  end-effects 

and  beam-front  transients  are  neglected.  As  the  beam 

electrons  propagate  into  the  drift  tube,  the  central 

electrons  are  axially  slowed  by  space  charge  effects 

more  than  the  edge  electrons,  so  that  the  beam  density 

n(r)  becomes  peaked  on  axis  while  the  relativistic 

factor  y(r)  becomes  hollowed  on  axis.  For  given 

values  of  y and  b/a,  the  injected  current  is  a 
o 

peaked  function  of  y * Y(r  3 0) , as  illustrated  in  Fig.  1 

c 


of  Appendix  C.  It  is  desired  to  know  the  value  of  the 
peak,  or  limiting  current  IL  which  may  be  injected  into 
the  drift  tube,  as  well  as  the  corresponding  factor  y 

C Lj 

of  the  center  electrons,  as  a function  of  the  parameters 
Y0  and  b/a. 

Previously  published  relativistic  analyses  have 
essentially  ignored  the  radial  dependence  of  the  axial 
density  compression  which  occurs,  and  have  considered 
only  the  special  cases  of  needle  beams  (In  b/a  > > 1)  or 
ultrarelativistic  beams  (y  > > 1) . The  interpolation 
formula  (from  reference  6)  which  is  widely  cited  is 


me  

e 1+2  In  b/a 


and  the  corresponding  value  of  yc  is  ycL  = Yq1^3*  The 
more  exact  analysis  which  we  have  performed  indicates 
that  this  interpolation  formula  underestimates  the  true 
limiting  current  by  over  20  percent  in  some  cases. 
Furthermore,  the  corresponding  expression  for  the  axial 
electron  energy  overestimates  the  true  "limiting  energy" 
by  more  than  an  order  of  magnitude  in  some  cases.  For 
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example,  in  the  limit  of  a fat  (i.e.,  3 + 2 In  b/a  < < 

Y02/3),  ultrarelativistic  electron  beam,  we  find 

ycL  - Jl  (in  Y0  + 1 + 2 in  |)  < < T,,17’. 

Since  this  limiting  current  problem  is  of  such 

broad  relevance  to  various  electron  beam  applications, 

we  feel  that  our  theoretical  analysis  under  this  AFOSR 

contract,  supplemented  by  numerical  calculations  and 

ultrarelativistic  theory  under  contract  DASG60-76-C-0045 , 

will  be  considered  an  important  contribution  to  this 

field  of  research.  For  example,  schemes  have  recently 

been  proposed  by  Miller,  et.  al.  at  AFWL , Kirtland  AFB , 

New  Mexico  for  the  collective  acceleration  of  ions  which 

rely  on  the  controlled  motion  of  a virtual  cathode  (where 

Y =1  and  some  center  electrons  are  reflected)  created 
c 

in  an  intense  relativistic  electron  beam.  Since  virtual 

cathode  formation  occurs  due  to  excess  space  charge  slowing 

soon  after  I > I , and  since  the  incremental  distance  to 
L 

formation  of  the  virtual  cathode  is  dependent  on  y , it 

CL 

is  desired  to  know  both  I and  y accurately  in  designing 

L CL 

the  motion  of  the  virtual  cathode.  We  believe  that  our 
calculations  fill  this  need,  and  we  have  in  fact  been 
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informed  by  Miller  that  good  agreement  has  been  observed 
between  our  predictions  for  the  limiting  current  thresh- 
old and  relativistic  particle  simulations  of  virtual 
cathode  formation.  Our  work  on  this  problem  appears  to 
complement  similar  work  at  AFWL  which  has  been  based  on  a 
different  approach,  and  has  focused  more  on  the  case  of 
hollow  electron  beams. 

A fourth  problem  which  we  have  pursued  under  this 
AFOSR  contract  is  the  radial  structure  of  linear  beam 
eigenmodes  in  waveguide  geometry.  We  have  begun  with  a 
theoretical  examination  of  the  influence  which  a radial 
dependence  in  the  relativistic  factor  y(r)  has  on  the 
structure  of  the  relativistic  electron  beam  cyclotron 
mode,  and  this  calculation  is  included  as  Appendix  D herein. 
Such  a radial  dependence  in  y(r)  may  be  created  by  the 
above  mentioned  space  charge  slowing  which  occurs 
non-uniformly  when  a beam  is  injected  into  a waveguide. 
Knowledge  of  the  radial  eigenmode  structure  is  important 
to  almost  any  application  which  would  exploit  the  cyclotron 
waves,  including  particularly  the  ion  acceleration  and 
microwave  generation  schemes  based  on  these  waves. 


It  is  found  in  Appendix  D that  a variational 
principle  can  be  formulated  which  allows  an  accurate 
determination  of  the  eigenvalues  of  kxa,  where  'a'  is 
the  radius  of  a beam  bounded  by  a conducting  wall. 

These  eigenvalues  are  found  to  increase  significantly 
as  y(r)  varies  from  being  uniform  (as  at  the  anode)  to 
having  the  asymptotic  radial  dependence  approached  far 
down  the  waveguide.  For  instance,  the  value  of  kxa  for 
the  lowest  (no  nodes)  cylindrically  symmetric  eigenmode 
increases  from  2.95  to  5.77.  Accompanying  this  increase 
in  kxa  is  an  evolution  in  the  field  profiles  toward  more 
peakedness  at  the  beam  edge.  This  is  particularly  true 
of  the  radial  electric  fields  and  the  perturbed  perpen- 
dicular electron  velocities.  However,  the  axial  electric 
fields  of  the  wave  do  not  have  a significantly  altered 
radial  structure  due  to  the  y(r)  dependence. 

One  reason  that  this  calculation  is  important,  is 
that  the  larger  the  eigenvalue  kxa,  the  closer  is  the 
cyclotron  motion  of  the  electrons  to  resonance  with  the 
oscillating  fields  of  the  corresponding  eigenmode.  Hence, 
the  y(r)  dependence  may  significantly  alter  the  threshold 
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for  the  onset  of  nonlinearities,  and  may  impose  somewhat 
more  stringent  restrictions  on  the  beam  "temperature" 
threshold  at  which  Landau  damping  effects  become 
important.  In  addition,  of  course,  for  many  applications 
the  eigenmode  structure  is  relevant  in  and  of  itself. 

Appendices  E and  F contain  material 
describing  research  on  relativistic  electron  beam-plasma 
heating,  which  was  conducted  under  our  previous  AFOSR 
Contract  F44620-72-C-0071  and  was  reviewed  under  this 
present  AFOSR  Contract  F49620-76-C-0002 . The  material  in 
Appendix  E is  some  previously  unpublished  lecture  notes 
summarizing  some  of  our  conclusions  on  relativistic 
electron  beam-plasma  heating.  These  notes  were  compiled 
for  a July  1974  presentation  by  Dr.  Drummond  at  the 
Fourth  National  School  on  Plasma  Physics,  Novosibirsk, 
which  was  mentioned  in  reference  1.  Appendix  F is  com- 
prised of  three  letters  in  which  this  beam-plasma  heating 
research  is  reviewed  and  in  which  our  conclusions  are 
compared  in  some  respects  to  those  of  Thode,  et  al.  (as 
described  in  several  references  listed  in  the  third 
letter  of  Appendix  F) . The  primary  result  is  that  we 
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seem  to  be  in  agreement  that  very  high  density  plasmas 
may  be  heated  with  reasonable  efficiency  by  high  quality 
(low  angular  scatter) , high-y  (ultrarelativistic)  electron 
beams,  via  the  primary  beam-plasma  longitudinal  streaming 
interaction. 

In  Appendix  G we  reexamine  the  problem  of  the  self- 
consistent  equilibria  which  exist  for  an  unneutralized, 
relativistic  electron  beam  injected  into  a cylindrical  wave- 
guide along  an  external  magnetic  guide  field.  One  objective 
of  this  equilibrium  analysis  is  to  determine  the  self- 
consistent  radial  profiles  of  all  equilibrium  quantities, 
including  the  electromagnetic  fields  and  the  number  density 
and  flow  velocity  of  the  electron  beam.  These  equilibrium 
properties  are  in  turn  needed  to  determine  the  linearized 
eigenmodes  (or  waves)  which  may  be  supported  upon  the  electron 
beam. 

In  order  to  self-consistently  compute  the  various 
equilibrium  quantities,  assumptions  about  the  nature  of  the 
beam  equilibria  are  required,  since  many  different  equilibria 
are  possible  in  general.  It  is  most  accurate  to  make  these 
assumptions  about  the  conditions  of  creation  of  the  beam  at 
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the  cathode,  with  conservation  equations  then  utilized  to 
complete  determination  of  the  equilibrium  quantities  in  the 
downstream  region  of  the  waveguide  where  the  beam  energy  is 
to  be  extracted  for  the  application  of  interest.  In  this 
way  the  problem  of  mating  the  beam  to  the  waveguide  is 
addressed,  and  the  accessibility  of  the  equilibrium  is 
included  in  the  calculation — which  would  not  be  the  case 
if  ad  hoc  assumptions  were  instead  made  about  the  downstream 
equilibrium. 

The  use  of  equations  for  the  conservation  of  particle 
number,  particle  energy,  and  particle  canonical  angular 
momentum — which  will  be  valid  for  all  steady  state,  0- 
symmetric  electron  beams — allows  the  downstream  equilibrium 
to  be  linked  to  the  known  cathode  conditions  while  avoiding 
the  chore  of  computing  the  actual  beam  trajectories  all  the 
way  downstream  from  the  cathode.  However,  even  with  the  aid 
of  such  conservation  equations,  the  complete  set  of  equations 
which  determine  the  beam  equilibrium  is  highly  nonlinear  in 
general  and  mathematically  tractable  only  to  numerical  compu- 
ter solution.  Fortunately,  there  exist  a number  of  small 
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parameters  which  are  exploited  in  Appendix  G to  obtain  an 
iterative  analytical  solution  for  the  equilibrium,  and  these 
parameters  remain  small  throughout  the  parameter  space  antici- 
pated for  many  electron  beam  applications.  These  same  small 
parameters  also  permit  the  simplification  of  the  linear  theory 
of  waves  supported  by  the  electron  beam  equilibrium.  These 
parameters,  and  the  underlying  physical  phenomena  which  they 
reflect  are  itemized  below. 


(1)  e, 

JL 


5 


I (kA) 

e 

B (kG)a(cm)yS 
z z 


2 


— Axial  Diamagnetism, 
Relativistic  Beam 
Frame  Precession 


Because  of  the  electromagnetic  self-fields  produced  by  the 
space  charge  and  current  of  the  unneutralized  electron  beam, 
an  axial  magnetic  guide  field  must  be  present  for  equilibrium 
and  the  beam  must  precess  in  this  field  with  an  angular 

velocity  $ fr)c  to  create  a B„B  pinch  force.  However, 

8 9 z 

this  8-precession  causes  a diamagnetic  reduction  in  the 
strength  of  the  B^  field  witin  the  electron  beam,  and  the 
relative  magnitude  of  this  reduction  is  characterized  by 
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t . Also,  in  the  beam  frame  of  reference,  8 (a)  = el/2, 

-L  6 B 1 

so  tnat  e|^2  characterizes  effects  of  relativistic  pre- 
cession in  the  beam  frame. 

-1/2 

1 

(2)  Ba(a)  = — x B Forces,  S « n Current 

9 y o — ofl 

. ^6  3 

Density  Sources,  — — 

r d0 

Convective  Derivative  Terms 

Because  the  beam  equilibrium  contains  a zero  order  9-velocity, 
the  linearized  equations  which  describe  the  waves  supported 
by  the  equilibrium  will  have  a number  of  complicating  terms 
traceable  to  this  velocity.  These  terms  reflect  such  physi- 

A 

cal  phenomena  as  perturbed  9 skin  currents  on  the  beam, 
and  they  can  permit  resonances  at  the  laboratory  frame  fre- 

A 

quencies  of  interest  which  involve  9 fluting  modes  (i.e., 

*.  * 0). 


(3)  e2 


I (kA) 

e 

0.7  B2 (kG) a2 (cm) yB 
z z 


— Centrifugal  and 
Coriolis  Forces 


Because  of  the  precession  of  the  electron  beam,  there  are 

outward  radial  centrifugal  forces  ym82/r  which  must  be 

0 
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contained  to  achieve  an  equilibrium,  and  in  the  linearized 


equations  of  motion,  there  are  both  centrifugal  and  Coriolis 
forces  which  influence  the  nature  of  the  waves.  Further- 
more e_  = 2uj2  /y 2ft2,  where  u>  is  the  relativistic  plasma 
2 pe  e pe 

frequency  and  0.^  is  the  relativistic  cyclotron  frequency. 

Consequently,  (e^/2)l/z  characterizes  the  ratio  of  the 

longitudinal  plasma  frequency  uj  /y  to  the  cyclotron 

pe 

frequency,  and  hence  characterizes  the  degree  of  coupling 
which  exists  between  these  linear  eigenmodes. 


(4)  y~2  — 8,1-0  Longitudinal  Beam 

Compression 


When  this  parameter  is  small,  the  beam  is  ultrarelativistic 

in  the  laboratory  frame  of  reference,  and  the  difference 

between  the  beam  velocity  and  the  speed  of  lignt  c 

may  be  neglected.  If  8X  is  also  small,  there  can  be  no 

<P 

longitudinal  compression  of  the  beam  within  linear  theory. 


(5) 


I (kA)  A8 
e z 

17  y8  ''  3 

z z 


ay 

Yj82 

z 


Space  Charge 
Limitations 
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These  parameters  characterize  the  relative  magnitude  of 

the  radial  variation  in  y and  6 of  the  beam,  which 

z 

arises  when  the  central  electrons  are  preferentially 
slowed  by  the  space  charge  forces  which  oppose  the  propa- 
gation of  the  beam  into  a conducting  waveguide.  If  these 
parameters  are  too  large,  it  is  even  possible  for  a virtual 
cathode  to  be  established  in  the  waveguide,  from  which 
electrons  may  be  axially  reflected. 

In  Appendix  G,  analytical  expressions  for  the  beam 
equilibrium  quantities  are  derived  through  an  iterative 
process  in  which  the  parameters  of  (1),  (2),  (3),  and  (5) 

are  treated  as  small  but  finite.  Expressions  are  obtained 
which  describe  many  of  the  features  discussed  above.  One 
interesting  result  is  a generalization  of  the  pencil  beam 
limiting  current  formula  to  the  case  in  which  the  magnetic 
guide  field  is  not  infinite: 


I.  = 


me 

e 


(Yo/3  ’ 


1) 


3/  2 


1 + 2 


P b 1 

in 7 

a 4 


el(Y 


1/3 
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This  result  suggests  that  the  limiting  current  should 
increase  somewhat  as  the  magnetic  field  is  decreased, 
which  agrees  with  reported  observations  from  numerical 
simulations  of  this  problem.  However,  so  long  as  the 
parameter  is  small,  this  effect  should  be  weak,  and 

consequently  the  infinite  magnetic  field  limiting  current 
results  of  Appendix  C should  retain  their  validity  so  long 
as  <<  1. 

The  phenomena  represented  by  the  above  five 

parameters  can  also  significantly  effect  the  linear  waves 

supported  by  an  unneutralized  electron  beam.  Two  limits 

which  have  often  been  invoked  in  prior  calculations  are 

the  infinite  magnetic  limit  (for  studying  plasma  waves) 

and  the  ultrarelativistic  limit  (for  studying  cyclotron 

waves) . In  the  infinite  magnetic  field  limit,  the 

parameters  e , 6 , and  e are  neglected,  and  the 
1 0 2 

perturbed  velocities  8 . 8 also  tend  to  zero.  Hence 

0 r 

the  perturbed  electron  motion  is  strictly  axial  in  this 

case,  and  plasma  modes  are  allowed. 

By  comparison,  in  the  ultrarelativistic  limit,  the 

parameters  y,  B , and  I are  tended  to  infinity  in 

z e 
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such  a way  that  the  ratio  of  any  two  of  these  quantities 

remains  finite.  In  this  way,  the  relativistic  cyclotron 

and  plasma  frequencies  (ft  and  u>  ) remain  finite,  but 

e pe 

longitudinal  mass  terms  and  precessional  terms  tend  to 

zero.  In  this  limit  the  parameters  e , 8 , e , and  y"2 

10  2 

all  are  neglected,  and  the  perturbed  velocity  8 also 

z 

tends  to  zero.  Hence  the  perturbed  electron  motion  is 
strictly  transverse  in  this  case,  and  cyclotron  modes  are 
allowed . 

In  either  the  ultrarelativistic  or  the  infinite 
magnetic  field  limits,  there  is  no  coupling  between  the 
cyclotron  and  plasma  modes.  Also,  in  neither  of  these 
limits  are  the  Ay  space  charge  effects  intrinsically 
small,  and  the  radial  dependence  which  they  entail  greatly 
complicates  the  mathematical  eigenvalue  problem  within 
linear  theory.  Consequently,  we  have  examined  these  space 
charge  effects  in  separate  calculations.  In  Appendix  D 
the  space  charge  effects  on  the  electron  cyclotron  wave 
are  considered  within  the  ultrarelativistic  limit,  while 
in  Appendix  I the  space  charge  effects  on  the  longitudinal 
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plasma  wave  are  considered  within  the  infinite  magnetic 
field  limit. 

In  Appendix  H,  the  linear  theory  of  collective 

waves  supported  by  an  unneutralized  relativistic  electron 

beam  is  examined  in  the  limit  that  the  parameters  e^, 

8.,  Ay/Yz  and  A8  /8  are  negligible,  but  the  parameters 
0 z z 

and  y~2  remain  finite.  This  allows  the  linear  theory 
to  be  solved  in  terms  of  relatively  simple  eigenfunctions 
(i.e.,  Bessel  functions),  while  retaining  the  effects  of 
longitudinal  compression,  centrifugal  and  Coriolis  forces, 
which  are  significant  for  many  applications.  In  this  way, 
the  dispersion  relation  for  plasma  and  cyclotron  waves  may 
be  derived  which  retains  the  coupling  effect  between  these 
waves.  Either  the  ultrarelativistic  or  the  infinite  magnetic 
limits  may  still  be  taken  if  desired,  but  at  a later  stage 
in  the  calculation.  For  both  plasma  waves  and  cyclotron 
waves,  the  coupling  to  the  other  mode  is  characterized  by 
the  e2  parameter.  For  example,  the  dispersion  relation 
for  the  negative  energy  cyclotron  wave,  subject  to  certain 
orderings  discussed  in  Appendix  H,  may  be  approximately 
given  by 
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where  k2  = k2  + k2.  Here  the  2u>2  /y2  correction  term 
l z pe 

reflects  the  centrifugal  and  Coriolis  forces,  while  the 

k2/y2k2  correction  reflects  the  longitudinal  beam 

compression.  These  longitudinal  plasma  wave  corrections 

are  in  opposite  directions,  and  upon  expanding  the  square 

root,  their  net  size  relative  to  the  cyclotron  frequency 

is  seen  to  be  about  — e„. 

4 2 

In  Appendix  I,  the  Ay  space  charge  effects  on 

the  longitudinal  plasma  wave  are  considered  in  the  infinite 

magnetic  field  limit.  Previous  calculations7  in  the  pencil 

beam  limit  y2/^  3 <<1  + 2 In  b/a  had  suggested  that  the 

normalized  plasma  wave  phase  velocity  $.  = ui/k  c could 

$ z 

only  be  made  small  (as  desired  for  ion  loading  in  a Conver- 
ging Guide  Accelerator)  by  increasing  the  electron  current 
extremely  close  to  the  limiting  current  which  can  be 


7.  R.  J.  Briggs,  Phys.  Fluids  19^,  1257  (1976). 
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propagated  without  virtual  cathode  formation.  In  this 

Appendix  I,  we  found  through  both  analytical  and  numerical 

calculations,  that  this  unfortunate  result  appears  to  also 

hold  for  more  realistic  nonpencil  beam  geometries  where 

the  radial  variation  in  the  equilibrium  quantities  Y,  3 , 

z 

is  significant.  The  plasma  wave  eigenfunctions  were 
also  examined,  and  found  to  be  somewhat  more  peaked  near 
the  beam  axis  when  the  limiting  current  conditions  were 
approached. 

Finally,  in  Appendix  J we  have  examined  the  effect 
of  a finite  thermal  spread  in  the  axial  beam  velocity  on 
the  dispersion  relation  for  the  electron  cyclotron  beam 
mode  in  the  ultrarelativistic  limit.  The  calculation  is 
performed  in  slab  geometry,  and  the  equilibrium  distribu- 
tion function  is  taken  to  be  a 6-function  in  energy  with 
a spread  in  transverse  momentum.  This  approximates  the 
thermal  spread  acquired  by  a relativistic  electron  beam 
upon  passage  through  an  anode  foil.  The  result  which  was 
found  is  that  if  the  transverse  velocity  scatter  v is 
low  enough  so  that 
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and  if  finite  Larmor  radius  effects  can  be  neglected 


then  the  electron  cyclotron  beam  modes  will  propagate 
undamped  by  kinetic  effects. 

We  believe  that  the  results  produced  under  this 
contract  have  been  substantial,  and  have  relevance  to  a 
moderately  broad  spectrum  of  problems  of  interest  to  the 
Air  Force.  We  believe  that  this  research  has  helped  to 
raise  the  level  of  understanding  of  beam  interactions  to 
a point  that  will  soon  permit  the  practical  and  efficient 
realization  of  a number  of  applications,  such  as  collective 


ion  acceleration  and  microwave  generation. 
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NONLINEAR  EFFECTS  IN  AUTO-RESONANT  ACCELERATION 


by 

W.  E.  Drummond 
J.  R.  Thompson 
H.  V.  Wong 
G.  I.  Bourianoff 

Austin  Research  Associates,  Inc. 

I.  INTRODUCTION 

In  recent  years,  a number  of  collective  ion  accelera- 
tion methods  have  been  proposed.  In  1973,  a passively  con- 
trolled method  (known  as  Auto-Resonant  acceleration)  of 
collectively  accelerating  ions  within  the  potential  troughs  of 

„ 1 r 2 

a negative  energy  mode  of  the  electron  beam  was  suggested, 
which  appears  to  afford  the  possibility  of  producing  up  to 
kiloamperes  of  ion  current,  at  relativistic  ion  energy,  with 
pulse  durations  which  might  be  an  appreciable  fraction  of 
the  electron  beam  pulse  duration.  Among  the  attractive 
features  of  this  concept  are  the  following: 

(i)  the  collective  feature  allows  fairly  large 
accelerating  electric  fields  of  order  MV/cm, 
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permitting  reasonably  short  accelerators; 
in  addition  the  collective  fields  of  the 
electron  beam  permit  substantial  ion  currents 
in  the  range  of  0.1  - 1 percent  of  the  electron 
current  to  be  produced  without  unacceptable 
performance  degradation; 

(ii)  the  use  of  a travelling  body  wave  of  the 
electron  beam,  rather  than  a localized 
potential  depression,  permits  the  acceleration 
of  a quasi-continuous  train  of  ions  throughout 
a major  portion  «of  the  electron  beam  pulse 
duration;  in  addition,  the  ion  energy  per 
particle  can  be  raised  to  a level  bounded 
theoretically  only  by  the  ion-to-electron  mass 
ratio  times  the  electron  energy  per  particle; 

(iii)  the  passive  method  of  controlling  the  accelera- 
tion of  the  wave  (and  ions)  avoids  technological 
switching  difficulties;  and 

(iv)  the  negative  energy  feature  of  the  wave  permits 
the  electron  beam  to  automatically  supply  power 
both  for  accelerating  the  ions  and  for  sustaining 
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the  electric  potential  of  the  travelling 
wave  as  it  accelerates  the  ions. 

Although  one  or  another  of  these  features  are  common  to 
various  other  schemes  of  ion  acceleration,  in  combination 
they  seem  to  be  unique  to  Auto-Resonant  acceleration. 

The  negative  energy  electron  beam  mode  which  was 
selected  for  Auto-Resonant  acceleration  was  the  lower  (slower) 
electron  cyclotron  eigenmode.  This  mode  was  selected  because 
it  alone  seemed  to  possess  the  desired  characteristics  of  a 
phase  velocity  which  could  be  varied  from  a value  much  less 
than  c (the  speed  of  light)  to  near  the  beam  velocity,  as 
well  as  a simple  method  of  passively  controlling  the  phase 
velocity:  flaring  the  guide  magnetic  field  to  spatially 

reduce  its  strength. 

Although  the  central  features  of  this  accelerator 
concept  are  embodied  in  a single  "accelerator  section" 
wherein  the  magnetic  field  strength  is  reduced  to  accelerate 
the  cyclotron  wave,  a number  of  other  sections  or  components 
will  be  necessary  to  complete  a practical  accelerator  system. 
These  components  may  be  conceptually  itemized  as  follows: 
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(1)  High-voltage  power  supply 


(2)  Electron-beam  diode 


(3)  Beam-compression/tailoring  section  (optional) 

(4)  Cyclotron  wave-growth  section 


(5)  Ion  loading  section 


(6)  Ion-wave  accelerator  section 


(7)  Accelerator  termination 


Since  these  components  are  intrinsically  interrelated, 
a program  of  theoretical  support  for  the  construction  of  such 
a system  must  encompass  each  of  these  areas.  We  describe 
herein  the  results  of  only  a small  portion  of  such  theoreti- 
cal research,  concerning  nonlinearities  in  the  cyclotron  wave 
amplitude,  which  may  affect  the  central  components  (4) , (5) , 
and  (6) . 

Naturally  much  of  the  early  calculations  of  the  Auto- 
Resonant  Accelerator  have  been  based  upon  linear  analysis  of 
the  electron  cyclotron  mode,  which  is  valid  when  the  mode 
amplitude  is  sufficiently  small.  For  example,  the  detailed 
investigation  of  various  means  of  growing  the  cyclotron 
eigenmode  has  been  conducted  within  linearized  theory,  as 


t ;>  ■ . 

_ . * : * 


has  the  analysis  of  ion  loading  and  acceleration.  However,  it 
is  of  interest  to  determine  the  mode  amplitude  level  at  which 
nonlinear  effects  first  become  important  for  a number  of 
reasons : 

(1)  the  successful  functioning  of  the  accelerator 
for  nonlinear  wave  amplitudes  has  not  yet 
been  theoretically  established  (or  disproven) ; 

(2)  nonlinear  effects  may  well  cause  saturation 

of  the  linear  growth  mechanisms,  establishing  an 
upper  bound  to  the  achievable  mode  amplitude; 

(3)  theoretical  knowledge  of  such  nonlinear  satura- 
tion amplitudes  and  mechanisms  is  needed  for 
comparison  with  computational  and  experimental 
investigations;  and 

(4)  knowledge  of  the  magnitude  and  parameter  scaling 
of  nonlinear  mode  amplitudes  is  a valuable  design 
tool,  since  this  information  bears  heavily  on  the 
required  length  of  both  the  growth  and  accelera- 
tion sections  of  the  accelerator. 

Therefore  we  have  examined  situations  where  the 
amplitude  of  the  electron  cyclotron  mode  increases,  such  as 
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for  an  electron-ion  instability  or  an  electron-waveguide 
instability,  to  determine  the  amplitude  level  at  which 
nonlinearities  become  important.  The  latter  case  may 
include  either  a resistive  liner  interaction,  where  the 
resistance  allows  the  negative  energy  cyclotron  mode  to 
grow  in  amplitude,  or  else  an  interaction  between  beam 
modes  and  modes  supported  by  some  form  of  surrounding  slow 
wave  structure.  Such  interactions  are  relevant  as  possible 
growth  mechanisms  for  the  cyclotron  mode  prior  to  ion  load- 
ing and  wave  acceleration. 

On  the  other  hand,  interaction  between  the  electron 
beam  and  the  subject  ions  is  also  of  interest  as  possible 
mechanism  for  both  wave  growth  and  ion  loading,  provided 
that  the  growth  saturates  via  ion  trapping.  This  possibility 
has  been  considered  by  Indykul,  et  al.,3  whose  analysis  was 
supported  by  a one-dimensional,  single  wave  computer  calcula- 
tion of  the  nonlinear  ion  dynamics.  In  addition,  electron-ion 
interaction  will  also  occur  during  the  process  of  accelera- 
tion, when  it  has  no  desired  purpose,  but  must  be  assessed  to 
determine  how  disruptive  the  interaction  is  on  the  electron 
beam  and  on  the  acceleration  process.  This  electron-ion 
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interaction  problem  has  been  selected  as  the  vehicle  of  our 
detailed  nonlinear  calculations,  although  some  of  the  results 
are  also  relevant  when  the  cyclotron  mode  is  excited  in  a 
different  fashion. 

The  influence  of  nonlinear  effects  may  be  broadly 
characterized  in  two  different  ways,  and  our  results  have 
been  divided  into  the  resulting  four  categories  as  shown  on 
the  chart  below: 


Electron 

Ion 

Effects 

Effects 

Violation  of  primary 
conditions  for 
linearity 

Section  III 

Section  IV 

Growth  saturation 
without  violating 
primary  linearity 
conditions 

Section  V 

Section  VI 

The  first  division  separates  nonlinearities  which  influence 
the  electron-wave  dynamics  from  those  which  influence  the 
ion -wave  dynamics . For  example,  the  Indykul  calculation 


X 


would  encompass  only  ion  nonlinearities  since  the  electrons 
were  "legislated"  to  obey  linear  theory;  likewise,  an 
electron -waveguide  instability  can  manifest  only  electron 
nonlinearities . 

The  second  division  separates  those  cases  where  it 
is  possible  to  construct  a nonlinear  theory  which  retains 
its  validity  to  the  point  of  growth  saturation,  so  that 
the  saturation  amplitudes  are  known  precisely,  from  cases 
for  which  linear  approximations  breakdown  before  a satura- 
tion mechanism  appears.  In  this  latter  circumstance, 
experience  frequently  suggests  a correlation  between  wave- 
growth  saturation  and  the  breakdown  of  linear  approxima- 
tions, which  can  often  be  substantiated  by  computer 
simulations  or  by  experiment. 

Before  describing  our  analytical  results,  we  first 
describe  briefly  our  approximations  and  the  linear  theory 
of  the  electron  cyclotron  wave. 
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II.  APPROXIMATIONS  AND  LINEAR  THEORY 


Below  are  definitions  of  the  symbols  and  notation 


which  has  been  used: 


(x,y,z) 


= coordinates  of  slab  geometry 


= unit  vector  in  the  axial  direction 

= unit  vector  in  the  radial  direction, 
normal  to  the  beam  and  waveguide 
boundaries 

= unit  vector  in  the  transverse,  ignorable 
direction 

= radius  (half-width)  of  the  conducting 
waveguide 

- speed  of  light 

•-  external  magnetic  guide  field 

= DC  self  electric  field  of  the  electron  beam 


B = DC  self  magnetic  field  of  the  electron  beam 

oy 

= DC  self  electric  potential  of  the  electron 
o . 

beam 

(E  ,E  ,E  ) [or (E  ,E  , E )]  = perturbed  components  of  the  wave 
X y 2 electric  field 

(B  , B ,8  )[or(B  ,B  ,B  ) ] = perturbed  components  of  the  wave 
x y 2 magnetic  field 


) 

w 


A- (or  A ) 


<l>(or  <Z) 


k (or  k ) 
z o 

K*  (or  k ) 

i 


uu  (or  (U  ) 
o 


- e , m 


u = p c 
o o 


v (x) 
oy 


= certain  linear  combinations  of  E , E , 

B , B X 7 

x y 

= perturbed  electric  potential  function 

= axial  component  of  the  wave  vector 

= radial  component  of  the  wave  vector 

= k2  + k2 
A o 

= frequency  of  the  electron  cyclotron  wave 

= electron  charge,  mass 

= mean  electron  number  density 

= electron  beam  radius 

= DC  electron  axial  flow  velocity 
, 2 -1/2 

= (1  - 3 ) = DC  electron  y-f actor 

o 

= mean  transverse  electron  drift  for  radial 
equilibrium 

l/  2 

= (4tt  n es/y  m)  = relativistic  electron 
plasma  frequency 

= (e  B /y  me)  = relativistic  electron 
cyclotron  frequency 

* 0uUp/k2c2  = width  of  cyclotron  resonance 
= perturbed  electron  density 
= perturbed  radial  electron  coordinate 
= perturbed  radial  electron  velocity 
= perturbed  axial  electron  velocity 
■ perturbed  electron  y-factor 
-10- 


e,M 


N.  (X,t) 
1 — 

n 

Ol 


a . 
1 


x . 

Ol 


V 

ox 


= ion  charge,  mass 
= total  ion  number  density 
= mean  ion  number  density 
= radius  of  ion  distribution 

= radial  coordinate  of  equilibrium  ion  motion 

= radial  velocity  of  equilibrium  ion  motion 

= (4rr  n e2/M)1/2  = radial  bounce  frequency  of  an  ion 
the  s'elf  electric  field  E of  the  electrons 


,JJBz'  ^Bx 


'D 


PI 

n. 

i 


ox 

= k ,ki[e$(0)/M]1/2  = axial,  radial  bounce  frequency 
of  ion  in  the  wave  electric  field 

= (4tt  n . e2/M)1/2  = ion  plasma  frequency 
oi 

= (e  B /Me)  = ion  cyclotron  frequency 
oz 


n^(or  n^)  = perturbed  ion  density 


x . 

l 


z 

J (or  J) 
p (or  o) 

r 

Ao) 


= perturbed  radial  ion  coordinate 

= perturbed  radial  ion  velocity 

= perturbed  axial  ion  velocity 

= perturbed  total  current  density 

= perturbed  total  charge  density 

=>  linear  growth  rate  of  electron  cyclotron  wave 

= nonlinear  frequency  shift  of  the  electron 
cyclotron  wave 


Perturbed  quantities  with  a tilde  (e.g.,  E^)  are  representa- 


tions in  the  usual  configuration  space,  while  those  without 
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the  tilde  (e.g.. 


E ) are  representations  of  the  Fourier 

coefficients  for  the  basis  functions  exp(-  iuut  + ik  z)  . 

z 

Also,  the  quantities  uu,  k , K+  are  the  frequency/ 

z 

wavenumber  components  for  a general  wave,  while  u , k , 

o o 

kx  are  the  particular  values  existing  for  a pure  electron 
cyclotron  mode. 

As  a basis  for  the  nonlinear  calculations  to  be 

described  later,  we  first  consider  a slab  model  equilibrium 

in  which  a cold  relativistic  electron  beam  of  uniform  density 

n propagates  with  uniform  velocity  u -*  c along  a uniform 
o o 

guide  field  3 in  the  z -direction,  and  interacts  with  a 
oz 

sparse  group  of  ions  whose  axial  velocity  is  zero.  For 
ready  reference,  a number  of  the  approximations  or  restric- 
tions which  underlie  this  calculation  are  listed  below. 


az  at 


0 for  equilibrium 


(1) 


ay 


0 for  both  equilibrium  and  waves 


(2) 


E (x)  = £ ( - x) 
z z 


(3) 
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r 


(n  ,/n  ) 
01  o 


1/  2 

8 (y  m/M)  (Q/'ju  ) 
o p 


< 1 


2uu 

B 

n 


1/2 

2 (Y  m/M)  ( uu  /n)  < 1 
o p 


n. 

— = (Y  m/M)  (Q/uu  ) < 1 

o, B  ° P 


(fi  ./n  ) 
oi  o 


8 (M/y  m) 
o 


1/  2 


(n3u)  A4c4) 

p 


1/  s 


< l 


r (axial  ion  response) 

r 


0.08  (y  m/M)1/2(nB 
o 

(n  ./n  ) 
oi  o 


( 


( 


( 


a4  /<d  c4  ) 

E 


(It  is  desired  that  these  inequalities  be  well  satisfied.) 


Both  the  electron  beam  and  the  ions  are  assumed  to  be 


T 


The  electrons  are  assumed  ultrarelativistic , according  to 
Eqs . (8)  - (10),  which  allow  neglect  of  the  transverse 
drift  velocity  v^  as  well  as  the  diamagnetic  corrections 
to  the  guide  magnetic  field.  The  electron  equilibrium  is 
assumed  uniform,  as  indicated  in  Eq.  (7). 

If  we  subject  this  equilibrium  to  small  amplitude 
perturbations,  electron  beam -cyclotron  modes  propagate  along 
the  beam.  These  cyclotron  modes  will  be  largely  decoupled 
from  the  longitudinal  plasma  modes  due  to  the  high-v 
assumptions,  and  they  are  the  subject  of  this  linear  analysis. 
Let  us  define  the  new  variables 

u 

A*  = E ± i E ± i — (B  ±iS) 
x y c x y 

From  Maxwell's  equations,  we  obtain  the  following  equation 
for  A-  . 


~±  / 5 

A = 4rr  ( — ± 
\ dx 


u 


£ )(Jx  * 1 


(19) 
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where  the  boundary  conditions  on  A~  at  x = ± a are: 


A+  « A" 


aA*  + 5A~ 
dx  9x 


u 


(20) 


If  we  restrict  the  analysis  to  perturbation  with  no 

variations  in  y(a/ay  = 0),  and  consider  perturbations  of 

the  form  exp(-  iuut  + i k z)  , then  from  Eq . (17)  and  Eq. 

z 

(19) , we  obtain: 

( s?  + k j a;  - ° (21) 

where 


tC 


±2 


u>3  (<u  - k u 

E 5_ 

c2  (tv  - k u 
z o 


± n) 


(22) 


Assumption  (8)  allows  neglect  of  the  perturbations  in  the 

axial  electron  velocity  because  of  the  large  longitudinal 

mass.  Likewise,  the  electron  high-v  assumptions  imply 

o 

that  p -(u  /c2)  J -»  0 and  may  be  neglected, 
e o ze 
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Solving  Eq.  (21)  subject  to  assumptions  (4),  (5), 

and  the  boundary  conditions  (20) , yields  the  dispersion 

relation  for  the  electron  beam  cyclotron  mode  (uj  = , 

o 

k = k , K+  = k,)  symmetric  in  E : 
z o J-  z 


!i>  = k u 

o o o 


n k2  c2 

(k2  C2  + «,*) 

p 


k2  = k2  + k2  (23) 

o -i. 

tan  k±a  = - k±a 

Here  assumption  (4)  for  nonrelativistic  wave  speed  tends  to 
be  justified  for  these  interactions  with  ions,  and  permits 
simplification  of  the  cyclotron  eigenfunctions  by  decoupling 
the  light  modes.  Assumption  (5)  implies  strong  decoupling 
of  the  positive  and  negative  energy  cyclotron  modes,  which 
also  simplifies  the  cyclotron  eigenfunctions. 

For  the  lowest  modes,  K*a  = k^a  « 2,  4.9,  8,  . . . 

The  eigenfunctions  are 


A+  (x)  = A sin  K + x 
o o 


A“  (x)  = A 
o o 


sin  '<*  a 
sin  K-a 


sin  K"x 


(24) 


as  A — sin 
o a 


K a 


i V 

Letting  A^  = |A^ | e , it  may  be  shown  that  the  perturbed 
electric  potential  function  is 


k c2 |A  ' 
1 o 


so  that 


$(x,t) 

Cm 

^ a2 

e5>(0) 

k . e A 
x o 

Y me2 

Y m n2 

o 

o 

aarkxx  cos  (k_z  - uuQt  + co)  (25) 


(26) 


and  the  corresponding  perturbed  electric  fields  are 


E 

z 


(27) 


As  indicated  by  assumption  (3)  we  are  considering  only  those 
waves  symmetric  in  (x) . Qualitatively  similar  results 
may  easily  be  derived  for  the  antisymmetric  waves. 
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If  a small  fraction  of  ions  with  mean  density 

n , is  introduced  into  the  electron  beam,  electron-ion 
01 

streaming  instabilities  occur  due  to  coupling  of  electron 
beam-cyclotron  modes  to  the  ion  motion.  Assumption  (11) 
insures  that  the  equilibrium  is  not  appreciably  affected 
by  the  ions , and  we  obtain  for  A : 


/ S2  ±2 \ ± 9 

( + K ) A - 4ne  ni(x) 


(28) 


where  n^ (x)  is  the  perturbed  ion  density. 

Equation  (28)  may  be  solved  by  first  order  perturba- 
tion expansion  about  the  cyclotron  mode.  Let 


A (x)  = A + A + 
o 1 


IJJ  - W + (JU_ 

o 1 


k = k + k 
z o 1 

where 

AT  < A“  , uj.  < ui  , k.  < k , 
1 O 1 O 1 o 


and  A~  is  given  by  Eq.  (24) . 

The  most  stringent  condition  for  the  convergence  of 
this  perturbation  expansion  is  obtained  from  the  perturbation 
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/ %r.  ** 


of  the  quantity  K+2,  defined  in  Eq.  (22).  This  convergence 

condition  eventually  leads  to  Eq . (15)  that  the  growth  rate 

F be  less  than  6'u  3 Q'ju2/k2c2  . 

P 

By  subtracting  Eq.  (21)  from  Eq.  (28),  and  forming  a 
Rayleigh -Ritz  integral  to  avoid  solving  for  the  first  order 
eigenfunctions  A^,  we  obtain 


+a 


I 


-a 


dx  A (x)  uu  - — K + k T r—  K A (x) 
o I 1 Sou  1 9k  I o 

o o 


+a 

w 4rre  f dx  Ta+  (x)  + A-  (x)  1 — n.  (x) 
« L o o J Sx  i 

-a 


(29) 


To  complete  the  dispersion  relation,  we  need  the  perturbed 
ion  density  n^(x). 

Since  the  ions  will  be  non -relativistic  in  the 
laboratory  frame  of  the  calculation,  the  dominant  forces  on 
the  ions  are  those  due  to  the  electric  fields,  and  in  the 
limit  of  Eq.  (5),  K+2  a k2  > k2,  and  the  E^-component 
is  the  largest.  Thus,  we  obtain  for  the  ion  dynamics: 
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dx 

dt 


dv 


= v 


(30) 


x 

dt 


*bx  + it  Vx'z>t) 


where  the  term  - u^x  is  the  radial  restoring  force  due  to 

the  self -electric  field  E of  the  electron  beam,  and  the 

ox 

perturbed  electric  field  is  of  Eq.  (27) . 

For  sufficiently  high  ion  density  and  for  low  magnetic 

field  strength,  the  growth  rate  of  these  ion-electron  cyclotron 

interactions  can  exceed  the  perpendicular  bounce  frequency 

a;  and  these  equilibrium  oscillations  may  be  ignored. 

However,  for  electron  beam  and  magnetic  field  parameters 

typical  of  state  of  the  art  experimental  facilities,  this  may 

require  ion  loading  corresponding  to  kiloamperes  of  output 

ion  current.  Hence,  we  shall  not  discuss  this  regime  herein, 

but  will  address  the  more  likely  regime  of  lower  ion  density, 

where  the  growth  rate  is  below  rU  as  in  Eq.  (12),  and  the 

B 

perpendicular  ion  oscillations  are  very  important.  However, 
it  will  also  generally  be  true  that  2 ajg  < n as  in  Eq.  (13), 
so  that  the  cyclotron  wave  phase  velocity  is  non-relativistic 
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and  the  wave  number  is  not  much  perturbed  from  the  value 

Q/3  c.  In  addition,  Q . < -jj  as  in  Eq . (14),  so  that  the 
o i B 

ion  coupling  to  the  guide  field  is  truly  negligible. 

From  the  linearized  Vlasov  equation  for  the  ions,  we 
determine  the  perturbed  ion  density,  assuming  it  to  be 
dominated  by  the  transverse  ion  dynamics . After  substitution 
in  the  right  hand  side  of  Eq.  (29),  we  obtain  the  dispersion 
relation 


k u M [ (kacaj/(kac2  + <oa)l  n + Auu_  from  the  zero  order 
o o P B 

dispersion  relation  (23). 

An  equilibrium  ion  distribution  for  which  all  ions 
have  the  same  energy  and  the  same  radial  turning  points  is 
given  by 

9 (a)  - 6<a3  -.?)  (32) 

O TT  'U  1 


where  n is  the  average  ion  density,  and  a.  £ a is  the 
oi  1 

amplitude  of  the  equilibrium  ion  trajectory.  The  equilibrium 

ion  density  function  will  be  peaked  near  the  turning  points 

x = ± a . . 

i 

Prom  Eq . (31), 


'JU 


1 


k,  u » 
1 o 


2 i >M3  . <i)2  J (k  a ) J'(k  a ) 
pi  p j 1 i l 1 i 

(u  uu  Q k a (1  + cos2  k a) 

1 B x- 


■ - — J.Oc.a  ) J,' 

uu^  l 1 r l L r 


where 


2 l uu 


'JU 


£i E 


(JU 


B 


Q k^a 


(1  + cos  k a) 
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For  maximum  growth. 


0,  and 


F = Im 


= A 


/ 


j;(k  a 


i>] 


1/  2 


(33) 


An  alternate  choice  for  the  equilibrium  ion  distri- 
bution, for  which  the  ion  density  is  parabolic  in  x, 
peaked  about  the  beam  center,  is 

2 n . 

9o(a)  = 7~ft.  (34) 


where  © is  the  unit  step  function.  For  this  distribution, 
the  maximum  growth  rate  is  given  by 


r 


Vi 


-n  1/2 


(35) 


which  is  fairly  close  to  the  prior  result. 

The  assumption  that  the  interaction  of  the  ions  is 
dominated  by  the  transverse  ion  response  rather  than  the 
axial  response  is  quantified  in  Eq.  (16)  for  the  case 
a^  •*  a — which  insures  that  the  axial  response  yield  a smaller 
growth  rate.  This  inequality  bounds  the  magnetic  field 
strength  from  above  and  the  ion  density  from  below,  contrary 
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to  most  of  the  other  restrictions,  and  tends  to  require 

smaller  magnetic  fields  than  Eq.  (14) . 

Conditions  (5)  and  (13)  will  also  insure  that 

Q m a2  < 4cs  so  that  the  influence  of  v B forces 
B i ox  oy 

on  the  axial  ion  motion  will  not  affect  the  interaction 
appreciably. 

In  addition  to  the  ion-electron  cyclotron  inter- 
actions which  we  are  considering,  there  may  be  longitudinal 
ion-electron  plasma  interactions  and  transverse  Weibel-like 

interactions  which  can  occur  for  k -»  0 . However, 

z 

assumption  (10)  that  the  longitudinal  electron  plasma 
frequency  be  below  the  electron  cyclotron  frequency,  assures 
that  the  plasma  and  cyclotron  modes  are  well  decoupled,  and 
also  that  the  plasma  oscillations  occur  in  the  "magnetized" 
regime.  Consequently  the  large  transverse  k -component 
indicated  in  Eq.  (5)  is  sufficient  for  electromagnetic 


, effects  to  stabilize  the  fastest  longitudinal  plasma 

10 

oscillations  for  which  k 8 c =-  oj  / y . There  will  remain 

z o p o 

much  slower,  longer  wavelength  plasma  oscillations  with 


k 8 c »•  a)  . / but  the  growth  rate  of  these  modes  will  be 

7 c\  m J 


r 


r 

i 

i 

i 

i 

r 


below  the  cyclotron  mode  growth  rate  under  the  given 
assumptions . 

For  magnetic  field  strengths  large  enough  to  obey 
Eq.  (13)  well,  the  transverse  Weibel  interactions  also 
appear  to  be  stabilized  by  the  oun-  oscillations  of  the 
ions . Only  for  smaller  magnetic  fields  would  these 
interactions  be  unstable  with  significant  growth  rates. 

Therefore,  for  the  equilibrium  conditions  and 
assumptions  which  have  been  set  forth  herein,  the  fastest 
ion-electron  interaction  occurs  between  modes  of  transverse 
ion  oscillation  in  the  beam  electric  field  and  electron 
cyclotron  modes.  To  determine  the  wave  amplitude  at  which 
nonlinearities  become  significant,  we  shall  focus  on  the 

electron  cyclotron  dynamics  and  the  ion  uj  -oscillations . 

B 
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III.  BREAKDOWN  OF  LINEAR  ELECTRON  BEHAVIOR 


In  this  section,  a brief  survey  is  presented  of  a 
number  of  nonlinear  electron  effects,  in  which  one  or 
another  of  the  primary  conditions  for  linearity  is  violated. 
Rough  estimates  are  given  of  the  magnitude  and  scaling  of 
the  mode  amplitude  level  at  which  these  effects  become 
important.  Many  of  the  effects  may  be  seen  to  be  related 
to  one  another.  Although  care  has  been  taken  to  preserve 
the  numerical  coefficients,  which  are  presented  here  for  the 
case  of  slab  geometry,  the  very  nature  of  these  formulae  is 
merely  to  indicate  the  level  at  which  linear  approximations 
lose  their  validity,  and  this  occurs  gradually  rather  than 
sharply  as  a function  of  amplitude.  The  correlation  between 
the  breakdown  of  linear  approximations  and  wave  growth 
saturation  must  remain  inferential,  pending  the  accumulation 
of  additional  theoretical  or  computational  evidence. 

In  addition  to  the  cyclotron  eigenmode  structure, 
which  was  presented  in  Eqs . (24)  through  (27),  one  may  derive 
from  linear  theory  the  expression  for  the  perturbed  radial 
coordinate  of  an  electron: 
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X4 


i-.  m 


r-  - 


x(x) 


k2c2 


sin  k x 
x 


(36) 


e $ (0) 


which  is  the  magnitude  of  the  perturbation  amplitude — with 
the  2,  t phase  information  supressed.  Likewise  $(0)  is 
the  magnitude  of  the  perturbed  electric  potential  amplitude 
at  the  beam  axis  x = 0.  This  expression  characterizes  the 
electron  particle  dynamics  in  the  wave,  and  shows  that 
electrons  near  the  beam  edge  are  the  most  strongly  perturbed 
for  the  broadest  perpendicular  eigenmode  (for  which  k^  =“  2/a)  . 

The  various  conditions  wherein  nonlinear  effects 
become  significant  are  listed  below  for  ready  comparison 
with  one  another;  a discussion  of  the  different  nonlinear 
effects  follows . 


JL  » kac3  e j(0) 

n w2  y me2 

o p o 
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o o 
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(40) 
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ox 


k . a 21/3 
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Condition  (37)  prescribes  the  wave  amplitude  level  at 
which  the  electron  density  fluctuations  due  to  the  wave 
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become  comparable  with  the  mean  electron  density.  These 
density  fluctuations  will  be  strongest  on  axis  for  this 
equilibrium  where  nQ  and  yq  are  uniform  in  x,  and 
at  this  wave  amplitude  the  electron  density  will  vary  from 
near  zero  to  around  twice  the  ambient  density. 

Condition  (37)  also  indicates  that  these  density 

nonlinearities  will  be  come  important  at  an  amplitude  level 

smaller  by  u)2/k2c2  from  that  for  which  longitudinal 
P 

trapping  of  the  electrons  by  the  wave  potential  would  occur. 

Condition  (38)  indicates  three  other  nonlinear  effects 
which  enter  at  the  same  wave  amplitude  level.  The  first  of 
these  is  the  condition  that  the  radial  displacement  of  an 
electron  by  the  wave  be  as  large  as  the  scale  length  k"1 
characterizing  the  perpendicular  eigenmode  structure.  This 
is  judged  to  be  more  accurate  than  the  criterion  for  dis- 
placements as  large  as  the  beam  radius,  which  requires 
higher  wave  amplitudes.  Equation  (36)  indicates  that  the 
displacements  will  be  largest  at  about  0.8  times  the  beam 
radius.  One  physical  consequence  of  radial  displacements 
as  large  as  k”1  =**  a/2  is  that  wall  collisions  might  occur 
unless  a sufficiently  wide  vacuum  gap  was  available.  The 
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second  condition  in  Eq.  (38)  is  that  finite  electron 
gyroradius  effects  become  important  because  of  the  perpen- 
dicular kinetic  energy  excited  by  the  wave.  The  third 
condition  is  that  the  radial  motion  of  the  electrons 
becomes  nonlaminar  by  virtue  of  radial  "crossings"  of 
electron  trajectories.  This  may  be  easily  seen  to  occur 
when  |3x/dx|  exceeds  1,  which  occurs  first  for  particles 
on  axis.  These  last  two  conditions  in  Eq.  (38)  are 
symptomatic  of  the  breakdown  of  the  fluid  nature  of  the 
wave  due  to  nonlinear  turbulence. 

Condition  (39)  prescribes  when  the  perpendicular 
wave  electric  field  exceeds  the  self  radial  field  of  the 
beam,  which  tends  to  occur  first  on  axis.  This  requires 
a wave  amplitude  larger  than  that  for  density  nonlinearity 
by  only  k3/k3.  Condition  (39)  also  reveals  the  connection 
between  these  various  criteria  of  nonlinearity  and  the 

measurable  ratios  such  as  E /E  (a)  which  are  con- 

x max  ox 

veniently  extracted  from  computer  simulations  as  a diagnostic. 

One  physical  consequence  of  the  nonlinearity  suggested  by 

Eq.  (39)  is  the  possibility  of  electron  field  emission  at 

the  conductor  walls.  The  self -field  E is  typically  much 

ox 
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larger  than  the  threshold  amplitude  required  for  field 
emission,  but  is  in  the  wrong  direction.  One  might  think 
that  if  the  wave  fields  became  larger  than  the  self-field 
and  had  the  opposite  radial  sense,  then  field  emission 
might  result.  However,  this  is  a deception  of  the  linear 
theory  which  is  breaking  down,  since  Poisson's  equation 
shows  that  even  for  nonlinear  waves,  the  electrostatic 
portion  of  the  electric  field  must  always  point  toward  a 
region  of  net  negative  charge.  Hence,  the  only  possibility 
of  field  emission  is  due  to  the  inductive  or  electro- 
magnetic portion  of  the  electric  field,  which  is  small 
for  nonrelativistic  phase  velocity,  or  else  to  the  presence 
of  the  ions  whose  density  might  locally  exceed  the  nonlinear ly 
reduced  electron  density.  Although  this  criterion  for  field 
emission  is  not  quantified  herein,  we  anticipate  that  the 
required  ion  density  will  be  a significant  fraction  of  the 
electron  density,  since  the  electron  wave  fields  are 
geometrically  reduced  below  their  peak  value  at  the  walls, 
and  the  peak  value  does  not  exceed  the  self -field  until 
(n/nQ)  2 k2a2/kia  according  to  Eq.  (39) . 
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Condition  (40)  for  the  radial  electron  oscillations 


in  the  wave  to  become  relativistic  will  be  stronger  than 
condition  (38)  only  if  k^  exceeds  k^ . It  is  possible 
that  this  might  be  modestly  true  in  the  growth  section  and 
at  the  beginning  of  the  accelerator  section,  but  toward 
the  end  of  the  accelerator  section  k will  be  below  k 


z 


i 


and  g « 1. 
x 

The  significance  of  the  0^  ^ 1 nonlinearity  is 

that  it  implies  a significant  fluctuation  y in  the 

electron  y-factor  as  shown  in  condition  (41) . This  occurs 

because  the  electrons  are  perturbed  by  the  wave  to  move 

with  or  against  the  radial  self-field  of  the  beam.  At  the 

beginning  of  the  accelerator  where  I is  largest, 

E (a)/B  is  smallest  and  may  be  somewhat  less  than  1. 
ox  oz 

The  net  product  of  these  terms  yields  an  expression,  also 

shown  in  Eq.  (41) , which  indicates  that  y/y  is  generally 

o 

somewhat  less  than  n/n  . 

o 

t 

However,  condition  (42)  indicates  that  the  variation 

y/y  is  magnified  by  kac2/(jua  to  produce  a variation  in 
o p 

the  fluid  electron  cyclotron  resonance  |uu-k  c@  +n|  which 

Z Z 

may  be  comparable  to  its  ambient  value  (i.e.,  6uu)  . This 


-34- 


nonlinear  spreading  of  the  cyclotron  resonance  occurs  due 


to  the  y-dependence  of  Qj  it  is  sometimes  called  "resonance 
broadening."  It  is  also  symptomatic  of  the  breakdown  of 
fluid  theory  and  the  conversion  of  the  cyclotron  mode  to 
a kinetic  regime  due  to  nonlinear  turbulence.  Condition 
(42)  indicates  that  this  nonlinear  resonance  broadening 
effect  is  felt  at  wave  amplitudes  somewhat  below  those  for 
density  nonlinearity,  and  hence  is  the  most  sensitive 
nonlinear  effect  yet  discussed.  Experience  would  suggest 
that  this  resonance  broadening  would  reduce  the  wave  growth 
rate  when  the  wave  amplitude  reaches  the  level  given  by 
Eq.  (42). 

Condition  (43)  gives  the  corresponding  resonance 
broadening  which  occurs  due  to  changes  in  the  axial  electron 
velocity.  The  first  terms  of  Eq.  (43),  which  are  linear  in 
the  wave  amplitude  n/ nQ,  are  relatively  unimportant  as 
usual  for  Vq  » 1.  However,  the  last  term  of  Eq.  (43), 
which  is  quadratic  in  the  wave  amplitude,  may  be  seen  to  be 
more  (or  less)  important  than  the  y/y  resonance  broaden- 
ing, depending  on  whether  kaca/2ka  uua  a2  is  greater  (or 

z p 

less)  than  1.  This  quadratic  term  in  the  wave  amplitudes 
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reflects  the  perturbation  in  0 due  to  linear  perturbations 

z 

in  0 , according  to 

1 ~, 

0^  — 0 L + linear  terms 

+ other  quadratic  terms  (44) 

which  follows  immediately  from  a perturbation  expansion  of  y(0). 

The  square  brackets  in  conditions  (41)  through  (43) 

are  hybrid  representations  in  which  the  dominant  first  term 

reflects  the  influence  of  transverse  electromagnetic  fields 

on  the  electrons,  such  as  0 E , while  the  smaller  second 

X ox 

term  shows  the  influence  of  the  perturbed  axial  electric 

field  E . 

z 

The  mode  amplitude  at  which  these  nonlinear  effects 
appear  is  reduced  for  higher  values  of  k2 , and  consequently 
the  higher -k  modes  should  saturate  in  growth  at  lower 
amplitudes  than  the  lower-k  modes.  Hence,  when  exciting 
or  growing  the  cyclotron  wave,  care  should  be  taken  (as  in 
designing  computer  simulations)  that  the  lowest  k mode  is 
not  discriminated  against  by  any  frequency  or  wavenumber 
selection  which  may  be  imposed. 
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Although  this  list  of  electron  nonlinearities  includes 
L ^ a variety  of  effects,  it  is  by  no  means  complete.  As  just 

one  further  example,  the  possibility  of  a nonlinear  three - 
wave  decay  of  the  electron  cyclotron  wave  was  examined  and 
j found  to  be  disallowed  by  the  coupling  selection  rules. 


li 

li 

il 


IV.  BREAKDOWN  OF  LINEAR  ION  BEHAVIOR 


As  in  the  previous  section,  we  here  present  estimates 
of  the  perturbation  amplitude  level  at  which  certain  primary 
conditions  of  ion  linearity  cease  to  be  valid.  As  before, 
these  estimates  are  rough,  and  the  inference  of  wave  growth 
saturation  at  these  wave  amplitude  levels  will  again  be 
tentative  pending  additional  theoretical  or  computational 
evidence . 

The  radial  ion  motion  is  described  by  Eq.  (30) ; the 
zero  order  ion  motion  is  seen  to  be  simple  harmonic  oscilla- 
tion at  frequency  <ju  , with  uniform  amplitude  x . = a. 

B or  i 

for  the  distribution  of  Eq.  (32).  Any  ion  distribution 

symmetric  in  x will  be  carried  by  the  zero  order  particle 

motion  into  an  equivalent  state  during  one -half  bounce 

period,  or  ttAu  . If  this  time  also  equals  an  integral 
B 

number  of  wave  periods,  then  the  perturbed  wave  quantities 
will  have  also  returned  to  their  initial  state.  Conse- 
quently Am  are  the  resonant  frequencies  of  the  linear 
B 

ion-electron  cyclotron  interaction,  where  f * 2,  4,  6,  ...  . 
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Equation  (30)  may  also  be  solved  for  the  perturbed  ion  motion, 
and  the  result  for  the  simplest  1 = 2 interaction  is 


— - x . 
tjjg  i max 


. 2 2 
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JL 


e $ (0) 


The  axial  ion  motion  is  described  by  dv  /dt  =- 

z 

(e/M)  E^x,  z,  t)  since  conditions  (5)  and  (13)  render  the 
Vox  Boy  forces  unimportant.  This  equation  may  be  solved 
to  yield  the  perturbed  axial  ion  motion 


V e $ (0) 


n e $ (0) 


Based  upon  an  examination  of  the  ion  dynamics , the 
following  conditions  for  primary  ion  nonlinearity  may  be 
derived. 
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where  according  to  Eqs . (33)  or  (35),  the  magnitude  of  the 
1 = 2,  a^  = a growth  rate  is  about 
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(51) 


The  first  two  conditions  (47)  and  (48)  relate  to  the 
transverse  ion  dynamics,  while  the  latter  conditions  (49) 
and  (50)  relate  to  the  axial  dynamics.  In  addition,  condi- 
tions (48)  and  (50)  are  sensitive  to  the  ion  density  through 
r,  while  conditions  (47)  and  (49)  are  not. 
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Condition  (47)  prescribes  the  wave  amplitude  level 
at  which  the  radial  ion  bouncing  near  the  axis  is  influenced 
by  the  wave  fields  as  much  as  by  the  DC  self -fie Id  of  the 
electron  beam.  Quantitatively  this  condition  is  exactly 
the  same  as  condition  (39) . Since  the  ion-electron  inter- 
action is  strongly  affected  by  the  resonance  with  DC  ion 
bouncing,  one  suspects  that  this  wave  amplitude  would  com- 
pletely alter  the  coherence  of  the  interaction. 

Condition  (48)  is  that  the  radial  displacement  of 
an  ion  be  as  large  as  the  scale  length  k^1  of  the  perpen- 
dicular eigenmode.  At  this  amplitude  level,  which  may  be 
fairly  low  when  r is  small,  one  world  expect  the  growth 
rate  to  decline  as  it  might  when  the  ions  had  a transverse 
" temperature . " 

Condition  (49)  is  that  some  ions  be  longitudinally 

trapped  by  the  wave.  This  is  the  growth  saturation 

mechanism  appealed  to  in  the  one -dimensional  calculation 

3 

of  Xndykul,  et  al.  However,  it  may  be  seen  here  to  require 
'tirlv  large  wave  amplitudes,  particularly  when  k^  is 
• •"'Ms  is  because  at  low  ion  density,  the  frequency 


«sl 


of  the  unstable  wave  cease  to  depend  on 


n . through  w . and  are  instead  determined  by  the 
01  pi  J 

frequency  <0g  of  radial  bouncing.  Hence  the  wave  speed 

does  not  approach  closer  and  closer  to  the  ion  speed  as 

H . is  reduced.  Condition  (49)  will  occur  before  or  after 
01 

condition  (47)  according  to  the  ratio  k /k^ . 

Condition  (50)  prescribes  the  wave  amplitude  level 

at  which  the  perturbed  axial  velocity  of  the  ions  broadens 

the  resonance  at  uu  « 2oj  by  an  amount  comparable  to  its 

B 

ambient  value  (i.e.,  P)  . This  effect  will  be  felt  at  a 

wave  amplitude  lower  than  that  required  for  longitudinal  ion 

trapping  by  P/m  . However,  it  will  not  occur  before  condi- 

B 

tion  (48)  so  long  as  the  usual  ordering  k < k is  obeyed. 

z i 

Compared  to  some  of  the  electron  nonlinearities, 
these  ion  conditions  are  not  quite  as  sensitive  to  high 
k-values,  but  two  of  the  ion  conditions  are  sensitive  to 
low  ion  density  and  hence  can  occur  in  some  cases  well 


before  the  electron  nonlinearities. 


V.  GROWTH  SATURATION  VIA  NONLINEAR 
ELECTRON  FREQUENCY  SHIFT 

In  this  section,  we  discuss  the  effects  of 
nonlinearities  in  the  electron  response  which  produce  a 
nonlinear  shift  in  the  wave  frequency  of  the  electron 
beam -cyclotron  mode.  We  expect  nonlinear  frequency  shifts 
will  influence  the  time  evolution  of  the  beam  cyclotron 
mode  since  growth  is  determined  by  resonant  coupling  to  the 
bounce  motion  of  the  ions . 

It  is  assumed  throughout  that  a fluid  description 
of  the  electron  beam  is  adequate.  In  evaluating  the  non- 
linear frequency  shift,  the  procedure  is  to  calculate  the 
nonlinear  generation  of  the  zeroth  and  second  harmonic 
components  of  the  fundamental  mode.  The  nonlinear  inter- 
action of  the  fundamental  with  these  harmonics  then  leads 
to  corrections  to  the  linear  dispersion  relation. 

The  details  of  this  calculation  will  not  be  presented 
here.  The  result  obtained  for  the  frequency  shift  Acju  is 


2 


n k2c2 

A(U  « - Q 2 

8 Wp 


e k,  A 
-i  o 

my  ft2 

o 


k2c2 

[4.!5  + 3.33  -Sr-  ] 
P 


s - A 


(52) 


where  the  fundamental  is  taken  to  be  the  lowest  beam 
cyclotron  mode  with  k^a  « 2.  In  deriving  Eq.  (52),  which  is 
insensitive  to  the  method  of  exciting  the  electron  cyclotron 
mode,  we  have  assumed  the  primary  linearity  condition 
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(53) 


as  well  as  the  frequency  ordering 

1 Aoo  | < Scu  = Q <ju2/k2c2  . (54) 

P 

Condition  (53)  is  just  the  reverse  of  the  nonlinear  condition 
(37) , while  condition  (54)  is  a time-ordering  criterion 
invoked  to  allow  calculation  of  the  "slow"  shift  in  the 
wave  frequency. 

If  we  denote  the  complex  mode  amplitude  Aq  by 
A a a elcp,  and  insert  this  nonlinear  electron  frequency 
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shift  into  the  equations  describing  the  growth  of  the 
ion-electron  cyclotron  instability,  then  we  obtain  the 
following  equation  for  the  time  evolution  of  the  mode 
amplitude 


£ a eitp  = Ta  / a eicp  dt'  + i 4 0s  e1® 


where  r is  the  linear  growth  rate  of 


Eqs.  (33)  or  (3  5), 


and  A is  the  frequency  shift  coefficient  of  Eq.  (52) 
Decomposing  this  equation  into  its  real  and 
imaginary  parts,  we  obtain: 


2 fS  + a|l|  , 34(ja  dg 

at  dt  dt3  dt 


The  first  integral  of  Eq.  (57)  is 


^ A Cl4  = C s constant 


Asserting  the  initial  condition  dco/dt  = 0 when  a = 0,  we 
obtain  C = 0.  Finally,  substituting  dcp/dt  from  Eq.  (58) 
into  Eq.  (56)  yields 


= p2  Q _ J_  .2  qS  _ 

dt2  1 a 16  A a 


Eq.  (59)  is  the  equation  of  motion  of  a "particle"  in  a 
"potential  well"  V^a),  illustrated  in  Figure  1 below,  it 
describes  the  time  evolution  of  a from  given  initial 
conditions  <2  (0)  and  dQ/dt  (0)  . 


Figure  1.  "Potential  Well"  V, (a) 


For  small  a[i.e.,  a < (T/a)1/2]/  0 contains  a 

component  growing  exponentially  with  the  linear  growth  rate 
P.  However,  when  A2  0*716  =“  r2  , da/dt  = 0,  and  growth 
ceases . Thereafter  the  amplitude  decreases  and  then  may 
continue  to  oscillate  in  time,  depending  upon  the  precise 
choice  of  initial  conditions.  The  maximum  amplitude 
attained  is  given  by 

1/  2 


which  corresponds  to 


E 

x max 

k2  0 
x max 

32 

1/  2 

c 

rill 

1/  3 

E (a) 
ox 

k2  4rr  n ea 
o o 

(4.15+3.32 

&u; 

P 

. n J 

P 


where  the  growth  rate  has  been  inserted  for  the  interaction 

with  a.  = a “ 2/k  and  i = 2 = id  /iu  • It  may  be  observed 
i -i-  o B 

that  this  maximum  wave  amplitude  is  lower  than  the  lowest 
amplitude  for  electron  resonance  broadening  [given  in  Eqs . 
(42)  and  (43)]  by  about  (2  to  6)  vF/ 6uu  , which  tends  to  be 
small  due  to  assumptions  (15)  and  (54) . 
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VI.  GROWTH  SATURATION  VIA  PHASE  MIXING  OF 
NONLINEAR  ION  TRAJECTORIES 

In  this  section,  we  discuss  the  effects  of  nonlinearities 
in  the  ion  response  on  the  time  evolution  of  the  unstable 
electron  beam -cyclotron  mode. 

We  again  treat  the  ions  as  a small  perturbation  and 
evaluate  the  nonlinear  trajectories  of  the  ions  in  the  wave- 
fields  of  the  electron  beam -cyclotron  mode.  Substituting 

k® 

E (x,  z,  t)  = - rr  a sin  k x cos  (k  z + cp  - uu  t) 

X K Z O 

o 


in  Eq.  (30),  and  defining  the  new  variables  a/  %,  t. 


x - a sin  (uj  t + ?) 


Vx  ” a 005  ^tUBt  + ^ 
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a dt ' 


we  obtain  from  Eq.  (3  0) 
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where 


d£  _ 3H 
dr 


■ Jt(^)  OJS  (q  + Cp) 


Ocx  a)2 


q = kQZ  + £? 


uu  - lm  , and  i is  even . 

O B 


In  Eq.  (63),  we  have  neglected  terms  which  are  oscillatory 
in  t.  These  equations  describe  the  slow  variation  of  a 
and  ? due  to  the  resonant  interaction  of  the  ions  with  the 
oscillatory  electric  field. 

Assuming  that  cp  is  independent  of  t [see  Eq.  (67)], 

H is  a constant  of  the  motion.  Since  p and  q must  vary 

such  that  J («/2p  ) and  m (q  + cp)  never  pass  through  zero, 
z 

the  variations  of  p and  q are  bounded. 
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Figure  2 shows  a sketch  of  the  ion  trajectories  in 

the  polar  coordinate  phase  space  (V 2p , q)  due  to  resonant 

coupling  to  modes  with  frequency  uo  = 2uo  . The  points 

o B 

given  by  j'(k,a  ) = 0,  an  (q  + co)  = 0,  (i.e.  k,a  « 

2 -L  S s 1 S 

3,  6.7,.  . . ; qs  + cp  = 0,  rr)  are  stable  stationary  points 
The  dashed  circles  have  radii  k,x  given  by  jr  (k  x ) = 0 
i.e.,  k x « 5.1,  8.4,  ....  The  ion  trajectories 

cl 

indicated  by  the  solid  lines,  are  restricted  to  regions 
bounded  by  the  dashed  circles.  The  initial  equilibrium  ion 
distribution  of  Eq.  (32)  would  be  a circle  of  radius 
V 2p  = k±a^  • while  that  of  Eq.  (34)  would  be  the  interior 
of  the  same  circle.  The  conducting  boundary  is  a circle  of 
radius  V"2 p = k±a. 

If  all  the  ions  initially  have  energy  such  that 
k,x.  < k,a,  5.1: 

(a)  The  ion  trajectories  eventually  intersect  the 
wall  due  to  interaction  with  the  lowest  beam 

« 

cyclotron  mode  kxa  = 2; 

(b)  A small  fraction  of  the  ion  trajectories 
intersects  the  wall  when  k±a  = 4.9; 
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(c)  The  ion  trajectories  do  not  intersect  the  wall 
when  k^a  =8,  or  for  higher  k±  modes. 

From  Eq.  (63),  we  obtain  for  p: 


1/  2 

3?  = 4 - QSJ 


(64) 


where  Q = constant.  This  equation  can  be  solved  by 
quadrature  with  the  initial  conditions 
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p = p0  * 


, q = q = k z + IS 
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at  "time"  t - 0. 

Since  g = g [a  (x,  v ) ] is  the  initial  ion  distri- 
o 00  x 

bution  of  the  ions,  the  ion  density  at  time  t is  given  by 


+00 

N . (x,  z,  t)  = f dv  g fv  (x,  z,  v , t)  , x (x,  z,  v , t)  1 
1 J x oL  ox  x o x J 

—00 

(65) 


where  (x  , v ) is  the  point  in  phase  space  at  time  t = 0 
o ox 

which  arrives  at  (x,  v ) at  time  t. 
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The  contribution  of  the  perturbed  ion  density  to  the 


right  hand  side  of  Eq.  (29)  is  given  by 
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where 


P = P(P  / q / t)  / q = q(p  / Q , r) 
o o o o 


are  the  solutions  of  Eq.  (63) . J dt/T  is  a time  averaging 

o 

operator,  where  1/r  > T > l/a)  • 

o 

Equating  the  real  and  imaginary  parts  of  Eq.  (29) , 
using  Eq.  (66) , we  obtain 


“t?  - 0 


(67) 
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where 


k (ju 

G ( t)  = - f da  a g (a  J f dq  J (k  a)  sin  (q  + op) 

2 n . J oooo  J o 4 x 


Equation  (68)  is  the  equation  of  motion  of  a "particle" 
moving  in  a "potential  well"  (t) , and  it  describes  the 
time  evolution  of  Q,  which  is  proportional  to  dT/dt. 

For  small  values  of  r < 1, 


J (k  a)  sin  (q  + cp)  » J (k  a ) sin  (q  + cp) 
i x l x o o 


J . <k  a ) J ' (k  a ) 
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and  it  follows  from  Eqs . (68)  through  (70)  that  the  amplitude 
a grows  exponentially  with  the  linear  growth  rate. 
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It  is  apparent  from  Eq.  (63)  that  nonlinear  modifica- 


tions to  the  exponential  linear  growth  of  the  instability 
occur  when  &q  > tt/2  , that  is 
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(72a) 
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for  k^a  ■ k^a^  = 2 * i . These  conditions  (71)  and  (72)  will 

occur  for  wave  amplitudes  smaller  than  those  given  in  the 

primary  electron  nonlinearity  conditions  (37)  through  (39) 

by  the  factor  r/iui  < 1-  However,  the  wave  amplitudes  of 

B 

(71)  and  (72)  are  only  higher  by  a numerical  factor  of  about 
tt/2  from  the  nonlinear  ion  condition  (48)  , since  this  con- 
dition addresses  the  same  physical  effects, 
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A detailed  description  of  the  nonlinear  time  evolution 
of  <2  requires  the  evaluation  of  (t) , which  is  analytically 
arduous.  However,  a representation  of  the  nonlinear  time 
evolution  of  Q which  has  immediate  physical  significance 
and  can  be  related  to  the  nonlinear  ion  trajectories  of 
Figure  2,  can  be  obtained  as  follows. 

We  begin  by  verifying  that  the  linear  dispersion 
relation  for  unstable  growth,  as  well  as  the  generalization 
to  include  the  effects  of  nonlinear  ion  dynamics,  is 
consistent  with  energy  conservation: 


The  electron  contribution  is 


The  ion  contribution  is 
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where  oj  = , and  i is  an  even  integer, 

o B 

The  field  contribution  is 
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V is  the  volume  of  integration  in  the  above  expressions. 
Thus  from  Eq.  (73) 
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which  is  identical  to  Eq.  (68) . 

Eq.  (74)  can  be  integrated  in  time,  using  the  version 
(74a)  , to  yield 
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where  0(0)  - 0 , and 


(75) 
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Consequently,  the  increase  in  the  wave  amplitude  is  directly 
related  to  the  transverse  energy  gained  by  the  ions  in  the 
interaction . 

We  expect  the  "coarse  gained"  ion  distribution 
function  will  evolve  to  a state  where  it  is  flattened  (i.e., 
made  uniform  by  phase  mixing)  in  the  region  of  phase  space 
(k^a,  q)  traversed  by  the  ions.  For  example,  in  the  case 
where  the  ion  trajectories  do  not  intersect  the  wall,  and 
with  reference  to  Figure  2,  if  the  initial  ion  distribution 
is  that  of  Eq.  (34)  , where 
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with  k a < 5.1,  <j,  » 2m,  k,x  =■  5.1,  k a * 8,  then  the 
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predicted  final  state  distribution  would  be 
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and  substitution  in  Eq.  (75)  yields  an  estimate  of  the 
final  wave  amplitude 
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where  r*  is  the  linear  growth  rate  of  Eq.  (3  5)  . This  differs 

from  the  prior  estimate  of  Eq.  (72a)  only  by  1/2tt  times  the 

square  bracket  term,  which  may  vary  from  0 to  1.03  depending 

on  the  size  of  a./x  . Here  the  final  state  amplitude 

i a 

approaches  zero  as  a.  -»  x simply  because  it  is  not  possible 

X 21 


-61- 


for  the  ions  to  increase  their  net  energy  by  much  if  they 

are  initially  filled  approximately  to  the  boundary  of 

allowed  phase  space  trajectories.  As  remarked  earlier, 

the  RHS  of  Eq.  (77)  contains  the  small  factor  T/2ijj  which 

3 

may  produce  a result  below  the  various  other  conditions  of 

nonlinearity.  If  not,  then  such  other  nonlinearities  might 

prevent  the  approach  to  this  phase-mixed  final  state.  In 

addition,  it  must  be  recognized  that  this  "final"  state  may 

also  be  unstable,  since  perturbations  of  the  outermost  ions 

may  place  them  on  new  trajectories  outside  the  "boundary 

at  x . | 

a • 
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VII.  RESULTS  OF  TWO-DIMENSIONAL  SLAB 
COMPUTER  CALCULATIONS 


During  the  past  year  or  so,  we  have  been  running  and 
improving  large  computer  simulation  codes,  which  are  fully 
electromagnetic  and  fully  relativistic,  and  operate  in  a 
two-dimensional  slab  configuration  space.  We  have  applied 
these  codes  to  the  problems  of  ion-electron  streaming 
interactions  and  resistive  growth  interactions  in  an  external 
magnetic  field,  and  have  observed  the  growth  of  the  electron 
beam  cyclotron  mode.  For  these  runs  in  which  an  instability 
is  computed  from  the  inrtial  noise  levels,  periodic  boundary 
conditions  are  applied  in  the  axial  dimension  and  conducting 
boundary  conditions  are  applied  in  the  radial  dimension. 

The  gross  charge  is  non-neutral,  since  even  for  the  cases 
of  ion-electron  interaction,  the  ion  density  is  much  below 
that  of  the  relativistic  electron  beam.  For  a number  of 
reasons  described  below,  we  have  not  yet  been  able  to  put 
the  previously  described  theory  of  nonlinearities  to  a 
good  test,  and  so  we  do  not  claim  that  this  theory  is 
substantiated  by  the  computer  results.  Nevertheless,  we 
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wish  to  describe  some  of  the  significant  general  features  of 


the  computer  runs,  in  which  behavior  characteristic  of  the 
transverse  dynamical  effects  which  we  have  analyzed  is  in 
fact  manifested. 

In  order  to  achieve  a good  simulation  test  of  the 
nonlinear  theory,  it  is  necessary  to  design  a run  in  which 
most  of  the  essential  conditions  (1)  to  (16)  are  obeyed, 
together  with  a number  of  practical  constraints  for  com- 


pleting a numerically  stable  and  accurate  simulation  within 
an  affordable  period  of  time.  Although  we  have  accumulated 
about  a half  dozen  large  and  lengthy  (i.e.  many  thousand 
time  steps  with  many  thousand  particles)  simulation  runs 
during  the  past  year,  most  of  these  failed  to  satisfy  some 
of  the  important  conditions  (1)  to  (16) . This  was  partly 
because  the  initial  computer  runs  were  designed  before  the 
desired  theoretical  constraints  were  known.  In  addition, 
the  instability  domain  has  proved  to  be  rich  with  variations, 
so  that  in  several  runs  in  which  certain  inequalities  were 
only  weakly  obeyed,  several  competing  effects  are  believed 
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Y ~ 3),  then  some  longitudinal  plasma  mode  effects  will 
o 

be  coupled  into  the  electron  cyclotron  mode  behavior. 

It  has  been  necessary  to  devote  significant  attention 
during  the  simulation  effort  to  curing  certain  purely 
numerical  instabilities  that  are  unique  to  relativistic, 
electromagnetic  codes . Central  among  these  is  the 
numerical  Cherenkov  instability,  in  which  the  finite 
differencing  of  Maxwell's  equations  allows  light  modes 
with  wavelengths  on  the  order  of  the  cell  size  to  have 
phase  velocities  below  c,  so  that  they  may  interact 
unstably  with  electron  beam  modes . These  may  be  cured 
by  filtering  out  the  shortest  wavelength  modes,  and 
arranging  for  the  wavelengths  of  physical  interest  to  be 
many  cell  lengths. 

The  simulation  runs  which  have  been  completed, 
although  violating  some  of  the  theoretical  constraints, 
have  shown  the  following  common  behavior.  When  the  electrons 
and  ions  are  released  in  approximate  conditions  of  equili- 
brium, modes  at  the  predicted  axial  wavenumber  k ~ fi/c, 

z 

and  with  low  transverse  wavenumber  components , grow  from 
noise  for  a few  e -folds  at  a rate  generally  faster  than  the 
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theoretical  estimate  of  the  linear  growth  rate.  However, 

this  initial  growth  phase  terminates  with  the  mode 

amplitude  well  within  the  bounds  of  linear  theory  [e.g., 

~ few  x 10~3],  and  is  usually  followed  by  a phase  of 

growth  at  a rate  from  1 to  0.2  times  the  estimate  of  the 

linear  growth  rate,  with  the  growth  rate  declining  in  time 

in  several  cases.  In  some  runs  the  wave  growth  saturated, 

with  the  amplitude  in  one  case  ranging  up  to  <$/<£  ~ 7%, 

o 

while  in  other  cases  the  run  was  terminated  with  <$/<!>  a 

o 

few  percent  but  still  slowly  growing.  For  the  runs  which 
showed  amplitude  saturation,  rough  theoretical  estimates 


suggested  that  nonlinearities  such  as  electron  cyclotron 
resonance  broadening  should  have  been  felt  at  amplitudes 
about  a factor  of  two  below  the  observed  saturation.  For 
the  runs  which  did  not  saturate,  rough  theoretical  estimates 
indicate  that  nonlinearities  should  not  have  been  felt  for 
an  additional  1 to  2 e-foldings. 

Rather  than  attach  any  significance  to  the  degree  of 
quantitative  agreement  between  these  runs  and  the  only 
partly  relevant  theoretical  estimates  of  growth  rates  and 
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saturation  amplitudes,  we  would  emphasize  the  following 
qualitative  results: 


(1)  The  desired  electron  cyclotron  mode  was  excited 
to  significant  amplitudes,  and  the  expected 
transverse  electron  dynamics  were  visible  in 
the  phase  space  diagnostics — as  shown  in  Figure 
3 (where  y is  the  axial  coordinate) . 

(2)  The  radial  ion  bouncing  in  the  electric  fields 
of  the  beam  and  wave  clearly  dominated  the  ion 
dynamics,  and  no  significant  axial  motion  of 
the  ions  was  detectable.  In  no  case  were  the 
ions  "trapped"  longitudinally  by  the  wave. 

(3)  In  those  cases  where  the  wave  approached  the 
"nonlinear"  amplitude  level,  symptoms  of  the 
breakdown  of  the  electron  fluid  nature  were 
evident — such  as  trajectory  crossings  and  other 
significant  phase  space  distortions.  However, 
in  no  case  were  the  distortions  so  severe  that 
the  electrons  were  "ripped  apart,"  or  reflected 
in  the  axial  direction,  etc. 


We  are  presently  involved  in  constructing  runs  in 
which  a very  quiet  equilibrium  state  is  established,  with 
a number  of  diagnostics  added  to  verify  all  of  the  linear 
wave  characteristics.  Then  we  will  proceed  to  a run  in 
which  this  nonlinear  theory  may  be  tested. 
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VIII.  CONCLUSIONS 


We  have  calculated  herein  certain  nonlinear  effects 
which  appear  when  the  amplitude  of  the  cyclotron  eigenmode 
of  a relativistic  electron  beam  is  increased  sufficiently. 

We  have  used  as  the  vehicle  of  this  calculation,  the 
streaming  interaction  which  arises  between  relativistic 
electrons  and  a sparse  group  of  ions  which  are  axially  at 
rest  in  the  reference  frame  of  the  calculation.  Both 
electron  and  ion  nonlinear  effects  were  isolated,  and 
the  electron  nonlinearities  are  generally  relevant  for 
other  situations  without  ions,  where  the  cyclotron  wave  is 
excited  via  some  other  beam-waveguide  type  of  interaction. 

These  calculations  were  completed  subject  to  a 
number  of  parameter  restrictions,  including  those  of  con- 
ditions (1)  to  (16) , which  are  generally  appropriate  to  an 
Auto-Resonant  ion  acceleration  device.  Under  these  restric- 
tions, the  negative  energy  cyclotron  wave  is  well  decoupled 
from  the  longitudinal  plasma  modes,  light  modes,  and  the 
positive  energy  cyclotron  mode,  and  the  interaction  of  the 
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ions  with  the  cyclotron  mode  is  faster  than  any  other 

ion-electron  interactions.  The  ion  density  was  assumed 

small  enough  that  the  ion  coupling  to  the  cyclotron  wave 

could  be  calculated  via  perturbation  theory,  and  small 

enough  that  the  growth  rate  of  the  interaction  was  below 

both  the  perpendicular  ion  bounce  frequency  .jj  , and 

B 

the  width  of  the  electron  cyclotron  resonance 

6uu  - Q uu2  A2c3  . 

P 

The  estimates  of  nonlinear  effects  were  evaluated 

in  detail  for  the  specific  case  where  a single  unstable 

eigenmode  dominated  the  wave  spectrum,  with  the  axial 

wave  electric  field  symmetric  in  the  transverse  coordinate 

x,  with  the  lowest  perpendicular  wave  number  k =“  2/a  > k , 

■1  z 

and  with  the  lowest  resonant  wave  frequency  iju  = 2uj  < k c. 

o B z 

A hybrid  expression  which  summarizes  the  chief 
results  of  the  various  calculations  of  the  wave  amplitude 
required  for  strong  nonlinearity  is 


♦ (x 
max 

♦ (x  * 
o 


>1’  V^’^JleBssr  X. 
[V.’+D  ' 3 ”B 

L 4 uj3  J 

P 


lesser 


in  slab  geometry. 


where  [e  $ (a)/y  me3]  =-  cu  2a3/2c3 
o o p 

The  coefficients  D^,  , D 4 are  of  order  unity  and 

will  generally  vary  somewhat  for  different  eigenmode 

characteristics  k^ , &,  etc.  The  term  represents 

the  primary  electron  nonlinearities,  of  which  the  most 

sensitive  is  cyclotron  resonance  broadening.  The  D 2 

term  represents  the  nonlinear  electron  frequency  shift 

effect,  which  is  important  when  T/ 6uu  is  very  small. 

The  term  represents  the  ion  nonlinearities  due  to 

transverse  phase  space  trajectories  or  to  axial  resonance 

broadening  (depending  on  whether  k or  k dominates 

i z 

the  k"2  coefficient),  and  these  ion  effects  are 

important  when  r/uj_  is  very  small. 

B 

The  nonlinear  effects  become  important  for  short 
wavelength,  high-k  modes  at  lower  wave  amplitudes  than 
are  needed  for  the  low-k  modes . Hence  if  the  low-k  modes 
are  not  discriminated  against  in  the  parameter  selection, 
one  may  expect  them  to  grow  to  the  largest  amplitudes 
before  nonlinear  growth  saturation  occurs . 

Both  the  theoretical  analysis  and  the  two- 
dimensional  computer  simulations  indicate  that  these 
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ion -electron  cyclotron  interactions  saturate  at  wave 
amplitudes  too  low  for  longitudinal  ion  trapping,  and 
generally  before  there  is  catastrophic  disruption  of  the 
electron  beam.  Therefore  for  Auto-Resonant  ion  accelera- 
tion, it  appears  that  ion-electron  cyclotron  interactions 
are  (i)  probably  unsuitable  as  a wave  growth/ ion  loading 
mechanism,  but  (ii)  unlikely  to  disrupt  the  acceleration 
process  itself. 
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STABILITY  OF  SLAB  ELECTRON 
BEAM  EQUILIBRIA 

In  this  report,  the  stability  of  a certain  class 
of  slab  electron  beam  equilibria  is  discussed.  These 
equilibria  would  describe  relativistic  electron  sheet 
beams  propagated  between  parallel  conducting  plates  along 
a uniform  guide  field.  Such  equilibria  have  been 
utilized  in  some  numerical  simulations  to  investigate  a 
variety  of  electron  beam  phenomena.  A sufficient  condition 
is  determined  for  stability  of  the  slab  equilibria  to  two- 
dimensional  perturbations  which  vary  across  the  slab  and 
along  the  direction  of  propagation.  The  restriction  to 
two-dimensional  perturbations  is  a significant  limitation 
in  the  analysis,  but  this  limitation  is  relevant  for  the 
problem  of  constructing  stable  equilibria  for  two- 
dimensional  computer  simulation,  for  which  the  third 
coordinate  is  ignorable. 
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Equilibrium 


It  is  convenient  to  move  to  a frame  of  reference 
moving  with  the  electron  beam.  The  electron  slab  equi- 
librium is  assumed  to  be  uniform  in  z and  y, 
extending  from  z = -“  to  z = -h=°,  y = -°°  to  y = +°°, 
and  is  bounded  in  x by  conducting  plates  at  x = ±b. 

The  electrons  gyrate  around  an  externally  imposed  magnetic 
guide  field  in  the  z-direction. 

There  are  three  constants  of  the  electron  motion: 

1 ? 

H = — mv  - e4>  (x) 
o 2 o 

e 

P = mv  - — A (x) 
oy  y c oy 

P = mv 
oz  z 


where  the  motion  in  the  beam  frame  is  considered  to  be 
nonrelativistic.  v is  the  electron  velocity,  is  the 


scalar  potential,  A ^ is  the  y-component  of  the  vector 
potential. 


2 


The  distribution  function  f (x,  v)  for  the  slab 

o - 

equilibria  can  be  expressed  in  terms  of  the  constants  of 
the  electron  motion. 

The  stability  analysis  will  be  restricted  to 
distribution  functions  of  the  form 


f (x,  v)  = F (h  + ~ P 2) 
o o \ o 2m  oy  / 


where  a is  a parameter  which  characterizes  the  magnitude 
of  the  mean  electron  velocity  in  the  y-direction.  The 
equilibrium  charge  and  current  densities  are 


N (x)  = d3v  F 
o o 


= -e  / 


d3v  v F = 0 


= -e/c 


d3v  v F 
y o 


a e 

= -N  (x)  — A 

o (1  + a)  me  oy 

= -n  e 'y 

o N y' 


y.  <*/ 


$ and  A 
o oy 


satisfy  the  equations 
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The  equilibrium  electric  and  magnetic  fields, 
are  given  by 


E and  B , 
— o -o 


Stability 

The  equilibrium  is  subjected  to  time  dependent 
perturbations  which  vary  only  in  x and  z: 


E = E (x)  + E.  (x,  z,  t) 
~o  — 1 


4 


B = B (x)  + B_ (x,  z,  t) 

— — o —1 

f = f (x,  v)  + f (x,  Z,  V,  t) 

0—1  — 

In  terms  of  the  perturbed  scalar  potential 
<Mx,  z,  t)  and  vector  potential  A^(x,  z,  t)  : 


where 


1 8-l 

E,  = -7*,  - - — - i 

—1  1 c 9 1 


-1  = V X -1 


9$ 

1 1 

A + - = 0 

—1  c 9 1 


The  perturbed  distribution  function  f^  is 
determined  by  the  linearized  Vlasov  equation: 


df.  / v x B,  \ 9F 
1 e /_  . - -1 \ o 

dt  m\— 1 c / 3v 

9F  dA 

- -a  -2-  P -2  -iZ  + e E, 
me  oy  9w  dt  -1 
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where  d/dt  is  the  time  derivative  following  the 


equilibrium  trajectory  and 


a ■> 

W = H + — P 2 

o 2m  oy 

= -^m(v2+v2+  s2)  - e$ 

2 ' x z / o 

+ I 2 ®L_  A 2 

2 (1  + a)  me2  oy 


s = (1  + 


c)1/!(v  - -2-  2-  A ) 

\ y 1 + a me  oy/ 


Thus , 


Since  d3r  has  now  been  expressed  as  a total 

time  derivative,  it  is  possible  to  obtain  from  Maxwells 
equations  the  following  energy  constant  K: 


- (<>&-  ♦ 

J \8tt  8tt  / 


If  , !!o  S |2 

+ fd3r  d3v  -Li ^ iiLJW q_t  *1±!  2 I 

•/  2 (-3F  /3W) 
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rr  d3v  Trhr  fp-  V 32 


- fd3r  yr-2 — 
_/  (1  + a 


7 * A A,  V • E, 

) 4irmc  oy  ly  —1 


This  constant  K may  in  turn  be  bounded  below  by  the 
quantity 


7 


wmmmm 


and  all  the  perturbed  quantities  must  remain  inf initessimal 
in  order  for  K to  stay  constant  at  its  initially  infini- 
tessimal  value. 

Equations  (7)  and  (8)  are,  therefore,  the 

sufficient  conditions  for  stability  to  two-dimensional 

perturbations  of  slab  equilibria  described  by  Fq (W)  and 

bounded  by  parallel  conducting  plates. 

Examination  of  the  equilibrium  expressions  for 

<V  > reveals  that  the  RHS  of  Eq.  (8)  contains  a term 

proportional  to  N b/y  B , which  must,  therefore,  not 

oL  o oz 

be  too  large  if  this  equilibrium  condition  is  to  be 

guaranteed.  (Here  N is  the  beam  density  as  seen  in 

oL 

the  lab  frame  of  reference.) 
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Limiting  currents  of  an  unneutralized  magnetized  electron 
beam  in  a cylindrical  drift  tube 

J.  R.  Thompson  and  M.  L.  Sloan 

Austin  Research  Associates.  Austin.  Texas  78758 

(Received  30  Au|uil  1977.  find  manuscript  received  2 May  1978) 

Results  of  an  investigation  of  the  steady  state  injection  of  a uniform  unneutralized,  magnetized,  relativistic 
electron  beam  into  a cylindrical  dnft  tube  are  presented.  The  space-charge-litnited  current  and  the 
asymptotic  kinetic  energy  of  electrons  on  axis  it  determined  both  numerically  and  analytically  as  a 
function  of  the  input  kinetic  energy  (y,—  I)nic  1 and  of  the  ratio  of  beam  — to  — wall  radii.  A previously 
cited  "interpolation  formula"  is  obtained  in  the  pencil  beam  limit,  but  more  accurate  limiting  current 
expressions  are  developed  for  other  cases  (such  at  the  fat  beam  limit)  where  the  interpolation  formula  is  as 
much  as  20%  in  error.  The  corresponding  axial  electron  energy  is  also  found  to  be  significantly  smaller 
than  the  previously  cited  value  of  (yi”  — I )mc  1 except  in  the  strong  pencil  beam  limit. 


I.  INTRODUCTION 

Significant  advances  of  recent  years  in  high  voltage, 
pulsed  power  technology  have  stimulated  the  exploitation 
of  high  current  electron  accelerators  for  varied  applica- 
tions. For  many  such  cases,  the  determination  of  space- 
charge -limited  currents  is  of  high  interest  and  has  been 
repeatedly  discussed  in  the  literature  [cf.  Ref.  1 and  the 
bibliography  given  there].  We  consider  herein  a small 
but  relevant  subclass  of  the  general  limiting  current 
problem. 

Specifically,  we  restrict  our  attention  to  the  following: 

(i)  A cold,  cylindncally  symmetric,  unneutralized 
electron  beam  of  radius  a propagating  in  steady -state 
interior  to  an  evacuated  waveguide  of  radius  b along  an 
axial  magnetic  guide  field  which  is  sufficiently  large  to 
render  the  equilibrium  precessional  velocity  and  radial 
excursions  of  the  beam  negligible. 

(ii)  The  beam  is  injected  through  an  anode  which  is 
grounded  to  the  cylindrical  drift  tube.  At  the  anode 
plane,  the  electron  number  density  n0,  the  axial  velocity 
fa,  and  the  particle  energy  (y#-l)nic*  are  assumed  in- 
dependent of  radius  r. 

(iii)  The  ratio  of  the  length  to  radius  of  the  drift  tube 
is  very  large,  so  that  end  effects  are  negligible. 

Here  we  consider  only  those  limitations  on  beam  cur- 
rent which  derive  from  the  conditions  for  propagation  in 
a steady-state  equilibrium  subject  to  the  preceding  as- 
sumptions. We  do  not,  for  example,  consider  limitations 
derived  from  stability  considerations,  although  we  note 
that  previous  work1  has  indicated  that  such  an  electron 
beam  in  a cylindrical  waveguide  is  absolutely  ctable  for 
a sufficiently  large  magnetic  guide  field. 

The  exact  numerical  solution  of  this  limiting  current 
problem  was  obtained  long  ago’  for  nonrelativistic 
(0O«1)  beams,  but  relativistic  solutions  have  only  been 
published  for  the  special  cases  of  pencil  beams  (In b/a 
»1)  or  ultra-relativistic  beams  (y#>  1)  for  which  the 
radial  variations  in  the  downstream  quantities  y,  0,  n 
are  unimportant.  A generalization  [see  Eq.  (17a)]  of  the 
pencil  beam  result  has  been  presented4  as  an  “interpola- 
tion formula’’  for  the  limiting  current  and  has  been  wide- 


ly employed  in  related  calculations  in  the  literature. 
Although  this  formula  indicates  the  correct  zero-order 
parameter  scaling  of  the  limiting  current,  it  is  found  to 
be  quantitatively  in  error  by  over  2U%  in  some  cases. 
Furthermore,  the  expression  for  the  axial  electron  ener 
gy  [see  Eq.  (17b)]  which  is  commonly  associated  with 
this  interpolation  formula  is  grossly  incorrect  (i.e., 
by  more  than  an  order  of  magnitude)  in  some  cases. 

We  describe  relativistically  correct  analytical  ap- 
proximations to  the  limiting  current  problem  for  arbi- 
trary values  of  b/a  and  compare  the  results  to  exact 
numerical  solutions.  The  parameter  range  of  validity 
of  these  approximations  is  delineated. 


II.  GOVERNING  EQUATIONS  AND  NATURE  OF 
NUMERICAL  SOLUTIONS 

Adopting  the  usual  cylindrical  coordinates  r,  8,  z,  we 
consider  the  anode  injection  plane  to  be  located  at  z -0. 
The  injected  beam  current  density  is  then  given  by 

( -noe0„cz  , rsa  , 

J(r,0)  = U) 

( 0 , r>  a , 

where  e is  the  magnitude  of  the  electron  charge. 

Since  the  injection  anode  plane  is  grounded  to  the  con- 
ducting waveguide,  the  potential  function  <p(r,  z). attains 
the  boundary  values 

<t>(r, 0)*(p(6,*)»0  ; (2) 

we  note  than  4>(r,t)s  0. 

The  large  magnetic  guide  field  approximation  implies 
that  the  radius  of  each  electron  is  fixed.  For  steady- 
state  flow,  the  space  charge  in  the  drift  tube  will  create 
fields  to  reduce  the  energy  of  each  beam  electron,  such 
that 

y»y(r,z)*y0+ 0(r,a)  (3) 

with  a corresponding  beam  charge  density 

p(*',*)»-*«»*0ty/(/-'l),/i  , (4) 

where  y0»(l  -0S)'1'1. 
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Far  downstream  of  the  injection  plane  U »6),  the 
variation  in  <b  is  purely  radial  and  is  described  by  the 
radial  Poisson  equation,  which,  from  Eq.  (3),  may  be 
expressed  in  terms  of  y(r): 


l ±fdy.  ^ ' r*a  ’ (5a) 

o,  r >n  , (5b) 

with  boundary  conditions  that  y(0)  be  regular  and  y(b) 

3 y0.  The  parameter  v is  defined 


(6) 


which  differs  by  a factor  of  0 from  the  usual  Budker  pa- 
rameter, so  that  v is  constant  in  r,  z and  proportional 
only  to  the  current. 

By  integrating  Eq.  (5b)  through  the  vacuum  region,  and 
applying  the  boundary  condition  y(6)  -y0,  one  may  con- 
tend only  with  the  single  beam-region  Eq.  (5a)  subject 
to  the  mixed  boundary  condition 


For  given  values  of  y0  and  b/a,  it  is  desired  to  deter- 
mine the  largest  value  of  v (and  hence  /,)  for  which  a 
solution  to  Eq.  (5a)  exists,  subject  to  the  boundary  con- 
dition, Eq.  (7).  This  largest  value  is  called  i/,  and  the 
corresponding  electron  current  is  called  the  limiting 
current. 

The  increase  in  the  right-hand  side  of  Eq.  (5a),  which 
occurs  as  y is  reduced,  reflects  the  axial  density  com- 
pression due  to  electron  slowing.  Equation  (5a)  (with 
boundary  conditions)  requires  that  y(r)  be  monotonically 


n 

FIC.  1.  Dimensionless  current  t>  ss  a function  of  the  rela- 
tivistic factor  yc  of  center  electrons  with  y#1  b/a  fixed. 


TABLE  !.  Numerical  solutions  for  >’»/Ye|. 


7,-2 

To  - 5 

Vo*7 

y„= 100 

iml 

0.553 

3.  12 

4.  98 

97.06 

a 

1.12 

1.26 

1.31 

1.68 

a 

0.221 

1. 14 

1.27 

1.31 

2.05 

1.37 

40.61 

1.74 

*.s 

0.  118 

0.  G95 

1.  12 

22.  90 

a 

1.18 

1.41 

1.48 

1.  89 

-=20 

a 

0.069 

0.  407 

0.  659 

13.  75 

1.21 

1.  50 

1.  61 

2. 14 

increasing,  indicating  that  the  electrons  on  axis  are 
slowed  the  most  below  their  injection  velocity.  The 
previous  relativistic  treatments  of  this  limiting  current 
problem  have  essentially  ignored  the  radial  dependence 
of  this  axial  density  compression,  approximating  the 
right-hand  side  of  Eq.  (5a)  as  a constant. 

In  principle,  if  one  specifies  a value  for  the  y of  the 
electrons  at  the  center  of  the  beam,  ye«y(0),  then  Eq. 
(5a)  determines  y(r)  throughout  the  beam  as  a function 
of  yc,  v,  a.  Applying  the  boundary  condition,  Eq.  (7), 
then  allows  one  to  solve  for  v as  a function  of  ye,  y„, 
b/a.  The  resulting  dependence  of  v on  yc  for  fixed  values 
of  y0,  b/a  is  illustrated  in  Fig.  1.  The  maximum  value 
of  v is  the  desired  t/,,  and  the  corresponding  value  of  ye 
is  yeI.  The  value  of  u at  ye  =1  is  denoted  vfl.  One  may 
show  from  Eqs.  (5a)  and  (7)  that  ye,  > 1 and  v,  > v0. 

Thus,  Fig.  1 is  qualitatively  correct  for  all  cases. 

It  is  of  interest  in  some  applications  to  know  ycl  as 
well  as  vt,  with  a fair  degree  of  accuracy.  For  exam- 
ple, schemes  have  recently  been  proposed5  for  collec- 
tive acceleration  of  ions  which  rely  on  the  controlled 
motion  of  a virtual  cathode,  created  In  an  intense  rela- 
tivistic electron  beam  by  approaching  the  limiting  cur- 
rent conditions.  If  the  motion  of  the  virtual  cathode  is 
to  be  well  controlled,  the  limiting  current  conditions 
must  be  known  precisely. 

It  should  also  be  noted  that  both  the  accessibility  and 
the  stability  of  steady -state  injecticn  is  questionable  for 
the  branch  of  Fig.  1,  where  l < y„  < y„.  However,  the 
branch  yeI  < ye<  y0  may  be  traversed  by  slowly  increas- 
ing the  diode  current  from  zero  while  holding  the  voltage 
constant. 

A set  of  scaled  equations  based  upon  Eqs.  (5a)  and  (7), 
which  allow  very  efficient  numerical  determination  of  the 
limiting  current  parameters,  is  presented  m the  Appen- 
dix. A representative  spectrum  of  the  numerical  solu- 
tions obtained  via  these  equations  is  listed  in  Table  I. 

The  monotonic  dependence  of  u,  and  ycl  on  separate  vari- 
ations of  y0  or  b/a  may  be  observed. 

III.  APPROXIMATE  ANALYTICAL  SOLUTIONS 

In  this  section  we  describe  several  theoretical  ap- 
proaches to  the  limiting  current  problem,  which  yield 
approximate  solutions  that  are  useful  in  different  regions 
of  parameter  space. 
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A.  Iteration  of  integral  equation 

Equations  (5)  and  (7)  may  be  combined  into  a single 
integral  equation,  which  may  be  written  in  iterative 
form  as 


’It)  5 v0  -4v  1 


rfv  V 


d\<  v 


(8) 


where  r * r/a,  0*(l  -y',)utl  and  the  superscripts  in 
parentheses  indicate  the  order  of  the  iteration.  If  y,0,(x) 
is  chosen  to  be  constant,  the  iteration  may  be  readily 
carried  out  to  second  order,  yielding  the  following  equa- 
tions: 

y(0,lr)»y,  , (9) 

y(l,lr)3yc  + (y, -yjx*  , (10) 

\ [y0  -ytl,U))(y.  -yc)  = 2(y,0,-ye0e)lnj 


■y,*elny,y‘" U)ftJ7‘'U)  + i|-l  ’ (11> 


where 


_ yo  + 2y,  1 n(b/a) 
l+ilnWa)  ’ 


(12) 


and 


Equation  (17a)  is  the  previously  mentioned  and  widely 
cited  interpolation  formula.*  The  corresponding  value 
of  ycl  in  Eq.  (17b)  has  also  been  previously  published.1 

The  results  of  Eqs.  (14 )— (1 7 ) are  essentially  first- 
order  iterations,  since  no  use  was  made  of  the  second 
order  Eq.  (11).  Consequently,  they  can  be  expected  to 
have  only  limited  accuracy.  The  two  examples  may  be 
expected  to  bracket  the  correct  solutions,  since  they 
represent  the  two  extreme  choices  for  y„  and  indeed 
this  is  found  to  be  the  case. 


An  example  which  contains  second -order  accuracy  is 
to  let  y,  be  such  that  ya,(l)*yf.  From  Eqs.  (11)  and 
(12),  one  obtains 


(y«0,  -rs8e)  \1  +2  In (b/a)j 


and  hence 


(y  _ l)J/* 

1/0  [1+2  ln(f>/a)|  [y#  + 1 +4  ln(6/a)|1/z  ' 

Although  Eq.  (18)  is  sufficiently  complex  so  that  it  does 
not  yield  closed  form  algebraic  solutions  for  v,  and  ycl, 
it  is  still  fairly  useful  since  it  may  be  numencall  maxi- 
mized in  minutes  on  the  presently  available  scientific 
hand  calculators.  Also,  Eq.  (19)  provides  a lower 
bound  on  v,. 


• (13) 

Here  and  in  the  following  text,  the  y’s  and  0’s  with  like 
subscripts  or  superscripts  are  related  in  an  obvious 
manner. 

There  remains  the  freedom  to  specify  the  constant  yt. 
For  example,  if  one  selects  y,3y„,  then  Eqs.  (12)  and 
(13)  yield 


flo^o  -y,) 

1+2  \n(b/a) 


and  hence 


(14) 


Equation  (18)  may  further  be  used  to  determine  the 
bounds  of  validity  of  Eq.  (16).  If  one  assumes  the  usual 
pencil  beam  limit  y,-y(<--  yc  (i.e.,  the  fractional  varia- 
tion in  y across  the  beam  is  small),  Eq.  (18)  reduces 
to  Eq.  (16).  This  pencil  beam  ansatz  specifically  re- 
quires the  ordering  yJ/J  « 1 + 2 In b/a  for  validity.  There- 
fore, Eqs.  (1 6)— (17)  are  essentially  accurate  only  in  the 
pencil  beam  limit,  which  is  the  lower  left-hand  portion 
of  Table  I. 

Another  example  containing  second-order  accuracy 
is  to  let  y„  be  such  that  y(,,(0)  =ye.  Then,  Eq.  (11)  gives 
the  result  (also  obtained  independently  elsewhere)': 


v,  «!/„■ 


6o(ro  - 1 ) 

1 + 2 ln(b/a)  - 


r«iml 


(15a) 

(15b) 


‘■■(r.-r.t  [<r.  S.  -r.t.) -(nThfew) 


(20) 


which  is  a previously  published  result.* 

Another  possibility  is  to  select  ykmyt  to  obtain 


l+2ln(6/a) 

(16) 

and  hence 

i +2  ln(6/e)  • 

(17a) 

v mvl/* 

Ycl  Y 0 » 

(17b) 

o 

a 

© 

(17c) 

with  y,  specified  by  Eq.  (12).  Therefore, 

v0*[l  ♦2ln(6/a)l',(y0-l),^[yJ-l  +4(y#  - 1)  ln(6/a)ll/* 

+ ln(l  + 1^-1M(yt-1)ln(6/a)|1/,+y»-1)r  (21) 

V 1-2  tn(6/a)  /J 

and  as  before  i/t  provides  a lower  bound  on  v„  but  Eq. 
(20)  must  be  numerically  maximized  to  accurately  de- 
termine v,  and  yel.  Equation  (20)  may  also  be  seen  to 
have  the  pencil  beam  limit  of  Eq.  (16). 

If  Eq.  (20)  is  expanded  in  the  limit  y,  »yc»l,  so  that 
the  variation  in  y across  the  beam  is  large,  one  obtains 
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I 

I 
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and 


2yJ,  » ln(y-/ycl)  + 1 +2  In (b/a)*  • 


(23) 


C.  Ultrarelativistic  perturbation  expansion 

If  we  restrict  our  attention  to  the  case  of  a sufficiently 
highly  relativistic  beam  where  yB  » 1 and  2yj,  »t,  so 
that  we  may  use  the  expansion 


The  expansion  ansatz  specifically  requires  that 
1 ♦ 2 ln(6/n) « yj/s  (hereafter,  denoted  the  fat  beam  limit) 
which  is  the  upper  right-hand  portion  of  Table  I. 

In  summary,  Eqs.  (18)  and  (20)  are  accurate  second- 
order  representations  of  fairly  wide  validity;  both  equa- 
tions have  the  correct  pencil  beam  limit,  and  Eq.  (20) 
also  has  the  correct  fat  beam  limit. 

B.  Variational  principle 

A variation  principle  which  generates  the  differential 
Eq.  (5a)  and  the  boundary  condition  (7)  is 


1 /0  = r(y*  - 1 )-,/*  * 1 ♦1/2/*...  , 

Eq.  (5a)  may  be  solved  iteratively  to  yield 
ur1 


y(r)  -yc* 


-LmA+JHV 
2y*  V J 


(28) 


(29) 


Then,  the  boundary  condition,  Eq.  (7),  gives 
Yo  = Yc  *v[l  + 2 ln(h/o)| 


+ ^-[Wl  + — \ + 2l/  ln(- 

2V«L  V Yc)  v+Y'  \a, 


(30) 


fy(x)]  = J +4i/(y* -l)1/2j 


The  peak  value  of  v (i.e. , v,)  may  be  found  by  differen- 
tiating Eq.  (30)  with  respect  to  ye  and  setting  Bv/&yc  =0, 
to  obtain 


y(x)[ 2y„  - ypt)] 
2 In (b/a) 


with  respect  to  possible  trial  functions  y(x).  One  ex- 
pects that  fairly  high  accuracy  may  be  achieved  on  the 
limiting  current  with  only  a modestly  accurate  guess  at 
the  trial  function  y(x). 

A trial  function  of  modest  accuracy  is  that  of  Eq.  (10), 
which  contains  the  two  free  parameters  ye  and  y,.  These 
two  parameters  are  then  chosen  to  render  the  variational 
principle,  Eq.  (24),  stationary.  Carrying  out  this  pro- 
cedure and  using  the  fact  that  v,(ye,y0,  b/a)  is  itself  an 
extremum  at  yc=yel,  one  obtains  the  following  paramet- 
ric solution  for  the  limiting  current  parameters. 


»i“(Yh  -Yci?[y.,(Ym0. i ~Yel0cl) 
Qm  2‘'jly„0ll-Yg,0cl) 


J 

1 1 

2^|=ln(l  + M 

1 + 2 ln/-Y| 

1 , (24) 

V Yc) 

‘'i+Yc,  [ 

WJ 

1 Xal 

, Zt/'  y“ , 

I,  which  is  to  be  extremized 

(i'i+YciF 

* 

(31) 


Equations  (30)  and  (31)  may  now  be  solved  simultaneously 
to  determine  the  limiting  current  parameters  v,  and  ye). 

1.  Pencil  beam  case 

For  this  limit,  we  assume  the  ordering 
2lnU>/a)»y§/J»l 

and  make  the  ansatz  u,  «ycl.  We  then  obtain 


Yel*Y,o/,-iYi'/3-4 


-ZlL 


In  (b/a) 


1 2ln(6/a)  21n(6/a) 


(32) 

(33) 

(34) 


(25) 

(26a) 


2t/, 


~ 0.i(y„ -Yel)  ' 
To  * r.i  *Q(y„  -yc,)ln(6/a) 


(26b) 

(27) 


These  equations  may  be  solved  by  first  specifying  trial 
values  for  y„  and  y„,  subject  only  to  y„  > yel  > 1,  and 
then  solving  Eqs.  (25 )— (27)  sequentially  for  the  corre- 
sponding values  of  v„  b/a,  yB.  Trial  values  of  yel  and 
y„  which  result  in  b<  a must  be  rejected.  One  may 
again  be  forced  to  make  several  calculations  on  a hand 
calculator  to  generate  the  desired  values  of  y„  and  b/a. 
The  resulting  solutions  are  found  to  be  quite  accurate 
except  in  the  fat  beam  regime. 

2035 


which  agrees  with  the  pencil  beam  results  of  Eq.  (17), 
when  the  latter  are  expanded  under  the  pencil  beam  or- 
dering of  Eq.  (32).  This  ordering  is  also  sufficient  to 
validate  the  ansatz  i>,  «yc, . As  we  noted  previously,  the 
pencil  beam  regime  is  the  lower  left-hand  portion  of 
Table  I;  however,  Eqs.  (33)  and  (34)  also  require  yj/s 
» 1 for  validity,  implying  a shift  toward  the  lower  cen- 
ter of  Table  I. 

2.  Fat  beam  case 

For  this  limit  we  assume  the  ordering 
3+2  ln(6/a)  « yJ/J  (35) 

and  make  the  ansatz  vx  »ycl.  Then,  Eq.  (31)  becomes 
2y*,  «ln(i/,/ycl)  + l +2  ln(6/o)  + . . . (36) 

while  Eq.  (30)  implies 

1 • . <"> 

To  obtain  ycl,  Eq.  (37)  may  be  inserted  in  Eq.  (36), 
which  may  then  be  solved  iteratively.  Carrying  the  iter- 
ation to  second  order  for  improved  accuracy,  one  ob- 
tains 
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y ..  ' 


) l/z 

+ 1 >21n(6/a) -ln[l  ♦2ln(6/a)]|  . (38) 


tivislic  (at  beam  limit  illustrated  in  the  last  column. 
There,  the  expansion  formulae  of  Eqs.  (37)  and  (38)  are 
best,  with  the  second-order  iteration  Eq.  (20)  also  being 
fairly  accurate.  The  ultra-relativistic  pencil  beam  ex- 
pansions are  accurate  only  for  very  large  values  of  b/a, 
as  assumed  in  Eq.  (32),  and  so  they  have  been  omitted 
from  Table  n. 


Once  ytl  has  been  computed  from  Eq.  (38).  then  v,  may 
be  determined  from  Eq.  (37).  One  finds  that  the  order- 
ing of  Eq.  (35)  is  sufficient  to  validate  the  insatz  u, 

It  may  also  be  seen  that  Eqs.  (22)  and  (23)  agree  with 
Eqs.  (36)  and  (37)  for  the  fat  beam  limit. 

When  y,  »1,  the  value  of  yel  given  by  Eq.  (38)  for  the 
fat  beam  limit  may  be  seen  to  be  considerably  less  than 
the  value  given  by  Eqs.  (17b)  or  (33)  for  the  pencil  beam 
limit.  Consequently,  it  would  be  grossly  inaccurate  to 
assume  that  these  pencil  beam  equations  for  yel  have  a 
general  validity  as  has  been  ascribed  to  the  interpolation 
formula  for  v„  Eq.  (17a). 

D.  Accuracy  of  the  analytical  approximations 

Exact  numerical  solutions  which  have  been  obtained 
and  which  are  illustrated  in  Table  I may  be  used  to  gauge 
the  accuracy  of  the  various  analytical  approximations. 
Comparisons  for  several  test  cases  are  illustrated  in 
Table  II. 

It  may  be  seen  that  the  second-order  iteration  solu- 
tions of  Eqs.  (18)  and  (20)  bracket  the  numerical  results 
much  more  tightly  than  do  the  first  order  iteration  solu- 
tions of  Eqs.  (15)  and  (17).  In  particular,  the  "interpo- 
lation formula"  Eq.  (17a)  is  observed  to  significant  un- 
derestimate the  limiting  current,  while  the  correspond- 
ing Eq.  (17b)  grossly  overestimates  yel. 

The  results  of  the  variational  calculation  are  seen  to 
be  the  most  accurate  overall,  except  in  the  ultra-rela- 
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APPENDIX:  EQUATIONS  FOR  A NUMERICAL 
INVESTIGATION  OF  LIMITING  CURRENT 

A simple  computational  procedure  may  be  used  to  in- 
vestigate Eq.  (5a)  and  directly  determine  the  value  of 
v,t  Yd • More  precisely,  in  this  computational  procedure 
we  choose  ycl,  b/a  and  determine  i/(,  y0  therefrom.  Us- 
ing the  change  of  variables  y =(i>rl/az),  Eq.  (5a)  may  be 
written 


_i  „ = r 

9y  y 9y  (y2  - l),/z 


(Al) 


By  requiring  that  y be  regular  at  the  origin  and  attain 
the  value  y(0)  = yc,  Eq.  (Al)  may  be  integrated  numerical- 
ly to  generate  a single  parameter  family  of  solutions 
y(y,  ye).  Using  the  boundary  condition  Eq.  (7)  and  evalu- 
ating y at  r -a  (i.e. , y - u),  we  obtain 


Thus,  given  b/a  and  yc,  the  direct  integration  of  Eq. 
(Al)  with  the  subsidiary  condition,  Eqs.  (A2),  yields 
the  general  solution  y,(6/<j,  yc,  v). 


TABLE  II.  Test  of  accuracy  of  analytical  approximations. 


V,  * 5,  b/a  • 2 

y,*7,  b/a-  4.985 

Vo*  100,  b/a-1 

V|  Tel  T.I 

*<1  7el  T„ 

v,  T c 

Numerical  computer 


results 

1.273 

1.313 

e e • 

1.125 

1.460 

. e • 

97.06 

1.676 

Interpolation  formula,  Eqs. 
(17)  (lst-order  iteration) 

1.12 

1.71 

3.09 

1.03 

1.91 

3.12 

93.12 

4.64 

Equations  (15)  (lst-order 
Iteration) 

1.84 

1.00 

2.68 

1.41 

1.00 

2.42 

99.00 

1.00 

Equation  (18)  (2nd -order 
Iteration) 

1.31 

1.27 

2.83 

1.13 

1.47 

2.78 

96.28 

1. 15 

Equation  (20)  (2nd-order 
Iteration) 

1.26 

1.38 

2.90 

1.12 

1.53 

2.83 

96.96 

1.82 

Equations  (25) -(27)  (varia- 
tions! principle) 

1.28 

1.31 

2.90 

1.  13 

1.49 

2.83 

97.82 

1.33 

Equations  (37) -(38)  (ultra  - 
relativistic  fat  beam) 

1.35 

1.11 

... 

1.09 

1.39 

... 

97.  13 

1.59 
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In  order  to  determine  the  limiting  current  p,,  we 
must  determine  the  value  of  ye  such  that  Bu/9yc  = 0. 
However,  since  y0  - y^b/n,  >c,  p)  we  may  write 


0=  *Zi*L 

dye  8p  8yc 


(A3) 


Therefore,  the  condition  of  limiting  current,  8i//ayc  = 0, 
corresponds  to  a determination  of  the  parameters  b/a, 
yc,  v such  that  8y0/8 yc  = 0.  From  Eq.  (A2)  a straightfor- 
ward differentiation  yields 


Equations  (Al)  and  (A5)  may  now  be  simultaneously 
integrated  outward  in  y until  a value  y is  reached 
such  that  (8y0/8yc),  as  given  by  Eq.  (A4).  vanishes.  The 
limiting  current  v,  is  then  given  by  i/,  *y„  with  y0  de- 
termined from  Eq.  (A2).  The  initially  specified  quantity 
yc,  of  course,  corresponds  to  ycl  for  this  value  y=y, 

»i/,.  It  may  be  demonstrated  that  such  a point  y,  >0 
always  exists  for  each  choice  of  yel>  1. 


where 

T»9y{y,ye)/dye  = T{y,ye)  . 

The  equation  governing  T may  be  obtained  similarly  from 
Eq.  (Al)  by  differentiation  with  respect  to  yc: 

V-iVn  ’ (A5) 

where  T(y  =0)  = 1 and  T is  regular  at  the  origin. 
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RELATIVISTIC  ELECTRON  BEAM  CYCLOTRON  MODE 
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APPENDIX  D 

I 

INFLUENCE  OF  SPATIAL  VARIATION  IN  y | 

ON  THE  RELATIVISTIC  ELECTRON 
BEAM  CYCLOTRON  MODE 

I.  Introduction 

In  a cold  relativistic  electron  beam  moving  with 

velocity  u -*■  c in  a vacuum  along  a uniform  guide  field 

B in  the  z-direction,  electron  beam  cyclotron  modes  can 
oz 

propagate  with  frequency 


where 


Q 


eB 

oz 

m yc 
e 


/ 


4tt  N e 
o 


m Y 

e 


2 


J 

I 

I 


1 


-1/2 


k2  = k 2 + k 2 


I 

1 

I 

r 


and  the  beam  density  N together  with  y are  assumed 

o 

constant. 

In  cylindrical  geometry,  there  exists  (for  a given 

value  of  k and  mode  number  m)  an  infinite  denumerable 
z 

set  of  modes  having  different  values  of  k . For  the 

i 

cylindrically  symmetric  m = 0 mode,  the  lowest  value  of 
k^  is  k^a  » 2.95,  where  a is  the  beam  radius  and  the 
beam  is  bounded  by  conducting  walls. 

The  perpendicular  electron  velocity  perturbation 
v^  is  proportional  to  ^1/y(oj  - k^u  + ft  )j  and  remains 
finite  across  the  beam.  If  y varies  in  magnitude  from 
the  inside  to  the  outside  of  the  beam: 


spatial  variations  are  introduced  in  ^1/y(oj  - k u + ftjj 
which  may  not  be  negligibly  small.  In  particular,  for  the 
case  when  the  parameter  a - 1,  the  spatial  variations  are 
significant  and  must  be  included  in  the  analysis  of  the 
electron  cyclotron  beam  dispersion  relation. 

In  this  report,  the  propagation  of  electron  beam 
cyclotron  modes  in  the  presence  of  spatial  variations  in 
Y is  investigated.  The  eigenvalues  are  calculated  using 
a variational  approach.  This  obviates  the  necessity  to 
obtain  exact  solutions  of  the  eigenvalue  equations.  It  is 
established  that  in  comparison  with  constant  Y»  the 
corresponding  eigenvalues  (k^a)  increase  as  the 
parameter  a increases.  The  larger  the  eigenvalue  (k^a) 
the  closer  is  the  cyclotron  motion  of  the  electron  to 
resonance  with  the  oscillating  fields  of  the  corresponding 
eigenmode.  In  addition,  the  perpendicular  perturbed 
currents  and  electric  fields  become  more  sharply  peaked  on 
the  surface  of  the  beam. 
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II.  Electron  Beam  Cyclotron  Mode 


II 


The  electron  beam  dynamics  is  described  by  the 
cold  fluid  equations: 


|S  + 7 

3 1 


nv  = 0 


eE (v  x B) 

— c — — 


and  closure  is  obtained  through  Maxwell's  equations: 


V x £ = - - T= 


1 3B 
c 3t 


7 x B 


1 3E  4tt 


7 • E * 4rrp 


7 • B * 0 


n 

i 


**  •• 

Ski.. 


’I 


where 


s 


p = - en 


J = - env 


Y = 


The  velocities  and  electromagnetic  fields  of  the 


beam  equilibrium  are: 


v = ( 0,  v 

V 


2-rrN  er 
o 

oe  ” by2 

02 


°'  uo) 


E = ( E = -2irN  er , 0,  0) 
-to  \ or  o / 
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(°, 


B 


09 


2ttN  eur  1 

, B i 

c oz> 


Y * Y,Jl  <r*  - a2)(  >>  1 


J 
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The  magnitude  of  the  parameter  a is  a measure  of  the 

spatial  variation  in  y. 

Let  the  perturbations  be  of  the  form 

exp[-iwt  + ik  z + im0 ] . The  following  differential  equa- 
z 

tions  describe  the  electron  beam  cyclotron  modes: 


13  3A  (m  + 1) 2 _+ 

n — A 


(1) 


1 3_  . 3A_  (m  - 1) 2 - 

S H 5 H 52  A 


p2a2A  = 0 


(2) 


f 


2 


The  boundary  conditions  on  the  conducting  wall  at 
5=1  are 


A+  = a"  (3) 


i r 

! j' 

' 
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i 
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I 

I 

I 

I 
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V 


il 


Equations  (1)  and  (2)  are  solved  for  A+,  A , 

and  after  substitution  in  Eqs.  (3)  and  (4) , the  eigenvalue 

equation  for  Q is  obtained.  In  the  low  phase  velocity 

limit  k 2 > (jo2/cz,  the  frequency  oj  is  given  by 
z 


where 


k2 


4 

a (Q  - l)a2 


(6) 


The  smaller  the  magnitude  of  tu  2/k2c2,  the  closer  is 

po 

the  cyclotron  motion  of  the  electrons  to  strong  coupling 
or  resonance  with  the  mode.  In  fact,  it  may  be  seen  from 
Eqs.  (1)  and  (6)  that  the  radial  location  of  the 
resonance  is 


C2 


1 + 


4 

ak2a2 


which  is  just  outside  the  beam  edge. 


(7) 
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The  perturbed  electron  velocity  and  electro- 
magnetic fields  are: 


v = v ± iv 
r 9 


my  fui  - k u ± n 
e \ z 
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I 
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C(A+  + A~) 
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, r E ik  c 

> - - TT — r P f + -S-  (A 

r (kzu  - u)  La  5 z 


4*  — 

- A ) 


1 fie  3EZ  V , + - 1 

[t  ?r + t ( + A 


III.  Eigenvalues 


The  solutions  of  Eqs.  (1)  and  (2)  for  a = 0 have 
been  derived  in  previous  reports.  They  are  restated  here 
for  completeness  and  easy  reference. 

When  a = 0 , 


A « a J ( <a£ ) 
o m+l 


A = b I (PC) 
o m-l 


where 


-k  2 + lim  ■ 2 

z aQa 
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J , . and  I are  Bessel  and  modified  Bessel  functions, 

m+1  m-1 

respectively. 

The  eigenvalue  equation  for  ica,  obtained  after 
substituting  these  solutions  in  Eqs.  (3)  and  (4) , is 


WKa)  We*’ 

2 + <a  — ; r + pa  7— ■-  = 0 


J (<a) 
m+1 


I i (pa) 
m-l 


(8) 


where 


J '(tea) 
m+1 


t — J , . (tea) 
d<a  m+1 


Equation  (5)  determines  the  frequency  with  values  of  k 
given  by 


k2  =*  (k  2 + k2) 


If  pa  < 1,  m = 0,  Eq.  (8)  reduces 


J (tea) 

tea  » - 2 

^ (<a) 
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I 


so  that 


I 

I 

I 

I 

r 


<a  = 2.95,  5.84 


If  pa  < 1,  m = -1,  Eq.  (8)  reduces  to 


<a  = 


(<a) 

J (<a) 
o 


so  that 


<a  = 1.91,  4.60, 


There  exists  an  infinite  denumerable  set  of  eigen- 
values for  each  value  of  m. 

When  a i 0,  the  solutions  of  Eq.  (1)  are  not 
expressible  in  terms  of  simple  functions.  However,  a 
variational  approach  may  be  adopted  to  estimate  the  eigen- 
values, and  exact  solutions  of  Eq.  (1)  need  not  be  known. 

+* 

Let  Eq.  (1)  be  multiplied  by  £A  (£)  and 
integrated  from  £ = 0 to  £ = 1.  After  imposing  the 
boundary  conditions  at  £ = 1,  namely,  Eqs.  (3)  and  (4) , 
it  may  be  verified  that 


- lb  I2  I (pa)  |"2I  .(pa)  + pa  I (pajl  * A (A+  ,A+) 

1 o 1 m-i  l m—  1 m-l  j 

(9) 


where 
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I! 
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I 

i 


»(»**,  a+)  , f 5||f|  as  *j  |A*|2  de 


+ q2a2 


/ 5iA+ia«  - ; / 


+ l 2 


5 A . 

d5  ,10) 


and  the  solution  A = b I , (pa£)  of  Eq.  (2)  has  been 

o m- 1 

utilized. 

The  first  variation  of  the  functional  A with 
respect  to  A+  is  stationary,  that  is,  6 A = 0,  if  A+ 
is  a solution  of  Eq.  (1) 

Thus  by  selecting  a trial  function  A+  which 
satisfies  the  boundary  conditions  and  for  which  6A  = 0, 
an  approximate  dispersion  relation  can  be  obtained  from 
Eq.  (9). 

Consider  the  following  trial  function  which  is 
suitable  for  estimating  the  eigenvalues  of  the  lowest 
eigenmode : 


A+  - 5 lra+1l  I' g + g252  + g45 


*) 


From  Eqs.  (3)  and  (4)  , 
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32  " bo[-  2g  + 2 K6  + lm  + 1l)Vl,'al 


+ Pal  1 (pa) ! 
in-1  ) J 


94  ‘ bo[9  ' |!(4  + lm  + 1l)Im-l(pa)+  paIm-l(pal|] 


where 


g 5 g /b 

o o 


A is,  therefore,  a trial  function  with  one 


+ *f  * 

parameter  g.  Substituting  for  A in  A (A  , A ) , 


A = lb  2 A (m,  a,  qa,  pa,  Q,  g) 
I o o 


The  parameter  g = g is  chosen  such  that 


3AQ(g) 


(ID 


3g 
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The  eigenvalue  equation  for  Q is  then  given  by 

A 

(m,  a,  qa,  pa,  Q,  g) 

+ I (pa)  [21  (pa)  + pa  I (pa)l  = 0 . (12) 

m— 1 L m-1  m— 1 J 

Equation  (12)  is  solved  numerically  to  determine 
Q and  hence  the  eigenvalue  ka.  Some  representative 
solutions  for  the  cylindrically  symmetric  m = 0 eigenmode 
with  eigenfunction  A+  having  no  nodes  are  tabulated  in 
Figures  1,  2,  and  3. 

Figures  1 and  2 display  the  increase  in  magnitude 

of  the  eigenvalue  ka  for  given  values  of  w 2a2/c2 

po 

and  k^a  as  the  parameter  a increases. 

Figure  3 displays  the  increase  in  magnitude  of  the 

eigenvalue  ka  for  given  values  of  id  2a2/c2  and  a as 

po 

k a increases.  It  may  be  seen  that  k a = \/k2a2  - k 2a2 
z l z 

increases  strongly  with  increasing  k a for  a = 1. 

z 

The  accuracy  of  the  variational  calculation  may  be 
judged  by  comparison  with  analytic  solutions  of  Eq.  (1) 


15 


and  (2).  Thus  for  m = 0,  co  2a2/c2  <<  1, 

po 

w2a2/c2  <<:  k^2a2  <<  1,  and  a = 0.5,  the  solution  of 
Eq.  (1)  is 

A+  = aQC(Q  - ?2) 

and  the  corresponding  eigenvalue  is  ka  = 4.  The 

variational  solution  for  m = 0,  u>  2a2/c2  = 0.06, 

po 

kz2a2  = 0.01,  and  a = 0.5  is  ka  = 4.01.  The  eigen- 
values obtained  from  the  variational  approach  are  indeed 
acceptable. 

+ +*  + 
The  trial  function  A which  minimize  A (A  , A ) 

are  plotted  as  functions  of  5 in  Figure  4.  A+(5)  is 

normalized  so  that  A (5=1)  =1.  The  associated  fields 

E (5)  and  E (5),  computed  from  the  trial  function 

Z IT 

A (C) , are  plotted  in  Figures  5 and  6.  The  profiles  of 
A (5)  and  E^C?)  are  not  appreciably  altered  with 
changes  in  the  value  of  a.  However,  E^(5)  becomes 
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- r- 

r * 


m = 0 


more  sharply  peaked  at  £ = 1 when  a > 0.5,  its 
magnitude  at  £ = 1 being  proportional  to  k2a2  which 
increases  with  a. 

IV.  Discussion 

The  spatial  variation  in  y was  expressed  in 
terms  of  the  parameter  a: 


Y 


au> 

_ 29 
4c2 


(r2 


) 


where  a was  assumed  to  be  in  the  range 


a‘u) 

po 


A relativistic  electron  beam  of  constant  Yq  entering  a 
cylindrical  waveguide  bounded  by  conducting  walls  will 
produce  a beam  equilibrium  having  a spatial  y variation 
with  a - 1. 


22 


The  magnitude  of  corresponding  eigenvalues  ka 
increases  with  increasing  a.  For  m = 0,  the  value  of 
ka ^(u^o2a2/c2  < 1,  k^a  < lj  for  the  lowest  eigenmode 
(no  nodes)  increases  from  2.95  to  5.77  as  a 
increases  from  zero  to  unity.  The  larger  the  magnitude 
of  ka  the  closer  is  the  cyclotron  motion  of  the  electron 
to  resonance  with  the  oscillating  fields  of  the 
corresponding  eigenmode.  Larger  ka  values  impose  some- 
what more  stringent  restrictions  on  the  temperature  of  the 
electron  beam  in  order  to  avoid  significant  "Landau" 
damping,  since  Eq.  (7)  shows  that  the  spatial  location  of 
the  cyclotron  resonance  then  occurs  much  closer  to  the 
beam  surface.  Hence,  a finite  energy  spread  would  more 
easily  allow  a few  particles  to  be  in  resonance  with  the 
wave. 


Furthermore,  E^{£)  and  consequently,  the 
perturbed  perpendicular  electron  velocities  exhibit  a 
transition  in  profile  to  one  which  is  sharply  peaked  at 


J 


C = 1 when  a > j (m  = 0) . The  surface  perturbations  are 
then  larger  than  the  interior  perturbations.  The  effect 
of  nonlinearities  are,  therefore,  likely  to  be  more 
significant  on  the  surface  than  in  the  interior  of  the 
beam. 
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RELATIVISTIC  ELECTRON  BEAM  PLASMA  HEATING  (20-30  min.) 

A Synopsis 

• I.  Introduction  (5-7  min.) 

wherein  the  problem  of  heating  101H-i0ir  density  plasmas 
with  beams  is  posed--and  it  is  assumed  that  one  gives  a damn 
about  efficient  (or  rapid)  heating — as  for  single  pass,  0-pinch 
geometry — and  so  optimization  of  the  heating  interactions  is 
desired 

II.  Return  Current  Heating  (5-7  min.) 

wherein  it  is  argued  that  return  current  heating  goes  to  pot 
for  (iip/ng)  >100 

III.  Primary,  hydro  dynamical,  beam-plasma  interaction  (10-15  min.) 
wherein  the  scaling  obtained  from  computer  runs  is  discussed, 
and  it  is  argued  that  reasonable  deposition  efficiencies  (>107.) 
can  be  obtained  even  for  (np/nB)~105  provided  that  the  beam 
quality  assures  a hydrodynamical  interaction:  9rm3£  2(nfi/iip) 
This  quality  may  be  easier  to  achieve  than  the  MA/cm2  current 
density  required  for  return  current  heating. 
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RELATIVISTIC  E-BEAM  PLASMA  HEATING 


W.  E.  Drummond  and  J.  R.  Thompson 

✓ 

(Notes  for  a 1/2  hour  presentation  at  the  Fourth  National 
School  on  Plasma  Physics,  Novosibirsk,  July  1974) 

It  is  well  known  that  advances  in  high  energy,  pulsed  power 
technology  have  made  available  extremely  powerful  electron  beams, 
capable  of  delivering  .1-1  megajoules  of  directed  beam  energy  within 
.1  psec.  My  remarks  today  will  be  confined  to  the  prospects  for 
utilizing  such  beams  for  the  heating  of  target  plasmas  of  moderately 
high  density,  101*  -10 1 ® cm-*,  in  volumes  of  10s  -10s  cm*,  to  multi- 
kilovolt temperatures,  [i  will  not  discuss  pellet  heating,  and  these 
remarks  will  be  less  relevant  for  the  heating  of  lower  (<10l ’’cm” *) 
density  plasmas.] 

I will  further  presume  that  it  is  either  necessary  or  desirable 
that  the  heating  be  accomplished,  rapidly  and  with  a resonable  efficiency 
of  energy  deposition  from  the  beam  into  the  target  plasma.  For 
example,  this  would  be  the  case  in  any  single  pass  heating,  in  which 
the  lifetime  of  interaction  for  a beam  electron  would  be  limited  to 
its  transit  time  of  a few  nanoseconds.  In  this  case  it  is  important 
to  study  the  alternative  mechanisms  of  beam-plasma  interaction,  so 
that  parameters  might  be  selected  to  optimize  one  or  another  mechanism. 

(There  is  another  class  of  problems,  which  will  not  be  further 
discussed  here,  in  which  the  primary  concern  is  how  to  inject  and 
confine  a powerful  beam  in  a toroidal  device  for  many  thousands  or 


vwW  m 


perhaps  millions  of  plasma  periods.  For  this  problem,  the  extended 
time  of  interaction  would  permit  almost  any  of  the  interaction 
mechanisms  to  couple  much  of  the  beam  energy  into  the  target,  and 
the  question  of  optimizing  one  or  another  interaction  would  be  much 
less  important.) 

There  exists  an  increasing  body  of  experimental  data,  indicating 
the  successful  one  pass  heating  of  target  plasmas  with  energy  deposi- 
tion efficiencies  in  the  range  of  10%-50%,  for  the  case  where  the 

target  plasmas  were  no  more  than  a few  hundred  times  denser  than  the 
* % * 
electron  beams.  The  experiments,  in  agreement  with  theory,  have 

shown  increasing  energy  deposition  efficiencies  as  the  beam  to  target 

density  ratio  is  increased  toward  one.  Therefore,  in  response  to  the 

demand  for  heating  plasmas  of  higher  and  higher  density,  there  has 

been  a certain  trend  toward  increasing  the  beam  current  (if  necessary 

at  the  expense  of  particle  energy)  and  reducing  the  beam  size  to 

further  increase  its  density.  This  evolution  becomes  increasingly 

more  difficult  and  technologically  more  demanding  as  one  considers 

heating  the  higher  density  target  plasmas  which  may  be  desired  in 

future  controlled  fusion  experiments.  For  example,  to  heat  a 10 1 ‘cm' 

density  target  plasma  with  relativistic  beam  of  density  above  27»  of 

the  target  density,  one  must  be  able  to  propagate  current  densities 

above  a megampere/cm2  for  a distance  of  many  cm.  Such  a beam  would 

be  much  more  intense  than  those  commonly  generated  at  particle 

densities  of  10 12  -101*  cmr*,  or  current— densities  of -some  kA/cm2. 

It  is  an  appreciation  of  the  difficulty  of  this  chore  of  pro- 
pagating beams  of  such  high  current  density,  as  one  scales  up  to  dens 


fusion  experiments,  which  has  motivated  us  to  examine  the  density 

scaling  of  the  energy  deposition  efficiencies  for  various  mechanisms 

of  beam-plasma  interaction,  so  as  to  determine  the  minimum  required 

- beam  current  density  and  any  alternatives  to  very  high  beam  density 

which  may  still  allow  adequate  heating. 

Let  us  suppose  that  the  target  plasma  has  been  preionized,  or 

else  that  the  neutral  gas  pressure  is  sufficiently  high  (above  ~30  mT 

or  1015cm"*)  that  field  induced  cascade  ionization  produces  a fully 

ionized  target  early  in  the  beam  pulse  time.  The  induction  fields 

of  the  electron  beam  will  quickly  drive  a return  current  carried  by 

•2  2 2 

the  plasma  electrons,  provided  that  P,  - ^ nP  a - >>  1.  (Here  is 

c 3X10 11 

the  target  electron  plasma  frequency,  a is  the  beam  radius,  and  the 
condition  has  been  shown  by  Sudan  and  others  to  be  roughly  appropriate 
with  or  without  external  magnetic  fields . ) 

One  mechanism  for  depositing  beam  energy  into  the  target  plasma 
is  via  the  dissipation  of  this  return  current  due  to  micro  turbulent 
interactions  between  the  return  current  electrons  and  the  target 
plasma  ions.  These  interactions  extract  drift  energy  from  the  plasma 
electrons,  which  is  replenished  by  the  induction  fields  which  transfer 
energy  from  the  beam  electrons  to  the  plasma  electrons.  The  scaling 
of  these  microturbulent  interactions  is  illustrated  on  the  first 


slide. 


The  first  point  to  notice  is  that  the  return  current  drift  velocit 
scales  with  the  hearn  to  target  plasma  density  ratio,  and  that  heating 
rate  per  target  electron  scales  as  the  inverse  square  of  the  plasma 
density.  This  is  because  the  return  current  drift  energy  is  what 
directly  feeds  the  turbulence.  Consequently,  the  rate  of  return  current 
dissipation  falls  rather  dramatically  as  one  attempts  to  heat  target 
plasmas  much  denser  than  the  beam. 

The  second  point  to  notice  is  that  there  is  evidence  that  the 
"effective  collision  frequency,"  whose  average  value  influences  the 
L/R  current  decay  time,  is  not  constant  but  in  fact  decreases  as  the 
plasma  heats  up.  For  example,  we  and  others  (Morse,  Biscamp)  have 
performed  one-dimensional,  particle  computer  calculations  of  electron- 
ion  interactions,  in  which  the  initial  inductive  coupling  to  the 
beam  electrons  is  mocked  up  by  a constant  current  boundary  condition  j 
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i(E  + 4ffJ  ■ 4ffJ0) . The  resultant  electron  heating  is  as  shown  on 
the  slide: 

(i)  an  initially  rapid  rise  of  Te  to  e drift  “ ^ mvd2  ‘ 
due  to  the  two  stream  interaction,  if  Te<e^ 
to  start  with,  followed  by 

(ii)  a progressively  slower  increase  in  Te  above  e drift 
due  to  ion-acoustic-like  turbulence. 

In  this  latter  phase,  the  temperature  only-increases  logarithmically 
in  time  during  the  simulation,  corresponding  to  an  exponential 
decrease  in  veff  with  temperature.  This  is  a faster  decrease  than 
the  linear  ion  acoustic  growth  rate  displays,  and  the  implication  is 
that  only  a few  units  of  return  current  drift  energy  may  be  dissipate* 
within  —10s  upe’1-  Consequently,  when  nB/rip  is  very  small,  very 
little  of  the  beam  energy  can  be  dissipated  by  these  return  current  1 
interactions,  within  any  reasonable  time.  | 

The  difficulty  with  heating  very  dense  targets  via  return  current 
dissipation  may  also  be  seen  easily  from  an  elementary  examination.  • 
of  the  macroscopic  current  decay,  as  illustrated  on  the  next  slide: 


It  is  clear  that  if  the  return  current  decay  time,  which 


increases  with  the  target  plasma  density,  should  become  larger  than 
the  beam  pulse  time,  then  the  return  current  will  be  inductively 
* quenched  prematurely  when  the  tail  of  the  beam  passes  by.  This 
prevents  completion  of  the  resistive  decay  which  heats  the  plasma. 

The  effective  collision  frequency  of  the  dissipative  microturbulence 
tends  to  fall  from  .2  wpi  to  .01  wpi  as  the  heating  proceeds,  so  that 
for  dense  targets  and  centimeter  beams , the  decay  time  can  well  be 
microseconds  in  many  cases  and  may  exceed  the  pulse  time.  The  energy 
deposition  efficiency  for  return  current  interactions  is  indicated 
on  the  bottom  of  the  slide  for  a relativistic  beam,  and  the  penalty 
for  large  decay  times  is  evident.  Since  the  transit  time  is  only 
nanoseconds,  the  ratio  of  "times"  will  be  small,  and  large  v/y  beam 
parameters  are  essential  to  increase  the  efficiency.  (Here 
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, and  G is  a dimensionless  geometrical  factor 


relating  to  flux  linkage.) 

In  summary,  dense  target  plasmas  can  only  be  heated  effectively 
by  return  current  interactions  if  the  beam  current  and  current  density 
are  increased  to  the  point  that  vB/yB  is  large,  Tdecay^pulse  *-8  not 
too  large,  and  iip/ng  is  not  too  large.  In  short,  one  must  confront 
the  difficult  problem  of  propagating  ultra  high  beam  current  densities 


to  heat  dense  target  plasmas  in  this  way. 
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An  alternative  mechanism  for  depositing  the  beam  energy  in  the 
target  plasma  is  the  primary  beam  plasma  interaction.  This  class  of 
primary  interactions  is  richly  endowed  with  a large  parameter  space, 
which  encompases  a hydrodynamical  and  a kinetic  regime,  longitudinal 
and  transverse  waves,  electrostatic  and  electromagnetic  regimes, 
gyroeffects,  relativistic  effects,  effects  of  finite  geometry  and 
inhomogeneity,  and  the  various  nonlinear  effects  which  appear  as 
the  interaction  evolves.  Without  .attempting  to  survey  all  of  this 
parameter  space,  I would  only  point  out  some  of  the  interesting 
scalings  of  a well  studied  and  important  subclass  of  primary  inter- 
actions, as  indicated  on  the  next  slide. 

•■•  ••  if 


This  is  the  hydrodynamical  interaction  of  a sufficiently  mono- 
energetic  relativistic  beam,  generating  longitudinal  electrostatic  . 
turbulence.  • We  examine  this  interaction  in  one  dimension  for  the 
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case  where  the  beam  is  weak:  YfinB/np<  0*1  and  where  gyroeffects, 
geometrical  and  inhomogeneity,  effects  are  negligible.  The  "weak 
beam"  regime  is  appropriate  for  the  problem  we  have  posed  of 
- attempting  to  heat  plasmas  several  orders  of  magnitude  denser  than 
the  electron  beams,  and  it  implies  that  the  plasma  response  will  at 
first  be  linear  and  nonrelativistic,  while  the  beam  responds  non- 
linearly  to  the  narrow  spectrum  "single  wave"  which  the  interaction  pro- 
duces. We  are  interested  here  in  the  deposition  efficiency  from  the 
beam,  achieved  by  the  time  of  the  first  nonlinear  saturation  of  the 
interaction — as  defined  on  the  slide — and  in  particular  in  how  this 
efficiency':  depends  on  YBf  nB»  an<*  Hp* 

These  results  are  obtained  from  a variety  of  computer  calculations 
performed  at  the  University  of  Texas,  including  (l)"single  wave 
calculations"  in  which  the  linear  background  response  is  treated 
analytically  and  a single  wave  is  legislated  and  the  nonlinear  beam 
response  is  computed,  (2)  one  dimensional,  particle-in-cell,  relativistic 
calculations  with  two  and  occasionally  three  (including  ions)  species, 

' and  (3)a  very  few,  recent N two-dimensional,  relativistic,  electro- 
magnetic, particle  calculations. 

The  physics  of  this  interaction  is  by  now  very  well  known,  having 
been  studied  computationally  in  many  laboratories,  and  experimentally 
by  Ken  Gentle  at  Texas — among  others — for  the  nonrelativistic  version 
of  this  interaction.  A very  pure  electric  wave  is  observed  to  grow  up, 
propagating  barely  slower  than  the  beam,  and  trapping  the  beam 
electrons  when  it  attains  sufficient  amplitude.  As  the  beam  electrons 
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rotate  in  the  troughs  of  the  wave,  they  are  slowed  slightly  and 
energy  is  transferred  to  the  field  oscillations  and  to  initially 
coherent  kinetic  sloshing  oscillations  in  the  background  plasma 
. electrons.  It  was  early  pointed  out  by  Fainberg  and  coworkers 
that  for  a relativistic  beam,  even  a slight  slowing  in  velocity 
would  allow  an  appreciable  fraction  of  the  beam  energy  to  be  deposited. 

(Possible  slides  3.1,  3.2, on  computer  phase  space, etc.) 


The  computer  calculations  reveal  the  following  scaling  of  the 
energy  deposition  efficiency  with  Yg  and  njj/n_: 


‘it.'Xr.h.Ysl V? 


• M \v  h -W**  WjXC; 


lO« 

■m 

w: 

jaWN 

(1)  One  conclusion  is  that  for  a fixed  value  of  n p/ng  below 
a few  hundred,  the  efficiency  declines  with  Yg,  while  for  larger 
values  of  np/nfi  there  is  a weak  peak  in  e(yg). 

(2)  Second,-  for  fixed  choices  of  yB,  the  efficiency  may  be 
optimized  with  respect  to  ng/np.  The  efficiency  curves  are  more 
peaked  for  high  Yg  than  for  low  Yg,  and  they  are  steeper  for  ng/iip 
below  optimum  than  they,  are  above  the  optimum  points  (see  "arrows" 
on  slide  4) . 

(3)  The  peakedness  of  the  efficiency  curves  is  substantially 
reduced  for  ng/np  above  optimum,  if  efficiency  is  defined  from  the 
mean  energy  deposited  during  the  first  cycle  of  trapping,  rather  than 
from  the  peak  energy  deposited  (see  dotted  curve  on  the  slide). 

(4)  Consequently,  the  optimum  deposition  efficiency  can  be 
approached  provided  only  that  YbC1^/1^)  1/s£%»  and  this  optimum  will 
be  in  the  range  of  10%-40%.  This  means  that  higher-y  beams  are 
appropriate  for  heating  plasmas  much  denser  than  the  beam.  There 

is  no  requirement  that  Vg/yg  be  large  for  efficient  energy  deposition, 
as  there  is  for  return  current  dissipation. 

Finally,  if  we  assume  that  Yg  and  ng/np  have  been  matched  for 
optimum  efficiency,  this  energy  deposition  efficiency  may  be  plotted 
against  np/ng  as  in  the  next  slide: 
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For  the  smaller  values  of  rip/nB  where  the  efficiency  is  highest, 

the  primary  interaction  is  supplemented  by  reverse  current  dissipation 

\ 

both  the  background  and  the  beam  response  is  highly  nonlinear,  and 
the  turbulent  spectrum  is  rather  broad  from  the  start.  The  chief 
conclusion  we  wish  to  draw  is  that  for  the  larger  values  of  rip/ng 
within  the  "single  wave  regime,"  the  deposition  efficiency — although 
lower— is  still  appreciable.  The  efficiency  curves  are  "flattening- 
out"  verses  np.  This  means  that  extremely  dense  plasmas  can  be 
heated  with  much  less  dense,  relativistic  beams,  if  this  interaction 
can  be  achieved. 
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The  primary  difficulty  is  to  assure  that  the  beam  is  sufficiently 

monoenergetic  that  the  interaction  is  at  least  marginally  hydro- 

dynamical  rather  than  kinetic  in  nature.  This  means  that  the  angular 

1 11 

scatter  of  momenta  must  not  be  too  high—roughly  that  ~K2(ng/nj 

Achieving  this  beam  quality  is  also  relatively  more  difficult  for 
dense  target  plasmas,  but  it  might  prove  easier  than  propagating 
megamperes/cm2 . The  point  is  that  there  is  an  alternative  way  to 
go  in  trying  to  heat  dense  fusion  plasmas  with  relativisitc  electron 
beams.  One  can  either  improve  the  quality  or  else  the  current  density 
of  the  beam.  Although  there  lias  not  been  a large  experimental  effort 
to  improve  beam  quality,  steps  such  as  introducing  guide  magnetic 
fields  into  the  diode  appear  very  promising. 

Although  this  is  the  chief  point  of  this  talk,  I should  touch 
briefly  on  a few  other  points  that  have  been  neglected: 

(1)  Regarding  gyroeffects — it  will  generally  be  technologically 
unfeasible  to  magnetize  (i.e.;  Oe>a)pe^  the  target  plasma  at  the 
densities  we  have  been  considering  (1014  to  101'  cm_s).  However, 
it  is  entirely  possible  to  magnetize  the  beam,  such  that  B2/8TT>njj(YB-l)m! 

. (2)  In  the  absence  of  a magnetic  field,  the  beam  plasma  interaction 
with  the  fastest  growth  rate  is  the  electrostatic  transverse  beam 
instability,  with  kj.>k„  * u>p/0c.  However,  when  Yg(nB/np)  '*<  1,  the 
transverse  interaction  saturates  without  creating  so  much  momentum 
scatter  that  the  longitudinal  interaction  is  cast  into  the  kinetic  regim 
It  is  also  possible  to  strengthen  the  guide  magnetic  field  to  the 
point  [(ne/a>p)  « 1/S]  that  the  transverse  interaction  is 
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suppressed  while  the  longitudinal  interaction  is  unaffected.  In 
any  case,  our  preliminary  two  dimensional  computer  calculations, 
which  allow  both  transverse  and  longitudinal  interactions  with  B=0, 

. show  adequate  energy  deposition  not  much  below  the  one  dimensional 
predictions.  (Similar  results  have  been  reported  from  NRL.) 

(3)  The  effect  of  collisions  between  the  target  plasma  electrons 
and  ions  is  rather  weak  on  the  hydro  dynamical  primary  beam-plasma 
interaction,  but  it  can  lead  to  suppression  of.  the  kinetic  beam  plasma 
interaction  (where  the  beam  momentum  scatter  is  too  great.)  This  is 
why  the  beam  quality  must  be  assured  for  efficient  heating.  The  > : 
target  collisions  will  also  lead  to  randomization  of  the  coherent 
oscillations  excited  in  the  plasma  electrons--generally  on  usee  time 
scales--and  wilL  lead  to  energy  equilibration  with  the  ions --generally 
on  msec  time  scales. 

(4)  Finally  we  have  not  discussed  the  later  time  nonlinear  effect 
such  as  parametric  instabilities  or  trapped  particle  instabilities 
which  can  serve  to  further  enhance  the  energy  deposition  efficiency, 
to  randomize  the  turbulence  (spread  spectrum)  , and  to  heat  the  ions 
more  quickly,  provided  that  there  is  time  for  these  interactions  to 
occur.  Since  even  the  earliest  interaction  provides  resonable  de- 
position efficiencies  when  the  necessary  beam  quality  is  achieved, 
the  later  interactions  will  be  "icing  on  the  cake." 


(Possible  conclusion  by  showing  some  of  the  slides  from  a two 
dimensional,  B-0,  Y-q”3,  nn/n.i-8  beam  plasma  rinteraction  with 
.e~25%) 
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APPENDIX  F 


LETTER  SUMMARY  OF  BEAM-PLASMA  HEATING  RESEARCH 


AND  COMPARISON  OF  RESULTS  WITH  THOSE  OF 


THODE,  ET  AL. 


(1)  Letter  from  Brendan  B.  Godfrey  to  James  R.  Thompson 

(2)  Letter  from  James  R.  Thompson  to  Richard  L.  Gullickson 

(3)  Letter  from  James  R.  Thompson  to  Brendan  B.  Godfrey 


(Appendix  F contains  14  pages) 
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Dr.  J.  R.  Thompson 
Austin  Research  Associates,  Inc. 
1901  Rutland  Drive 
Austin,  Texas  78757 

Dear  Bob: 


I have  been  reading  over  some  of  the  work  you  and  others  at  ARA 
have  done  on  relativistic  electron  beam-plasma  interactions.  The  research 
seems  very  interesting,  but  is  not  described  in  sufficient  detail  in  the 
reports  I have,  namely 


"Turbulent  Heating  of  Plasmas  by  High  Energy  Particle  Beams" 


I-ARA-7  6-U-80 / ARA-23  6 
I-ARA-7  5-U-139/ARA-195 
I-ARA-7 4-U-13 9 /ARA-160 
I-ARA-73-U-88/ARA-106 


(Aug.  76) 

(Aug.  75) 

(Aug.  74) 

(June  73)  (Ch  III  only) 


I would  certainly  appreciate  any  additional  material  you  can  send  me. 


Enclosed  are  several  reprints  by  Les  Thode  which  you  may  find  of 
interest. 


Brendan  B.  Godfrey 


BBG : new 
Enel:  a/s 
cc:  B.  B.  Godfrey 
(2)  ISD-5 
(2)  T-15  File 
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1901  Rutland  Drive 
Austin,  Texas  78758 
(512)  837-6623 


August  22,  1977 


Captain  Richard  L.  Gullickson 

Program  Manager 

Directorate  of. Physics 

Air  Force  Office  of  Scientific  Research 

Bolling  Air  Force  Base,  D.C.  20332 

Dear  Dick: 


A couple  of  months  ago,  while  I was  writing  up  some  of  the  material 
for  the  renewal  proposal  which  we  sent  you  on  June  17,  1977,  I received 
from  Brendan  Godfrey  a request  for  elaboration  on  some  of  the  old  beam- 
plasma  interaction  research  which  we  had  performed  and  described  in  our 
prior  AFOSR  reports.  Apparently,  a question  had  arisen  during  an  AFWL 
meeting  on  a proposed  beam  plasma  heating  experiment  about  how  our 
research  compared  with  that  of  Les  Thode,  et  al. 


I was  too  busy  to  respond  to  this  request  for  some  time,  but  I 
finally  telephoned  Brendan  on  August  1 and  we  discussed  the  research. 

I understand  that  he  had  already  submitted  a preliminary  comparison  of 

°U"  r\Search  W±th  LeS  Thode's  research,  and  that  any  urgency  associated 
with  the  question  had  largely  dissipated,  inasmuch  as  AFWL  made  a 
commitment  to  go  ahead  with  their  experiment.  Nevertheless,  I have 
taken  the  time  to  make  a few  detailed  comments  on  this  research  in  a 
letter  to  Brendan,  and  since  the  work  was  largely  AFOSR-sponsored,  I 
have  enclosed  a copy  of  the  letter  herewith,  for  your  information. 

Despite  some  possible  minor  differences,  which  may  be  largely 
semantic  in  nature,  the  main  conclusion  which  struck  me  as  I looked  over 
some  1975-76  publications  by  Thode,  et  al.  (which  Brendan  had  sent  me) 
was  the  near  total  agreement  with  my  own  conclusions  in  1973-74,  when  I 
was  last  active  in  this  area.  This  is  not  to  belittle  the  very  valuable 
research  which  has  been  accomplished  in  the  meanwhile,  including  a 
number  of  two-dimensional  particle  simulations  to  probe  transverse 
dynamical  effects  and  magnetic  field  effects.  Nevertheless,  my  early 
cone ius ions  that  very  high  density  plasmas  might  be  heated  with  reasonable 

tHrrteTXCl  7 h±8?  qU*llty  (low  “S"1"  matter),  high-y  (ultrarelativistic) 
electron  beams,  via  the  primary  beam-plasma  longitudinal  streaming 
interaction,  seem  to  be  in  agreement  with  Les  Thode's  conclusions.  I 
believe  that  our  prescriptions  for  the  optimum  beam  parameters  are  also 
in  quantitative  agreement,  although  I have  not  seen  any  papers  by  Les  on 
the  AFWL  experiment  in  question.  Assuming  this  to  be  the  case,  I am 
wholeheartedly  in  favor  of  the  experiment  being  performed. 


Capt.  R.  L.  Gullickson 
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I will  be  interested  to  hear  how  the  Air  Force  Summer  Study  came 
out,  and  what  conclusions  were  reached.  Perhaps  I’ll  give  you  a 
call  in  a week  or  so  to  discuss  it. 

Also,  I presume  that  our  AFOSR  proposal  is  still  on  schedule 
in  the  reviewing  process.  Let  me  know  if  I can  be  of  any  more  help 
there. 

Finally,  the  typing  is  now  about  complete  for  submission  of  our 
publication  on  limiting  currents.  I should  be  sending  you  a copy  in  a 
week  or  so. 


Best  regards, 


James  Robert  Thompson 


JRT:sjb 
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Dr.  Brendan  Godfrey 
Group  T-15,  MS-531  . 

University  of  California 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87545 
Dear  Brendan: 


Following  our  telephone  conversation  of  August  1 regarding  beam-plasma 
heating,  I have  secured  a copy  of  some  lecture  notes  which  were  compiled 
for  a July,  1974  presentation  by  Dr.  Drummond  at  Novosibirsk,  which  summa- 
rize our  understanding  of  the  most  important  features  of  the  problem  as  of 
that  time.  These  notes  are  enclosed  herewith  (they  are  somewhat  rough,  as  the 
presentation  was  only  verbal.) 

Inasmuch  as  I spent  the  better  part  of  five  years  of  my  professional  life 
(i.e.,  1968  - 1973)  working  on  this  problem — mostly  supported  by  AFOSR — 
plus  another  year  (i.e.,  1973  - 1974)  trying  (unsuccessfully)  to  scare  up  sup- 
port for  an  experiment  to  test  some  of  our  conclusions,  I thought  I would  take 
the  time  to  clarify  our  chief  results  somewhat. 


The  historical  sequence  of  "our"  own  research  is  that  we  first  applied  a 
concentrated  theoretical-computational-experimental  attack  on  the  nonrelativ- 
istic,  one-dimensional  (strong  magnetic  guide  field),  weak  beam-plasma  inter- 
action, and  we  found  a gratifying  amount  of  agreement  on  the  results. We 
subsequently  engaged  in  a period  of  protracted  theoretical  research  designed 
to  extend  understanding  to  more  and  more  effects:  relativistic  beams,  electro- 
magnetic effects,  transverse  waves,  kinetic  (or  angular  scattering)  effects, 
gyroeffects,  finite  geometry  effects,  ion  effects,  trapped  particle  effects, 
parametric  effects,  etc.  Simultaneously,  our  computer  simulation  capability 
(via  Adam  Drobot)  was  slowly  improved  to  handle  1-D,  relativistic,  single  wave 
calculations  (a  la  Matsiborko,  et  al.  ),  then  1-D,  relativistic  particle  simu- 
lations, and  finally  2-D,  relativistic,  electromagnetic  particle  simulations. 

At  the  culmination  of  this  research  effort  in  1973,  the  outstanding  facts  which 
struck  me  were: 
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(1)  Technological  developments  in  high  voltage  pulsed  power  seemed 
very  promising,  with  relativistic  electron  beams  already  capable 
of  delivering  tens  to  hundreds  of  kj.  The  future  growth  potential 
for  engineering  desirable  beam  properties,  increasing  the  beam 
energy,  etc.  seemed  very  bright. 

(2)  Despite  the  large  amount  of  attention  then  being  given  to  return 
current  interactions,  partly  because  of  their  relevance  to  a 
number  of  existing  experiments  with  np/ng  < 100,  these  inter- 
actions held  little  theoretical  promise  for  heating  plasmas  at 
much  larger  density  ratios. 

(3)  However,  based  only  on  the  energy  coupling  obtainable  from  the 
primary,  longitudinal,  beam  plasma  interaction  in  the  hydro- 
dynamical  regime,  it  already  seemed  possible  to  heat  extremely 
dense  target  plasmas  to  fusion  temperatures  if  the  beam  density, 
energy,  and  quality  (i.e.,  angular  scatter)  could  be  properly 
controlled. 

(4)  Although  the  transverse  waves  and  the  beam  cyclotron  modes  were 
not  yet  as  well  understood  as  the  longitudinal  plasma  wave  inter- 
action, it  seemed  that  they  could  only  enhance  the  energy  coupling 
and  should  not  quench  the  longitudinal  interaction.  Even  without 
this  complete  understanding  of  the  transverse  interactions,  the 
preliminary  evidence  (discussed  later)  was  already  adequate  to 
justify  designing  an  experiment  to  heat  with  a high  quality,  highly 
relativistic  beam,  via  the  primary  longitudinal  interaction. 

(5)  Nevertheless,  the  thrust  of  the  existing  experimental  effort  seemed 
to  be  in  the  wrong  direction,  from  the  point  of  view  of  the  appli- 
cation of  single  pass  heating  of  very  dense  (i.e.,  > 101S  electrons/ 
cm3)  target  plasmas,  which  I viewed  as  one  of  the  most  promising 
fusion  goals  to  which  the  e-beam  technology  could  be  applied,  as 
well  as  an  area  where  it  should  have  a distinct  competitive  advan- 
tage. That  is,  the  existing  effort  seemed  aimed  at  producing  higher 
v/y  beams — giving  up  y for  v if  necessary,  and  producing  high  current 
density  beams.  There  was  almost  no  effort  to  produce  higher  - y 
beams,  or  to  improve  (or  even  measure)  the  beam  quality.  (In  fairness, 
I should  concede  that  a justifiable  motivation  for  some  of  the  emphasis 
on  high  current  density — small  cross  section — beams  was  the  pellet 
fusion  problem.) 

(6)  This  lack  of  interest  in  high  - y,  high  quality  (low  9)  beams  existed 
despite  the  likelihood  that  the  poor  coupling  efficiencies  observed 
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• occasionally  for  very  dense  target  plasmas  were  probably  due  to 

interactions  occurring  in  the  kinetic  regime  where  they  were  sub- 

Iject  to  collisional  quenching.  In  short,  there  was  not  much  appre- 
ciation for  the  dividends  which  would  accrue  from  improved  beam 
quality. 
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(7)  Finally,  again  from  the  point  of  view  of  relevance  to  a high  density 
heating  experiment,  it  seemed  to  me  that  there  was  a bit  too  much 
theoretical  attention  devoted  to  the  slower,  late-time,  nonlinear 
homogeneous  theory — and  perhaps  not  enough  attention  devoted  to 
effects  of  inhomogeneities,  finite  pulse  time,  and  finite  target  size. 

It  is  the  above  points  which  I tried  to  make  in  my  AFOSR  reports  covering  the 
1972  - 1974  research  period,  and  particularly  in  my  talk  at  the  Working  Group  on  the 
Interaction  Between  Plasmas  and  Intense  Laser  and  Electron  Beams,  at  the  International 
Centre  for  Theoretical  Physics,  Trieste,  Italy  in  August  1973,  and  in  my  papei?  at 
the  1973  Philadelphia  APS  meeting.  This  paper  in  turn  was  the  basis  of  the  enclosed 
notes  which  were  written  up  for  Dr.  Drummond's  1974  Novosibirsk  lecture.  If  this 
material  has  left  the  impression  that  we  were  only  concerned  with  one-dimensional 
effects,  or  that  we  had  not  done  any  work  on  transverse  interactions  or  gyroeffects, 
that  was  unintentional.  It  is  just  that  there  were  larger  points  which  I was 
trying  to  make  without  getting  bogged  down  with  subtleties,  and  the  computer 
evidence  on  transverse  interactions  was  still  too  preliminary  for  any  definitive 
statements  to  be  made. 

During  1973  - 1974,  I tried  to  promote  interest  in  a beam  heating  experiment 
of  a dense  target  plasma  with  the  high  - y beam  at  the  Aurora  facility,  Harry 
Diamond  Laboratories,  with  the  payoff  to  be  a plasma  radiation  source.  I first 
had  a long  battle  to  convince  them  that  collisional  stabilization  of  the  beam- 
plasma  interaction  could  be  avoided  if  the  beam  quality  could  be  controlled. 

Then  we  performed  a study  for  them  which  included  a rough  experimental  design. 

We  proposed  to  arrange  yB  (ng/iu,)1^3-  0.5,  with  a strong  magnetic  guide  field 
extended  into  the  diode  to  give  us  a better  chance  of  holding  the  angular  scat- 
ter to  0^3  Yg1  = 2(nB/np)1/3.  This  experimental  proposal  ultimately  floun- 
dered, not  because  of  skepticism  regarding  beam  heating,  but  because  the  expected 
pulse  duration  of  the  emitted  radiation  was  somewhat  longer  than  desired. 

Next,  let  me  touch  briefly  on  my  1973  understanding  of  some  of  the  trans- 
verse effects.  First  of  all,  for  an  infinite  magnetic  guide  field,  the  hydro- 
dynamical,  electromagnetic  dispersion  relation  is 
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and  the  most  unstable  wave  occurs  for  kx  = 0 (when  allowed  by  infinite  transverse 
boundaries)  and  is  electrostatic.  However,  in  a waveguide,  the  minimum  allowable 
kx ^ (2.4/a),  and  stability  can  result  if  the  wall  radius  is  small  enough  that 


this  exceeds 


At  the  other  extreme,  when  there  is  no  magnetic  field  the  hydrodynamical, 
electromagnetic  dispersion  relation  is 


Yb3(oj  - kz6Bc) 


<l)  2 + 


+ k2c2  - w2  - 


k12SB2c2o>B2ulp2 
Ybu)2(o)  -kz6Bc)2 


This  is  always  unstable,  and  includes  a longitudinal  and  a transverse  regime, 
defined  according  to  whether  kx  is  small  or  large.  The  longitudinal  regime  is 
electrostatic,  while  the  transverse  regime  is  electrostatic  for  large  kz  and 
electromagnetic  for  small  kz.  The  case  of  E0  = 0,  kx  -*■  »,  kz  -*■  0 is  sometimes 
referred  to  as  the  Weibel  instability.  Simulations  showed  that  rather  than 
transferring  much  beam  energy,  it  tended  to  cause  the  beam  to  filiment  in  the 
transverse  plane.  The  fastest  growth  for  B0  = 0 occurs  for  kz  = tUp/Bgc,  where 
the  mode  is  electrostatic,  and  the  growth  rate  of  this  streaming  interaction  is 
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I usually  define  the  longitudinal/transverse  wave  regimes  via  this  expression. 

A transition  region  extends  from  kx  ■ kz/Yg  to  kx  ■ kz,  with  kx  < kz/Yg  being 
longitudinal  and  kx  > kz  being  transverse,  and  the  middle  region  being  a mixture. 


Therefore,  for  an  optimized  interaction  with  yb  (nB/np)1/3  - 0.5,  I would 
expect  the  longitudinal  wave  regime  to  be  characterized  by  kz  - u)p/c  and 
0 < kx<^  kz/YB  - 2(ng/np)1,/3  k?.  Therefore,  the  "typical"  case  that  you  mentioned 
of  < 20°  angles  in  k-space  could  be  within  what  I term  the  longitudinal  regime 
unless  Yg  >>  cot  20°  = 2.75,  while  the  strongly  transverse  regime  corresponds  to 
angles  > 45°.  Clearly,  the  large  grey  area  invites  semantic  debate.  In  fact,  I 
now  suspect  that  some  of  our  apparent  differences  during  our  phone  conversation 
were  merely  semantic,  since  I find  the  above  mentioned  range  of  kx  cited  in  the 
very  nice  series  of  1975-76  papers  8-13  which  you  sent  me.  Of  course,  in  a 
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waveguide  it  is  theoretically  possibly  for  the  lowest  allowable  value  of  kx  ^ (2.4/a) 

k_  ^p 

to  fall  outside  of  — = — -r- — , thus  eliminating  the  longitudinal  regime,  but  for 

YB  YB^Bc 

the  problem  of  heating  very  dense  plasmas  this  will  not  occur.  In  any  case,  it 
has  always  been  clear  that  the  eigenmode  structure  for  a longitudinal  interaction 
will  necessarily  have  a non-zero  kx  component,  so  that  two  dimensional  computer 
simulations  are  necessary  for  a good  mock-up  of  the  physics.  Nevertheless,  one 
expects  that  the  basic  processes  observed  in  the  1-D  simulations  (i.e.,  longitudinal 
trapping,  etc.)  would  continue  to  dominate  the  longitudinal  regime  of  2-D  simula- 
tions, and  we  and  others  confirmed  this  to  be  the  case. 

A more  crucial  question  is  what  happens  to  the  strongly  transverse  waves  with 
kx  > k||  , whose  growth  rate  exceeds  that  of  the  longitudinal  waves  by  Y^3.  One 
might  expect  these  modes  to  saturate  at  relatively  low  amplitudes  due  to  transverse 
velocity  scatter  since  kx  Is  large,  and  indeed,  the  Soviets  (particularly  Rudakov) 
made  calculations  of  this.  Then  the  question  is  whether  the  associated  axial 
velocity  spread  due  to  this  scattering  would  throw  the  coexisting  longitudinal 
interaction  into  the  kinetic  regime.  The  answer  was  that  for  yb  (nfi/np)1/3  < 1, 
it  would  not — the  longitudinal  interaction  would  remain  hydrodynamical . My  own 
calculations  tended  to  confirm  these  results,  although  admittedly,  when 
Yb  ,(ni)/np)1/3  = 0*5  f°r  an  optimum  interaction,  one  might  expect  the  hydrodynamic 
criterion  to  remain  only  marginally  obeyed.  My  experience  suggested  to  me  that 
this  would  probably  be  good  enough,  although  the  proof  of  the  pudding  is  in  the 
2-D  simulations.  Our  own  2-D,  B0  - 0 simulations  had  just  become  available  in 
1973,  and  the  preliminary  results  supported  the  above  picture  of  the  interaction, 
with  continued  longitudinal  trapping,  and  with  field  amplitudes  comparable  to  the 
1-D  prediction  at  the  time  of  initial  saturation  of  growth  of  the  field  amplitudes. 
Lampe,  et  al.  reported  similar  preliminary  results  at  the  same  1973  APS  meeting. 

This  was  essentially  the  state  of  affairs  when  I ceased  active  research  on  this 
problem.  During  our  telephone  conversation,  I got  the  Impression  that  perhaps 
subsequent  2-D  simulations  had  led  you  to  disavow  the  above  described  picture. 
However,  in  looking  over  the  references'-*1^  that  you  sent,  I find  general  agree- 
ment and  even  affirmation  of  these  arguments  in  several  places.  Perhaps  you 
could  point  out  whether  and  where  any  disagreement  remains. 

Next,  consider  the  questions  of  gyroeffects.  For  this  problem  of  heating 
a very  dense  plasma  with  rip  > 104ng,  technological  limitations  will  insure  that 
fte  « (u_  for  the  target  plasma,  so  the  magnetic  corrections  may  be  neglected  for 
the  streaming  interactions  with  the  target  plasma  modes.  However,  the  magnetic 
field  will  allow  transverse  beam  cyclotron  modes  as  well  as  the  longitudinal  beam 
plasma  modes.  We  found  that  the  kx  -*■  kz  -►  0 Weibel  instability  may  be  stabi- 
lized if  S)e  > BgYg^ 2 (nj}/np)1<tt(1>p  • Th*  more  interesting  case  is  the  streaming 
Interaction,  where  peak  growth  occurs  in  the  electrostatic  regime,  for  k2  = (Up/6gC. 
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• The  electrostatic  dispersion  relation,  for  | oo | » fi  , is 
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which  yields  the  previously  quoted  growth  rate  when  Bc 


It  may  be  seen  that  now 


beam-cyclotron  modes  coexist  with  the  beam  plasma  modes,  and  each  may  interact  with 
the  background  plasma  mode.  The  beam-cyclotron  interaction  occupies  the  large-kx, 
transverse  region  of  phase  space,  where  the  beam  gyroeffects  modify  the  interaction 
from  the  above  discussed  case  for  B0  * 0) . When  (8e/u)p)  « 0.5  Y^Ang/np)2^3, 
the  gyroeffects  only  cause  a slight  slowing  of  the  usual  transverse  interaction,  but 
when  0.5  y^3  (nB/np)L/3«  fte/Wp  « 1,  the  transverse  growth  rate  is  slowed  significantly 

to 
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for  kx  » kz  = u)p/8Bc.  (Here  I should  apologize  for  my  use  of  the  word  "suppressed" 
on  page  13  of  the  Novosibirsk  notes.  This  is  misleading,  since  the  transverse 
interaction  is  not  stabilized,  but  is  only  slowed  by  the  gyroeffects,  and  it  will 
still  be  faster  than  the  usual  longitudinal  interaction  unless  the  magnetic  field 
is  raised  by  an  additional  factor  of  yB/3. ) The  more  relevant  question  is  again  the 
nonlinear  saturation  level  of  the  transverse  beam-cyclotron  interaction,  and 
whether  or  not  it  throws  the  coexisting  longitudinal  interaction  into  the  kinetic 
regime.  Our  estimation  was  that  the  transverse  wave  saturation  would  be  at  or 
below  the  level  occuring  for  B0  * 0,  since  additional  saturation  mechanisms  (e.g., 
finite  kp  - effects)  are  now  available.  Therefore  we  expected  the  longitudinal 
interaction  to  remain  hydrodynamical  under  the  previous  condition  that  YB  (nB/np)1/;^ 

It  was  not  clear  to  us  in  1973  whether  the  modification  in  the  interactions  due 
to  B0  + 0 would  have  a net  effect  of  increasing  or  decreasing  the  energy  coupling 
effeclency  from  the  beam,  although  we  expected  the  change  in  either  direction  to 
be  modest-to-small.  (Any  loss  in  transverse  Interaction  coupling  might  be  made  up 
in  a more  efficient  longitudinal  interaction.)  In  any  case,  we  again  deferred  to 
the  computer  simulation  evidence,  which  consisted  for  us  of  only  two  or  three 
2-D,  B0  + 0 runs  at  that  time.  These  indeed  showed  approximately  the  expected 
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. amount  of  energy  transfer  by  the  Initial  nonlinear  interaction,  and  continued 
to  show  the  trapping  dynamics  associated  with  the  longitudinal  interaction.  My 
own  feeling  was  that  the  importance  of  the  magnetic  field  was  not  likely  to  be 
associated  so  much  with  the  way  the  cyclotron  modes  directly  influenced  the  energy 
transfer,  as  with  the  way  the  field  might  improve  the  beam  quality  if  it  were 
carefully  designed  in  the  diode  and  injection  regions. 

During  our  telephone  conversation,  I got  the  impression  that  your  conclusions 
regarding  the  magnetic  field  effects  might  differ  from  ours  described  above  in 
two  respects.  First,  I thought  that  you  were  claiming  that  in  a parameter  range 
recommended  for  a heating  experiment,  the  transverse  beam-cyclotron  interaction 
itself  was  observed  to  stabllze  via  beam  trapping.  However,  after  looking  over 
your  papers,  I suspect  that  this  was  a misunderstanding  and  that  you  might  agree 
with  us  that  it  is  the  coexisting  longitudinal  beam-plasma  interaction  with  the 
background  plasma  mode  which  saturates  by  beam  trapping.  We  have  found  that  there 
are  almost  always  other  avenues  for  nonlinear  saturation  of  the  beam-cyclotron 
interaction  which  will  enter  prior  to  beam  trapping. 

The  second  possible  difference  is  that  I thought  I understood  you  to  say 
that  the  magnetic  field  had  a very  strong  influence  on  the  energy  coupled  out  of 
the  beam,  and  that  there  was  an  optimum  magnetic  field  strength  choice  which  would 
allow  significant  gains  in  the  energy  coupling  efficiency  over  a non-optimum  choice. 
However,  in  looking  over  one  of  your  1976  references1  , you  seem  to  agree  with  our 
own  preliminary  conclusion  that  the  beam  energy  loss  is  only  weakly  dependent  on 
the  magnitude  of  the  magnetic  field,  if  < w_.  Perhaps  you  could  clarify  whether 
subsequent  simulations  led  you  to  change  the  view  expressed  in  this  paper,  or 
whether  there  might  again  be  semantic  confusion  over  the  adjective  (strong,  weak) 
attached  to  your  observations. 

Before  leaving  the  topic  of  streaming  interactions  which  involve  magnetic 
field  modes,  let  me  mention  two  other  areas  not  discussed  above.  First,  when 
&e  <<  u , there  will  be  Whistler  and  electron  cyclotron  modes  in  the  background 
plasma,  in  addition  to  the  background  plasma  mode  mentioned  above,  which  are 
candidates  for  interaction  with  the  beam-plasma  and  beam-cyclotron  modes.  The 
Interaction  between  these  and  the  beam-plasma  mode  will  generally  be  stable  for 
a high  velocity  beam,  since  the  beam  mode  will  then  pass  above  the  knee  of  the 
Whistler  mode.  However,  the  lower  beam-cyclotron  mode  may  be  expected  to  Intersect 
the  Whistler  at  low  frequency  and  long  wavelength,  where  ion  effects  and  electro- 
magnetic effects  are  important.  However,  we  tended  to  dismiss  this  interaction 
on  the  grounds  of  its  being  slower  and  probably  less  potent  than  the  longitudinal 
beam-plasma/target-plasma  mode  Interaction.  We  never  detected  it  in  our  simulations 
I wonder  if  you  have? 
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The  second  area  is  the  regime  » u)p,  which  I have  not  discussed  in  detail 
on  the  grounds  of  technological  infeasibility,  for  target  plasma  densities  above 
1015  cm"3,  or  so.  However,  in  looking  over  your  reference  on  linear  theory®,  1 
noticed  that  there  appeared  to  be  an  error  in  the  analytical  expression  presented 
for  the  frequency  and  growth  rate  of  the  interaction  between  beam-plasma  modes 
and  background  cyclotron  modes  in  this  parameter  regime.  This  is  the  portion  of 
your  table  I which  contains  in  your  Region  3,  and  the  growth  rate  expression 

appears  to  be  illustrated  by  curves  drawn  on  your  Figures  13  and  14.  The  problem 
is  that  your  expression  for  the  growth  rate  is  independent  of  the  magnetic  field 
strength,  whereas  the  growth  should  clearly  vanish  as  -*■  <®.  The  relevant  terms 
in  the  electrostatic  dispersion  relation  are 

a)  2sln  20  u 2 cos2 9 

0 - 1 E S 

o)2  - n2  yb3(w  - kz®Bc)2 
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where  tan  0 - kj./k  . In  the  regime  1 « — « — — - — 
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the  peak  growth  rate  is  approximately  given  by 
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where  k26Bc  = (fte2  + ou2  dn  2 0)1/a.  This  growth  rate  is  smaller  than  your  expression, 
and  is  nolonger  symmetric  about  (^45®  in  qualitative  agreement  with  your  plotted 
computer  results.  Other  small  corrections  and  angular  dependence  may  be  harvested 
if  the  small  neglected  terms  in  the  electrostatic  dispersion  relation  are  reintro- 
duced, or  if  electromagnetic  corrections  are  recovered.  However,  this  formula 
shows  the  major  effect. 

If  the  magnetic  field  is  strengthened  further,  to  the  point 
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then  (because  of  the  Increased  relative  significance  of  the  splitting  of  the 
beam  plasma  modes)  the  Interaction  transits  from  the  cubic  regime  to  a slower 
quadratic  regime,  where 
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and  peak  growth  occurs  for  k $gC  = (fte2  + u)p2  sin2 
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This  growth  rate  falls  off  with  increasing  BQ  somewhat  more  rapidly  than  in  the 
cubic  regime.  (Of  course,  there  will  still  exist  the  interaction  between  the  beam 
plasma  mode  and  the  background  plasma  mode,  with  peak  growth  for  0 ■ 0 but  finite 
growth  for  large  0,  which  will  compete  with  this  cyclotron  mode  interaction.) 
Unfortunately,  you  did  not  run  fle/o)p  high  enough  in  your  calculations  for  the  paper 
that  this  flaw  in  your  formula  would  have  caught  your  attention. 


I only  point  this  correction  out  for  academic  reasons,  since  I believe  this 
strong  - B0  regime  to  be  irrelevant  for  the  beam-heating  problem  of  Interest. 

(While  I am  at  it,  I presume  you  have  noticed  the  1/2  -*•  1/3  typo  in  the  next-to- 
last  paragraph  of  page  346.) 

I had  thought  of  elaborating  on  kinetic  effects  somewhat,  but  I think  I will 
pass  over  that  since  this  letter  is  already  so  long  and  we  seem  to  be  pretty  much 
in  agreement  there  anyway.  Let  me  just  wholeheartedly  endorse  any  experimental 
efforts  to  improve  the  beam  quality. 

One  final  area  where  we  might  have  some  small  quantitative  differences  is  in 
the  energy  coupling  efficiency  to  be  expected.  Back  in  1973-74,  there  was  fairly 
general  agreement  that  efficiencies  in  the  range  of  20%  - 30%  were  reasonable 
expectations  for  a wide  range  of  beam  energy  and  density  ratios  above  ^ 10,  and 
allowing  some  late-time  nonlinear  Interactions — but  not  carrying  this  to  extremes. 
Some  of  our  own  predictions  are  depicted  graphically  in  the  Novosibirsk  notes, 
but  the  graphs  should  be  interpreted  with  some  care.  I am  intrigued  by  your 
statement  over  the  telephone  that  'V  60%  coupling  efficiency  may  be  expected  in 
the  optimum  case  for  target  plasma  heating  at  very  high  density  ratios,  like 
iip/ng  ^ 101* . I did  not  notice  any  coupling  efficiencies  that  high  in  your  refer- 
ences, and  I would  have  guessed  that  it  could  only  be  achieved  by  going  to  low 
density  ratios  of  ^ 10  or  so,  or  by  allowing  the  beam  to  coexist  with  the  target  for 
an  extraordinary  length  of  time  (or  space) . I would  be  curious  to  hear  if  perhaps 
there,  is  some  recent  evidence  which  gives  you  this  confidence,  and  if  the  evidence 
is  simulations,  whether  any  have  shown  this  result  at  the  very  high  density  ratios 
of  interest. 
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. Before  closing,  let  me  say  that  in  giving  all  this  discussion  of  my  under- 

standing of  some  aspects  of  this  problem,  I certainly  did  not  mean  to  imply  that 
I was  telling  you  anything  new,  since  I am  well  aware  that  you  and  many  others 
have  been  down  this  same  path.  Rather,  I sought  to  fill  some  of  the  gap  of 
elaboration  that  you  perceived  in  some  of  our  earlier  AFOSR  reports.  In  fact, 

I claim  no  originality  for  any  but  the  earliest  non-relativistic  work  that  I did, 
since  the  modus  operandi  for  beam  plasma  work  a few  years  ago  was  for  - ten  groups 
around  the  country  to  co-discover  things  more  or  less  independently. 

Let  me  also  say  that  in  dwelling  on  some  areas  of  possible  small  differences 
in  our  work,  I have  probably  misrepresented  the  dominant  conclusion  that  struck  me 
in  comparing  our  work — which  is  the  extraordinary  extent  to  which  we  are  in  agree- 
ment. I certainly  think  that  your  careful  exploration  of  two-dimensional  effects 
and  the  many  two-dimensional  simulations  were  well  worthwhile  and  represented  a 
major  advance  in  the  field.  I was  particularly  impressed  with  Les'  paper  surveying 
the  experimental  results.  It  seems  to  me  that  the  evidence  supporting  the  type 
of  high  - Y»  high  quality,  beam  heating  experiment  of  very  dense  plasma  targets 
which  I fought  for  in  vain  is  now  well  nigh  overwhelming,  and  if  I can  endorse 
one  of  your  proposals  for  such  an  experiment,  I will  be  only  too  Happy  to  do  so. 

If  you  do  become  actively  involved  in  pursuing  such  an  experiment,  I would 
be  Interested  to  receive  any  reports  which  might  be  issued  on  it.  I presently  have 
no  precise  idea  of  the  parameters  which  you  have  in  mind,  the  experimental  motiva- 
tion (e.g.,  radiation  source,  fusion),  and  so  forth. 


Best  regards. 


James  R.  Thompson 


JRT:tb 

Enclosure:  Report  I-ARA-74-U-114 
cc:  Capt.  Richard  Gullickson 


Dr.  Adam  Drobot 
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APPENDIX  G 


SELF  CONSISTENT  EQUILIBRIA  OF  UNNEUTRALIZED 
RELATIVISTIC  ELECTRON  BEAMS 


In  order  to  be  able  to  exploit  the  powerful  electron 
beams  which  technology  has  made  available  for  a number  of 
interesting  applications,  such  as  collective  ion  accelera- 
tion, it  is  first  necessary  to  understand  the  equilibrium 
properties  of  the  electron  beam  itself.  The  most  fundamental 
question  is  to  determine  the  conditions  under  which  an  equi- 
librium of  forces  can  be  arranged  such  that  the  electron 
beam  may  propagate  in  vacuum  away  from  the  diode  into  a drift 
tube  in  a grossly  stable  configuration.  If  such  a force 
equilibrium  is  not  achieved,  the  unneutralized  electron  beam 
will  rapidly  blow  itself  apart  by  the  powerful  self  electro- 
magnetic fields  generated  by  its  own  space  charge.  A second 
objective  of  such  an  equilibrium  analysis  is  to  determine 
the  self  consistent  radial  profiles  of  all  equilibrium  quanti- 
ties, including  the  electromagnetic  fields  and  the  number 
density  and  flow  velocity  of  the  electron  beam.  These  equi- 
librium properties  will  in  turn  be  needed  to  determine  the 
linearized  eigenmodes  (or  waves)  which  may  be  supported  upon 
the  electron  beam. 

Although  one  important  application  for  this  analysis 
is  collective  ion  acceleration,  the  contribution  of  the 
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subject  ions  (to  be  accelerated)  to  this  equilibrium  is 
neglected  in  this  calculation.  The  first  reason  for  this 
is  that  for  many  accelerator  configurations,  the  electron 
beam  must  propagate  stably  for  some  distance  in  vacuum 
prior  to  the  introduction  of  ions,  so  that  an  analysis  of 
pure  electron  beam  equilibria  is  required.  Second,  the 
accelerator  can  operate  in  principal  in  the  "test  ion" 
regime,  where  the  ion  contributions  to  the  equilibria  and 
waves  are  truly  negligible.  Finally,  it  is  anticipated 
that  the  highest  density  of  ions  which  may  be  accelerated 
will  still  be  much  less  dense  than  the  electrons,  so  that 
the  ion  effects  may  be  treated  as  a perturbation  on  the 
fields  supported  by  the  electron  beam. 

The  electron  equilibrium  may  be  examined  by  con- 
sidering the  single  particle  equations  of  motion  together 
with  the  Maxwell  equations  for  the  self-consistent  electro- 
magnetic fields. 
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V*B=0  = VxE  + — 7 — B 
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V • E = - 4iTen 
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Here  e = - | e | , m are  the  electron  charge  and  mass,  and 
£c,  (y-l)mc2  are  its  velocity  and  kinetic  energy.  Since 
the  single  particle  equations  of  motion  (1)  are  the  char- 
acteristics of  the  collisionless  Vlasov  equation,  the  equi- 
librium distribution  function  f (x,  B.)  may  be  constructed 
from  any  constants  of  the  single  particle  motion.  We  shall 
adopt  a cold  fluid  irodel  of  the  electron  beam,  in  which  the 
beam  pressure  is  neglected,  and  n,  B_  are  the  electron 
number  density  and  flow  velocity,  which  may  in  general  be 
functions  of  x,  t. 

The  existence  of  steady  state  beam  equilibrium 
solutions  is  by  no  means  guaranteed.  The  crucial  condition 


is  the  balance  of  radial  forces,  which  may  be  expressed  as 
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in  the  cylindrical  coordinates  r,  <p , z.  The  basic  assump- 
tions are  that  the  equilibrium  is  steady  state  (3/3t  = 0) 
and  cylindrically  symmetric  ( 3/ 3 4>  = 0)  throughout  the 
vacuum  waveguide.  It  is  further  assumed  that  there  is  a 
downstream  region  away  from  the  beam  diode,  fcr  which  the 
equilibrium  is  approximately  independent  of  axial  position 
(3/3z  = 0).  Practically  speaking,  this  will  be  true  even 
throughout  an  ion  acceleration  section  of  the  waveguide, 
wherein  the  waves  and  ions  are  accelerated,  if  the  z-varia- 
tion  in  equilibrium  quantities  is  small  compared  to  the 
r-variation,  and  if  the  z-variation  in  equilibrium  quantities 
is  small  (or  adiabatic)  compared  to  the  z-variation  in  wave 
quantities,  characterized  by  the  axial  wavelength. 

The  fields  E and  B are  the  self-fields  of  the 
r <(> 

electron  beam,  determined  in  the  downstream  equilibrium 
region  from  the  equations 
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It  may  be  seen  that  both  and  are  negative,  and 

that  E is  larger  than  B in  magnitude.  Consequently 
r <f> 

there  can  be  no  equilibrium  (in  vacuum)  for  which  8 = 0, 

<P 

since  the  sum  of  the  first  two  terms  of  equation  (5)  are 

always  positive,  evidencing  a net  outward  force  on  the 

electrons.  This  net  outward  force  must  be  overcome  by  the 

- 8X  B pinch  force,  without  the  final  centrifugal  force 
<p  z 

destroying  the  balance.  Here  B^  is  presumed  to  be  an 
externally  generated  magnetic  guide  field,  which  is  slightly 
reduced  by  the  beam's  diamagnetic  precession  in  $ accord- 
ing to 
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If  the  fields  E , B , B are  presumed  to  be  known, 

r <p  z 

equation  (5)  may  be  examined  for  the  conditions  which  these 

fields  must  obey  in  order  that  equilibria  exist  for  which 
d 

^r  = 0 = dt  ^r*  Th^s  analysis  is  most  easily  accomplished 

in  an  inertial  frame  of  reference  for  which  locally 

vanishes  (notice  that  in  general  3 will  vary  with  r) . 

z 

In  such  a "beam"  frame,  the  equilibrium  condition  (5)  may 
be  written  as 
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As  a function  of  B.,  the  left-hand  side  of  equation  (9) 

9 

increases  monotonically  from  zero  and  diverges  as  6^ 

approaches  one,  while  the  right-hand  side  is  a linear 

function  of  BA  (if  the  diamagnetic  B^ -dependence  of  B 
<p  <p  z 

is  ignored)  with  a positive  slope  and  negative  intercept 

at  B.  = 0.  There  will  obviously  be  two  solutions  if  the 
<P 

functions  intersect  (i.e.,  two  distinct  equilibria)  and 
no  solutions  if  they  do  not. 
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A detailed  algebraic  analysis  of  equation  (9)  reveals 
the  following  necessary  and  sufficient  condition  for  the 
existence  of  equilibria: 


where 
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and  the  function  Q(x)  is  sketched  in  Figure  (1)  below. 


Figure  1.  Function  Q(x) 


Here  E and  n are  the  radial  electric  field  and  the 
rB  B 

electron  density  as  seen  from  the  local  beam  frame,  and 

B (r)  is  the  total  axial  magnetic  field — including 
z 

diamagnetic  corrections. 

Equation  (10)  immediately  implies  two  necessary 

conditions  for  equilibrium.  First  is  that  the  argument  of 

Q,  - E /B  must  be  less  than  one,  which  is  obvious  from 
rB  z 

equations  (9)  or  (5)  since  otherwise  the  net  forces  are 
manifestly  positive.  Second,  since  Q £ 1,  the  coefficient 
of  Q in  equation  (10)  must  be  less  than  one.  This  second 
condition  ensures  that  the  centrifugal  forces  do  not  destroy 
the  equilibrium. 
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In  order  to  express  these  two  necessary  criteria 


in  terms  of  quantities  observed  in  the  laboratory  frame  of 
reference,  it  is  necessary  to  make  a Lorentz  transforma- 
tion, with  the  transformation  velocity  c 6 (r) . The 

z 

results  are 
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where  Y_(r)  = [1  - 82(r)]-1/2.  Both  conditions  (11)  and 
i z 

(12)  must  be  obeyed  for  all  r.  Since  the  numerator  of 
condition  (11)  increases  roughly  as  r,  while  the  denomina- 
tor of  (11)  increases  only  weakly  with  r,  condition  (11) 
will  be  most  stringent  near  the  beam  edge  r = a.  The 
numerator  of  condition  (12)  is  roughly  constant,  while  the 
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denominator  increases  weakly  with  r,  so  that  condition 
(12)  will  be  most  stringent  near  the  beam  axis  r = 0. 

For  virturally  all  cases  of  interest,  both  condi- 
tions (11)  and  (12)  will  be  well  satisfied  and  the  equi- 
libria may  be  readily  analyzed  in  terms  of  the  corresponding 
small  parameters.  In  this  case,  the  two  equilibrium  solu- 
tions for  8^  will  be  well  separated;  the  "slow"  equilib- 
rium (for  which  8,  is  lowest)  is  the  one  desired  for  the 

<P 

accelerator  applications. 

In  order  to  self-consistently  compute  the  various 
equilibrium  quantities,  assumptions  about  the  nature  of  the 
beam  equilibrium  are  required,  since  many  different  equi- 
libria are  possible  in  general.  It  is  most  accurate  to 
make  these  assumptions  about  the  conditions  of  creation 
of  the  beam  at  the  cathode,  with  conservation  equations 
then  utilized  to  complete  determination  of  the  equilibrium 
quantities  in  the  downstream,  z-independent  region  of  the 
waveguide.  In  this  way,  the  accessibility  of  the  equilibrium 
is  included  in  the  calculation,  which  would  not  be  the  case 
if  ad  hoc  assumptions  were  instead  made  about  the  downstream 
equilibrium. 
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Therefore,  we  shall  assume  that  the  electron  beam 

is  born  at  the  cathode  with  B^c  = 0,  yc  = 1,  and  uniform 

current  density  i = - ecn  8 throughout  the  initial 

C o oz 

beam  radius  ' a . The  cathode  is  assumed  to  be  at  a uniform, 
o 

negative  electric  potential  4>c  < 0.  The  beam  is  assumed 
to  flow  from  the  cathode  through  a diode  acceleration  region 
where  it  is  raised  to  ground  potential  $ = 0 and  then 
injected  into  a cylindrical  waveguide  whose  conducting  walls 
are  also  at  ground  potential. 

The  initial  beam  properties  which  are  assumed  above 
may  be  related  to  the  corresponding  beam  properties  in  the 
downstream,  z-independent  equilibrium  through  the  use  of 
three  conservation  equations.  In  this  way,  the  chore  of 
computing  the  actual  beam  trajectories  all  the  way  from  the 
cathode  to  the  downstream  equilibrium  region  may  be  avoided. 

The  first  equation  is  the  conservation  of  canonical 
angular  momentum,  which  follows  from  the  assumption  of 
azimuthal  symmetry.  It  may  be  written  as 
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where  r is  the  initial  radius  of  an  electron  at  radius 
o 

r in  the  downstream  equilibrium,  and  B is  the  external 

zo 

guide  magnetic  field  in  the  cathode  region,  which  is  assumed 
to  be  known. 

From  the  assumption  that  the  equilibrium  is  steady- 
state  in  time,  the  continuity  equation  implies  conservation 
of  particles  within  annular  rings: 


n 

o 


8 


oz 


n 8 


dr2 
z dr2 


(14) 


while  the  force  equations  imply  an  energy  constant: 
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(15) 


where  (y  - l)mc2  is  the  final  energy  given  to  the 
o 

electrons  during  acceleration  within  the  diode  region. 

For  the  case  in  which  a foil  anode  is  used,  it  may  be 
considered  to  be  held  at  ground  potential  * 0,  so 
that  yq  may  be  interpreted  as  the  y-factor  of  electrons 
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at  the  instant  that  they  penetrate  the  anode  foil.  However, 
the  application  of  equation  (15)  is  not  restricted  solely 
to  the  case  of  foil  anode  configurations. 

Now  the  solution  for  the  equilibrium  may  be  obtained 
by  iteration.  Let  us  assume  that 

n = n + 6n  (16) 

c 

e =*  0 +60  (17) 

2 C 2 


where  n , 0 are  the  values  of  n,  0 at  the  center  of 
c c z 

the  beam,  and  the  ansatz  is  made  that  1 6 n | <<  n^, 

1 60  I <<  0 . Next  one  may  define 
z c 
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where  y = (1  - 82)- 1/2  and  it  is  assumed  that  e,  <<  1, 
c c 1 

e2  <<  1* 

From  equations  (6)  - (8)  the  electromagnetic  field 
profiles  are  determined 


B a $ e * - 2Tte  n 0 r 
<p  c r c c 


B ~ 
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B (l  * £,  4)  . 

zc  \ 1 a 2 / 


and  from  equation  (5)  the  equilibrium  precessional  velocity 
is  determined 


.1/2 

1 r 


1 

4 £ 2 


Y a I , r‘ 


r 


14 


l 


It  is  heuristically  worthwhile  to  evaluate  the 
necessary  equilibrium  conditions  (11)  and  (12)  at  this 
point.  To  lowest  order,  condition  (11)  becomes 


1 > 


(a) 


B 

z 


to 2 a 

= el/2 

2 Y fl  c 1 
c c 


I (kA) 

e 

5 B (kG)a(cm)y  6 
z z 


(25) 


and  condition  (12)  becomes 


2 w2  2 oj2  I (kA) 

pB  _ pc  _ e / 26 \ 

1 > Q2  y2  ft2  “ G2  0.7  B2 (kG)a2 (cm)Y6 
B c c z z 

where  y and  5 have  not  been  subscripted  in  the  final 

z 

expressions  since  we  have  not  yet  taken  account  of  their 

radial  variations.  It  may  later  be  seen  that  y82  0 r = a 

is  appropriate  for  condition  (25)  while  y3  § r = 0 is 

z 

appropriate  for  condition  (26). 

It  is  now  apparent  that  e!/2,  which  reflects  the 
% X 


first  necessary  condition  for  equilibrium,  also 


characterizes  the  strength  of  the  axial  diamagnetism 

(which  varies  as  e^)  and  of  the  beam-frame  processional 

velocity  [B  (r = a)  = e*^2).  Likewise  e_,  which  reflects 
<pB  1 2 

the  second  necessary  condition  for  equilibrium,  also  char- 
acterizes the  strength  of  centrifugal  effects.  Therefore, 
for  the  majority  of  cases  where  and  e2  are  quite 

small,  these  complicating  effects  of  diamagnetism,  rela- 
tivistic precession,  and  centrifugal  effects  may  be  treated 
as  small  corrections. 

If  the  diamagnetic  corrections  are  neglected  in  equation 
(24),  then  8.  “ r and  the  electron  beam  precesses  as  a 
rigid  body.  Given  the  precessional  velocity  from  equation 
( 24) , then  r2/r^  may  be  determined  from  equation  (13) . 

One  finds 


dr2 

dr2 
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zo 
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zc 


(27) 
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where  B zq  is  assumed  to  be  uniform  across  the  cathode 
surface.  To  lowest  order,  equation  (13)  implies  that 
the  electron  beam  propagates  on  magnetic  flux  surfaces; 
the  small  corrections  due  to  the  mechanical  angular  momen 
turn  acquired  by  the  electrons  are  the  ^ e terms  in 
equations  (27)  and  (28)  . 

From  equation  (14) , the  downstream  equilibrium 
particle  density  may  now  be  expressed  as 


I 

I 

I 


Finally,  the  radial  dependence  of  y may  be 
determined  from  solving  Poisson's  equation  in  the  down- 
stream equilibrium  region  where  the  z-variation  is 
negligible  compared  to  the  r-variation.  The  beam  extends 
from  r = 0 to  r = a,  with  a vacuum  gap  from  r = a 
to  the  conducting  waveguide  walls  at  r = b.  The  third 
conservation  equation  (15)  may  be  used  to  define  y in 
terms  of  the  electric  potential  <J>  for  all  r in  the 


range  0 < r < b.  Hence  the  boundary  condition  4>(r  = b)  = 


is  equivalent  to  y(r  = b) 


= y . One  thereby  arrives  at 
the  following  differential  equation  and  boundary  condition 
for  y ( r ) s 


Id  dy  4-rre2 

j = 2 n 

r dr  dr  me' 


(30) 


|y  + In 


r=a 


(31) 


where  the  Poisson  equation  has  already  been  integrated 
through  the  vacuum  gap  to  yield  the  modified  boundary 
condition  (31)  at  the  beam  edge.  Equations  (30)  and  (31) 
now  determine  the  amount  of  y-reduction  (i.e.,  slowing) 
which  the  electron  beam  undergoes  because  of  the  space 
charge  in  the  drift  tube. 

Equations  (29)  - (31)  may  be  recognized  as  specifying 
the  limiting  current  phenomena  which  we  previously  investi- 
gated in  the  limit  of  an  infinite  magnetic  guide  field 

. In  that  limit  however,  ■*■0, 


(i.e.,  see  Appendix  C) 


e -*■  0 , B -►  B , a -*■  a and  equation  (29)  is  simplified 
2 zo  zc  o 

considerably.  In  this  previous  limiting  current  calcula- 
tion, the  parameter  space  of  y and  b/a  was  explored 

o 

exhaustively  for  expressions  which  would  accurately  predict 

the  limiting  current  threshold.  However,  for  the  present 

calculation  we  will  assume  that  the  application  is  such  that 

operation  of  the  electron  beam  diode  near  the  limiting 

current  threshold  is  not  required.  Then  extreme  precision 

as  to  the  determination  of  the  limiting  current  should  not 

be  necessary,  and  the  first  order  iteration  which  we  have 

begun  should  be  of  sufficient  accuracy. 

Inserting  equation  (29)  into  equation  (30),  with 

6 - 6 , one  then  obtains  the  results 

z c 


where 
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and  we  shall  assume  that 


1 » 

yc 


(>  - 5 ■ J 


y b 

C C 
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2 ei(X  + 4 El) 


(35) 
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Here  condition  (35)  is  imposed  purely  for  theoretical 
simplicity,  since  a wide  range  of  interesting  practical 
applications  are  still  possible  which  obey  this  restriction 
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f * 





to  moderate  electron  currents.  Other  cases  could  also  be 
examined,  but  we  shall  not  do  so  herein.  If  one  is  within 
the  pencil  beam  regime,  for  which  y^3  <<1  + 2 In  b/a, 
then  condition  (35)  is  valid  up  to  the  limiting  current 
threshold — as  discussed  below.  Otherwise,  condition  (35) 
may  still  be  valid  if  the  beam  current  is  sufficiently 
below  the  limiting  current. 

If  equation  (33)  is  solved  for  v(y  ),  v is  found 

c 

to  have  a maximum  value  when 


= y. 


1/3 


'1L 


<To  - 


6(2  *>  r * 1) 


'cL 


(36) 


1/3 

where  e,  is  e,  evaluated  for  Y * Y - Y • The 
1L  1 c cL  o 

e correction  in  equation  (36)  reflects  the  y dependence 

XL  C 

of  in  equation  (33) . We  shall  assume  that 


2 / 3 


1 >> 


1L 
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2 in  - + 1 
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(37) 
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so  that  this  correction  is  small.  When  y = y , then 

c cL 

v = v = v , and  the  corresponding  maximum  value  of  I 
L max  e 

is  called  I , the  limiting  current.  It  is  the  largest 
L 

electron  current  which  may  be  propagated  within  the  drift 
tube,  and  is  given  by 


XL  “ 


(y2/3  - 1) 3/2  mc3/e 
o 

2 In  - + 1 - j e 

a 4 1L 


(38) 


Equations  (36)  and  (38)  are  expressions  of  first  order 

accuracy  within  the  pencil  beam  regime  for  the  limiting 

current  parameters  y „ and  I, , which  reflect  the 

cL  L 

finite  magnetic  field  influence  through  the  parameter 

e . In  the  infinite  magnetic  field  limit,  e 0 and 

XIj  1L 

these  expressions  reduce  to  the  usual  pencil  beam  inter- 
polation results.  It  may  be  noted  that  as  the  magnetic 

guide  field  is  reduced,  e,,  increases,  y , decreases, 

1L  cL 

and  I increases— in  qualitative  agreement  with  some 

L 

reported  numerical  simulations  of  this  problem. 
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(24) 


Using  the  expressions  for  0^  and  y in  equations 

<P 

and  (32),  one  may  solve  for  0^  from 

0 = (1  - y~2  - 02)l/2.  The  result  is 

z <p 


0 

Z 


0 

c 


1 + 


(39) 


which  exhibits  the  corrections  to  0 . From  condition  (35) 

c 

one  may  observe  that  the  e./2  y2  02  correction  is  smaller 

1 c c 

than  the  Ay/y 3 02  correction  by  a factor  of  t Hence 

c c 2 4 2 

the  dominant  radial  dependence  of  0^  derives  from  the 

space  charge  slowing,  which  affects  the  center  electrons 

relatively  more  than  those  on  the  edge  of  the  beam.  The 

ansatz  1 >>  1 50  1/0  » A y/y 3 0*  follows  from  the  assump- 

z c c c 

tion  Ie  <<  1^,  but  the  ansatz  may  also  be  obeyed  fairly 

well  even  when  I < I . 

e L 

If  equation  (39)  is  inserted  in  equation  (29),  an 
improved  expression  for  the  radial  dependence  of  the  beam 
density  is  obtained: 
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(40) 


Again,  our  previous  assumptions  imply  that  the  e /2  y2  6 2 

1 c c 

correction  is  smaller  than  the  Ay/Y 3 6 2 correction,  and 

c c 

that  the  ansatz  1 >>  |6n|/nc  is  obeyed.  However  the 

relative  size  of  the  diamagnetic  correction  and  the 

Ay/y 3 B2  space  charge  correction  is  not  determined  from 
c c 

the  preceeding  assumptions. 

If  the  results  from  equations  (39)  and  (40)  are 
inserted  into  equations  (6)  and  (7) , improved  expressions 


Equation  (23)  for  B is  still  correct  to  lowest  order  in 

z 

the  small  parameters. 

At  this  point,  one  could  continue  with  a second 
major  iterative  cycle  through  the  equations  to  obtain  still 
higher  order  corrections.  However,  since  we  expect  our 
various  assumptions  of  smallness  to  be  fairly  well  obeyed 
for  most  applications  of  interest,  we  will  stop  after  list- 
ing one  further  updated  expression  for  the  equilibrium 
precessional  velocity: 
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If  the  radial  dependences  herein  derived  are 
inserted  into  the  necessary  condition  (11)  for  the  existence 
of  equilibria,  there  results 
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Since  our  previous  assumptions  insure  that  [ (Ay/Y  ) + 

c 

(Ay/2  y3  S2)  “ (e,/4  Y2  62)]  < the  right-hand  side  of 

c c 1 c c 3 

equation  (44)  is  largest  at  r = a as  we  previously 

hypothesized.  The  small  correction  terms  to  the  zero  order 
1 /2 

result  1 > are  only  present  because  we  chose  to 

define  in  terms  of  on-axis  quantities  in  equation  (20) . 

Likewise,  condition  (12)  is  found  to  be  most  restric- 
tive at  the  center  of  the  beam,  and  it  reduces  simply  to 
1 > e2  since  was  originally  defined  in  terms  of  on-axis 

quantities. 
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APPENDIX  H 


LINEAR  THEORY  OF  COLLECTIVE  WAVES  SUPPORTED  BY 
AN  UNNEUTRALIZED  RELATIVISTIC 
ELECTRON  BEAM 


(Appendix  H Contains  34  Pages) 
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Among  the  important  applications  which  have  been 
suggested  in  recent  years  for  high  energy  relativistic 
electron  beans  are  the  use  of  such  beams  as  a medium  for 
the  collective  acceleration  of  ions.  An  important  subclass 
of  these  collective  acceleration  schemes  are  based  upon  the 
utilization  of  self-consistent  collective  eigenmodes,  or 
waves,  supported  by  the  electron  beam — within  which  waves, 
the  subject  ions  are  to  be  confined  during  the  process  of 
ion-wave  acceleration.  For  example,  the  Auto-Resonant 
acceleration  concept  seeks  to  utilize  the  negative  energy 
electron  cyclotron  wave,  while  the  Converging  Guide  acceler- 
ation concept  seeks  to  utilize  the  negative  energy  longi- 
tudinal plasma  wave.  We  will  examine  herein  the  derivation 
of  the  dispersion  relation  for  these  waves,  paying  attention 
to  certain  effects  of  centrifugal  forces,  Coriolis  forces, 
and  longitudinal  compression  which  are  negligible  when  the 

electrons  are  ultrarelativistic  (i.e.,  y”’2  <<  1),  but 

o 

which  must  be  retained  for  the  case  of  lower  energy  electrons. 
Likewise,  the  effects  due  to  an  axial  magnetic  guide  field 


I 

I 

I 

I 

I 


of  finite  strength  are  retained  in  general,  although  these 


effects  become  negligible  in  the  infinite  magnetic  field 


limit  which  is  often  adopted  when  plasma  waves  are  examined. 
The  more  general  analysis  which  is  presented  here  will 
permit  an  appraisal  of  the  validity  of  the  "ultrarelativistic 
and  the  "infinite  magnetic  field"  limits  which  are  frequently 
taken,  and  will  display  the  coupling  effects  between  the 
plasma  waves  and  the  electron  cyclotron  waves. 

The  mathematical  complications  which  arise  within 
linear  theory  may  be  traced  to  the  complications  inherent  in 
the  equilibrium  of  such  unneutralized  relativistic  electron 

beams.  The  self-fields  E and  B of  the  electron  beam 

r 4> 

combine  to  produce  a net  outward  force  - e(E  -3  B ) = - 

r z <t> 

e y“2  e^  acting  upon  the  electrons,  which  can  only  be 

counterbalanced  if  the  electrons  have  a precessional  velocity 

B^c  in  an  axial  magnetic  guide  field,  so  that  the  magnetic 

pinch  force  - e 0 B , reduced  somewhat  by  th«i  outward 

<p  z 

centrifugal  force  y me2  02/r,  restores  the  radial  force 

o 4> 

equilibrium.  This  precession  of  the  electrons  in  turn 
produces  a diamagnetic  reduction  in  the  axial  magnetic 
field  within  the  beam,  and  slightly  alters  the  radial  dis- 
tribution of  the  electrons  relative  to  that  at  the  cathode 

2 


before  the  precession  is  established.  Finally,  the  injection 
of  the  electron  beam  into  an  evacuated  waveguide  is  accom- 
panied by  some  r-dependent  slowing  of  the  beam  due  to  space 
charge  forces.  This  produces  a slight  outward  radial 
gradient  to  the  axial  beam  velocity,  and  corresponding 
r-dependence  for  the  beam  density  and  relativistic  y-factor. 
All  of  the  above  mentioned  effects  produce  complicated  radial 
dependence  in  the  beam  density  and  velocity  parameters,  which 

in  turn  cause  the  self-fields  E and  B and  the  diamag- 

r <J> 

netic  correction  to  the  axial  magnetic  field  B z to  have 
complicated  r-dependence. 

These  complicating  features  in  the  self-consistent 
equilibrium  quantities  produce  a corresponding  complication 
in  the  radial  eigenvalue  problem  which  must  be  solved  to 
yield  the  linear  wave  dispersion  relation.  The  extent  of 
the  complication  is  apparent  from  examining  the  linearized 
version  of  the  relativistic  force  equation: 


— (E  + 0 x B) 
me  — — — 


(1) 


3 


where  e = - |e|,  m are  the  electron  charge  and  mass, 


Y = (1  - 32)-1/2,  3 = 3 z + 8.  $ is  the  normalized, 

o-o—o  z <p 

zero-order  velocity  vector  of  the  electrons,  and  E = E 

-o  r 

B = B z + B d are  the  zero-order  electromagnetic  fields 
—o  z <p  3 

In  zero-order  (i.e.,  equilibrium),  the  radial  component  of 


the  force  equation  (1)  is 


3 B + 3 B ) 

z <p  <J>  z 


(2) 


while  in  first-order,  the  force  equation  may  be  written  as 


- — (£  + 3 x B + 8 x B)  (3) 

me  — — — o o — 

where  a tilde  denotes  first-order  (i.e.,  wave)  quantities. 
The  r,  %,  and  z components  of  equation  (3)  are  as 
follows : 
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Here,  the  downward  pointing  arrows  indicate  terms  which  are 
present  because  of  the  equilibrium  precession  of  the  beam, 
while  the  upward  pointing  arrows  indicate  terms  which  are 
present  because  of  the  self-fields  of  the  beam  equilibrium. 
There  are  seventeen  such  terms  in  these  equations,  all  of 
which  have  complicated  r-dependence  and  all  of  which  are 
missing  for  the  case  of  a charge-  and  current-neutralized 
beam,  for  an  ultrarelativistic  beam,  or  for  a beam  in  an 
infinite  magnetic  guide  field. 

It  is  clear  that  the  arrowed  terms  vanish  individ- 
ually for  a neutralized  beam.  The  ultrarelativistic  limit 

may  be  taken  formally  by  tending  y , B , n (hence  E , B ) 

o z r <p 

to  infinity  in  such  a way  that  the  ratio  of  any  two  such 
quantities  remains  finite.  In  this  way,  the  relativistic 
plasma  and  cyclotron  frequencies  remain  finite,  but  longi- 


tudinal mass  terms  and  precessional  terms,  which  vary  as 


tend  to  zero.  Thus  the  0.  terms  of  equations 


(4)  - (6)  tend  to  zero  and  the  remaining  self-field  terms 


cancel  to  with  0(y~2)  and  may  be  omitted.  One  also  finds 

o 


that  0 « y” 3 -*■  0,  so  that  there  is  no  longitudinal  beam 

z o 

compression  in  the  ultrarelativistic  limit. 


For  the  infinite  magnetic  field  limit,  (3,0, 

<p  r 

and  (3  tend  to  zero  as  B 1 , and  only  the  z-equation 

<p  z 

(6)  is  required  to  relate  £3  to  the  perturbed  wave  fields. 

z 

Hence  the  perturbed  electron  motion  is  strictly  axial  in 
this  case,  and  longitudinal  plasma  waves  are  allowed.  By 
comparison,  in  the  ultrarelativistic  limit,  the  perturbed 
electron  motion  is  purely  transverse,  and  electron  cyclotron 
waves  are  allowed. 

Although  neither  of  these  limits  is  appropriate  in 

general,  one  can  still  define  parameters  e^,  and 

1/2 

3(a)  - e /y  as  in  Appendix  G,  which  are  small  for  most 
<v  l o 

cases  of  interest  and  thus  simplify  the  theory.  (Notice 
that  each  of  these  parameters  would  vanish  within  the  ultra- 
relativistic limit,  or  within  the  infinite  magnetic  field 
limit.)  Of  the  above  three  small  parameters,  the  largest 
for  many  cases  of  interest  is  e reflecting  the  centrifugal 

and  Coriolis  force  effects.  The  parameters  e,  and  3 (a) 

1 4) 

represent  effects  of  axial  diamagnetism,  3 . * B forces, 

<P  ~ 

0.n  current  density  sources,  and  relativistic  beam-frame 
$ 

precession  which  have  a relative  magnitude  of  only  a few 
percent  or  less  for  many  cases  of  interest.  We  shall  neglect 
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these  latter  effects,  but  retain  the  large  centrifugal 

effects  within  the  linear  theory.  In  addition,  it  is 
necessary  to  retain  effects  such  as  longitudinal  compres- 
sion which  have  a relative  magnitude  of  y-2,  and  are  not 

o 

at  all  small  for  electron  beams  which  are  not  ultrarela- 
tivistic . 

Finally,  there  are  the  small  parameters  Ay/y  and 

c 

Ay/y3S2,  also  discussed  in  Appendix  G.  These  parameters 
c c 

reflect  the  radial  variation  in  beam  energy  due  to  effects 
of  space-charge  slowing.  They  are  small  only  if  the 
electron  beam  current  is  well  below  the  Alfven-Lawson  limits, 
which — fortunately — is  often  the  case.  Nevertheless,  these 
corrections  can  sometimes  be  as  large  as  ten  percent  or  more, 
in  which  case  they  should  be  included  in  the  analysis.  In 
addition,  these  parameters  do  not  formally  vanish  in  either 
the  ultrarelativistic  or  the  infinite  magnetic  field  limits, 
and  they  consequently  have  relevance  for  many  accelerator 
applications.  (In  fact,  for  the  converging  guide  accelerator 
these  parameters  may  have  to  be  appreciable  in  order  for  the 
plasma  wave  phase  velocity  to  be  somewhat  low  at  the  front  of 
the  accelerator.)  However,  because  these  effects  produce 
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radial  dependence  which  greatly  complicates  the  mathematical 


eigenvalue  problem  within  linear  theory,  we  shall  neglect 
them  for  the  present  calculation.  This  will  allow  the  linear 
theory  to  be  solved  in  terms  of  relatively  simple  eigenfunc- 
tions (i.e.,  Bessel  functions),  and  will  allow  us  to  examine 
the  modifications  which  arise  due  to  the  retention  of  the 

e_  and  y-2  terms  which  are  omitted  for  ultrarelativistic 
2 o 

beams.  Also,  the  coupling  between  the  plasma  and  cyclotron 

waves  may  be  examined  within  such  a calculation. 

Elsewhere,  the  effects  of  the  radial  variation  in 

Y on  the  electron  cyclotron  wave  are  considered  within  the 

ultrarelativistic  limit  in  Appendix  D,  and  the  effects  on 

the  longitudinal  plasma  wave  are  considered  within  the 

infinite  magnetic  field  limit  in  Appendix  I. 

The  approximations  discussed  above  for  this  linear 

theory  calculation  imply  that  the  electron  density  n,  the 

axial  velocity  8 c,  and  the  kinetic  energy  (y-l)mc2  may 

z 

be  treated  as  constants  with  respect  to  r,  as  may  be  seen 

from  the  expressions  derived  in  Appendix  G.  Hence  we  shall 

denote  them  as  n = n , 8 = 8 , y * Y . In  addition  we 

o z o o 

may  let  y = (1  - 8Z)”I/2  since  the  8.  corrections 
1 o o <t> 

9 


....  y-  .. 


---  '.'.aw-  TT 


to  Yo  will  be  negligible  under  these  approximations. 
Furthermore , the  axial  magnetic  field  B may  be  treated 
as  uniform  in  r,  and  we  may  define  the  usual  character- 
istic frequencies 


Y me 
o 


4 tt  n e2 
o 

Y m 
o 


and  the  small  parameters 


1 \ 2 Y ftc  / 

' o ' 


<<  1 


/ 2 0)2  \ 

•M 


<<  1 


The  equilibrium  self-fields  of  the  beam  are  then  given  by 


10 


(11) 


B = 8 E » - 2 Tie  n 8 r 
0 o r o o 


and  the  equilibrium  velocity  of  precession/c  is 


8 


♦ 


(12) 


Among  the  approximations  invoked  for  this  calculation  are 
the  orderings 

el  <<:  e2  <<  1 (13) 

and 


e,  <<  0A(a)  <<  y~2,  1 
i <p  o 


(14a) 


No  assumption  is  invoked  about  the  size  of  y-2,  so  that 

o 

the  calculation  may  be  relevant  even  for  a nonrelativistic 
beam.  Assumption  (14a)  may  also  be  written  as 


1 >:>  e!/2  * Y2  6.  (a)  ■ 

o 1 o <|> 


Er(a) 

B 

z 


(14b) 
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There  will  exist  some  cases  of  interest  for  which  condition 


(13)  is  fairly  well  obeyed,  but  condition  (14)  is  not. 

However,  for  these  cases,  y”2  <<  1 and  the  violation  of 

o 

condition  (14)  merely  means  that  the  y”2  effects  are  no 

o 

larger  than  the  neglected  8 and  e,  effects,  and  so  the 

<P  1 . 

— 2 

yQ  effects  should  be  neglected  as  well. 

Returning  again  to  the  linear  theory  of  waves 
supported  by  such  an  equilibrium,  the  force  equations, 
continuity  equation,  and  Maxwell's  field  equations  are 
linearized  and  transformed  with  respect  to  the  basis  func- 
tions exp  i(!<t>  + k z - oot)  . It  is  convenient  to  define 

z 

the  following  wave  field  variables 


A+  = E + i E + i 8 B - 8 B (15) 

r <p  or  o 0 


A~  i E - i E - i 8 B - 8 B . (16) 

r <}>  or  o <j> 


For  this  case  in  which  8 is  constant,  Maxwell's  equations 

o 

+ 

imply  the  following  pair  of  equations  for  A" 


12 


1 J_  J_ 
r 3r  r 3r 


U t l)2 


- k2  + 


OJ 


■)a!  - 


4-rriAw  . z . 

-rr~  (J  ± i JJ 


4’  (£  * r)(5  * f 3z) 


(17) 


where  Aw  = w - k 0c  and  p , J are  the  bean  charge  and 

z o — 

current  density.  In  general,  there  will  be  boundary  condi- 
tions at  the  conducting  wall  at  r * b to  ensure  that  , 

E , B vanish  there,  and  in  addition  there  must  be  jump 

<P  r 

conditions  at  the  beam  edge  r = a to  properly  patch  the 
beam  interior  field  solutions  to  the  vacuum  gap  field  solu- 
tions. However,  a great  deal  of  this  complexity  may  be 
avoided  when  the  vacuum  gap  is  narrow  (i.e.,  b - a <<  a), 
since  then  only  the  beam  interior  solutions  are  required 
and  the  boundary  conditions  may  be  applied  at  the  beam  edge. 
We  shall  make  this  assumption,  in  which  case  the  boundary 
conditions  reduce  to 
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0 

J 

I 


r-. 


1 


B = 
z 


*(Vo-;)lr 


:r  r (A  - A ) + - (A  + A 
r d r r 


B = 


2 (k  0 - -) 

\ z c c / 


+ - 21  - 
ik  (A  - A ) + — E 
z r z 


B<0  = 


(k  B - ~) 

\ 2 O C / 


2i  77  + k2  (A+  - A"> 


The  perturbed  beam  charge  and  current  densities  are  related 
to  the  perturbed  velocities  and  number  density  as  follows 


p = - e n 


J = - ec  (B  n + n 0) 
z o o z 


J = - ec  n 6 
r o 


i 


* 


where  approximation  (14)  is  invoked  to  simplify  equation 

46  o 

(32)  , and  iu  h u - k S c - = Aw  - £cj  . Here 

z o r <p 

= 3^c/r  = — S2  is  the  precessional  frequency  of  the 
beam  equilibrium. 

Finally,  the  continuity  equation  is  employed  to 
relate  n to  6 as  follows 


n 


n 

o 


k c 

~~  0 

6w  z 


+ 


C 1 J_  ~ 
i6oj  r 9r  r r 


+ 


Ic 

6wr 


(33) 


It  may  be  observed  that  the  z^  effects  are 
represented  through  the  ^-portion  of  the  convective  deriva- 
tive terms  (i.e.,  j Iz  ^l)  and  through  the  centrigugal 

and  Coriolis  force  terms  in  equations  (30)  and  (31),  which 

2 C6d,  1 

are  the  terms  with  coefficients  = — e„.  These  latter 

r n 2 2 

effects  may  be  conveniently  retained  by  introducing  the 
modified  variables 


(34) 
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...  y 


T 


in  terms  of  which  equations  (30)  and  (31)  become 


i ^ 8 - 8 = 

5 r * 2 i 


(35) 


I 

1 

f 


5 + iis  5 = ft— -ft 
r a 41  2i  B 


(36) 


The  source  terms  + i in  equation  (17)  require  only 


A /V 


&r,  3^  which  may  be  readily  related  to  A via  the  r, 
force  equations  (35)  and  (36).  However,  the  source  terms 


3q 

P - — J require  y-2  n - n 6 3 , which  couples  the  con- 

c z o o o z 


tinuity  equation  and  the  z force  equation  to  Maxwell's 
equations. 

From  equations  (28),  (29),  (35),  and  (36),  one  finds 


0 

i 


4iri  Aw  i . 

(Jr  1 1 V 


u)2 

= _£ 
c2 


Aw 


6oj  ± 
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(37) 


I.  ■ - 


y 


while  from  equations  (20),  (26),  (27),  (32 
and  (36) , one  finds  that 


/ 0 

V 

( 

1 o 

C-t« 

II 

• 2 T ■ 

M 

1 3 

r 3r  r 

+ 

Z 

+ — 

A 

a" 

r 

1 + 5 „ 

1-6 

L z 

Z . 

B iw  1 

4-  — ■■  ■■ 

O 0 

1 3 

+ — 

k c 

■ - 0 CO 

r 3r  r 

(A  + A ) 

z 

o 

_A 
1 + 


where 


One  nay  now  insert  equations  (37)  and  (38)  for  the  self- 


consistent  beam  source  terms  into  the  pair  of  equations 

(17)  to  obtain  a pair  of  coupled,  second-order  differential 

+ 

equations  for  the  wave  field  variables  A . First  however, 

it  may  be  noted  that  the  complicating  effects  associated 

with  the  #,u)  terms,  which  are  present  also  in  6w,  6^, 

and  , will  be  removed  for  the  important  cases  of  1 = 0 

(i.e.,  ^-symmetric  sausage  modes  or  compressional  modes)  or 

w,  * — -*■  0 (e.g.,  the  infinite  magnetic  field  limit). 

9 4 2 

When  e_  <<  1,  these  wA  effects  will  tend  to  influence 
2 

the  plasma  and  cyclotron  waves  significantly  only  at  large 
values  of  l — where  some  fluting  perturbations  might  be 
anticipated  for  the  electron  beam. 

Let  us  briefly  consider  the  infinite  magnetic  field 

limit,  for  which  w -*-0,  6w  -*■  Aw,  6 „ -*•  0,  J •+•  0,  J -*•(). 

<p  x.  r 0 

In  this  limit,  the  pair  of  equations  (17)  may  be  expressed 
as 


20 


r 


[~  ~ r ~ It  t II  ! ,1 

lr  3r  3r  ^ -*z+Jt|a1 


(3r  ^ r)fr  57  r (A++A-,  + i (a+.a-( 

. 


where 


T s V%  - 0) 


2 Y~  A“2 
O 


y*  Au)2 

O 


The  boundary  conditions 


U8,  and  U9)  then 


(A+  - A")  = 0 . [l  3 + 1 

r*a  [r  i?  r ^ + A-J 

Jr 


Pair  of 


equations  (41)  may  be 


Put  in  the  form 


M1+2T)  i-i  3 £ 2 2 1 

( r i?*  - ;r- #♦  *r  ,*♦+*-,  .aiL  rA  + iJSj  . 

J r [ 9r  rJA~A) 


(44a) 


— r- — 


II 


I 

I 

1 

I 

f 
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from  which  it  is  apparent  that  the  modes  decouple  for 
1 = 0,  since  then  there  is  one  equation  and  one  outer 
boundary  condition  (in  addition  to  regularity  at  r = 0) 
for  each  of  the  variables  A+  + A and  A+  - A . Equation 
(44a)  is  seen  to  yield  the  longitudinal  plasma  wave  disper- 
sion in  this  case,  while  equation  (44b)  describes  electro- 
magnetic waves.  To  discover  the  decoupling  which  might 
occur  for  1/0,  one  is  led  to  define  the  variable 


13  4*  — 4.  __ 

R = - — r (A  + A ) + - (A  - A) 
r 3r  r 


since  this  combination  appears  on  the  right  hand  side  of 


equation  (41)  and  one  also  finds  that  E * r.  Sure  enough 


the  equations  (44a)  and  (44b)  imply  that 


while  the  boundary  conditions  (43)  imply  that  R(r  = a)  =0. 


It  follows  that  R(r)  = J^(k^r),  where  the  dispersion 
relation  is  given  by 


H-r)  - A) 


(47) 


and  the  are  the  zeroes  of  J^(A)  = 0 (e.g.,  XQ1  = 2.4) 

Equation  (47)  is  the  desired  dispersion  relation  for  longi- 
tudinal plasma  waves  in  the  infinite  magnetic  field  limit. 

One  may  also  define  a variable  conjugate  to  R which  may 
be  shown  from  equations  (44)  to  obey  an  equation  which 
describes  electromagnetic  waves.  However,  we  shall  not 
explore  these  waves  here. 

The  second  case  which  we  wish  to  examine  in  some 
detail  is  that  of  l =*  0,  ^-symmetric  sausage  modes  or 
compressional  modes  of  oscillation.  This  case  is  important 
for  the  applications  of  collective  ion-wave  acceleration 
since  it  is  the  l * 0 waves  within  which  the  ions  are  to 
be  trapped  and  accelerated.  For  £ = 0,  it  is  again  true 
that  6d)  -►  Aw,  and  we  shall  employ  the  variables  T and 
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il 
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X 


(48) 

u) 2 

A-y._  (49) 

C Aw  ± ft 

If  the  sources  of  equations  (37)  and  (38)  are  then  inserted 
into  equations  (17)  - (19),  there  results 

l 
I 

r 


/ a i a ,\  ± a l a / a*  . a~  \ 

(£  7 3?  r ■ p±)  A “ ■ T r s7  c U + « + i-s I 


(50) 


A ) 


r=a 


0 


(51) 


r 
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The  pair  of  equations  (50)  would  be  uncoupled  in  the 
ultrarelativistic  limit  where  T vanishes,  but  otherwise 


I 

t 


they  are  coupled  due  to  the  y~2  compressional  effects. 
This  coupled  set  of  second  order  equations  is  equivalent 
to  the  following  uncoupled  pair  of  fourth  order  equations 

J_  l ± _ P+tl  + «Hl-*+T) 

3r  r 3r  r ~ 1 - 6 2 + 2T 


_3_  1 _3_ 
3r  r 3r 


i 

I 

I 

! 

1 


p2  (1-<S)  (1+6+T)  + T2p2p2(l-62) 

1 - 6 2 + 2T  A = (1  - <52  + 2T)  2 A 

(53) 

which  nay  be  solved  in  terms  of  Bessel  functions  since  the 
coefficients  are  independent  of  r.  If  one  makes  the 
ansatz 


then  equation  (53)  reduces  to  a quadratic  equation  for 

k2  which  must  be  obeyed  by  both  k2  and  k2  , which  are 
1 xl  ±2 

therefore  the  two  roots  of  the  quadratic: 


k2 

i 


+ 


p2  (1  + 6)  (1  - 6 + T) 
1 - 6 2 + 2T 


k2 

l 


p2  (1  - 6)  (1  + 6 + T) 
1 - 6 2 + 2T 


t2  p+  p!(1  - 62> 

(1  - 6 2 + 2T) 2 


(56) 


If  equations  (54)  and  (55)  are  inserted  into  equations  (50) 
and  the  boundary  conditions  (51)  and  (52),  one  finally 
arrives  at  the  dispersion  relation,  which  may  be  cast  in 
the  form 


kna  Jo(kna) 
Ji(kna) 


+ G. 


kl2a  Jo(ki2a) 

Vkl2a) 


(57) 
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where 


G,  = 


(58) 


and 


S2  ' ' Gl(kh  * kI2> 


(59) 


The  familiar  infinite  magnetic  field  and  ultrarelativistic 
limits  correspond  in  this  notation  to  5 -*■  0 and  T -*■  0, 
respectively.  For  an  improved  analysis  which  retains  the 
effects  of  coupling  between  the  cyclotron  and  plasma  waves 
and  yet  exploits  the  smallness  of  the  e parameter  which 
characterizes  the  frequency  separation  of  these  two  modes, 
one  merely  treats  either  6 or  T as  a small  (but  nonzero) 
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parameter — depending  upon  which  mode  is  being  examined. 


We  conclude  this  calculation  by  considering  the  case  of 
the  negative  energy  cyclotron  mode,  for  which  we  may 
anticipate  that  Aw  - - ft,  hence  6 = - 1 and 
T = j <<  1.  Now  let  k*  be  that  solution  of  equation 
(56)  which  is  near  the  zero  of  the  first  square  bracket 


P*(l  + 6) (1  - 6 + T) 
1 - 6 2 + 2T 


and  likewise  k2  will  be  the  remaining  solution  which  is 
near  the  zero  of  the  second  square  bracket 


p2  (1  - 5)  (1  + 5 + T) 
1 - 6 2 + 2T 


Referring  to  equation  (49),  one  notes  that  p^  < 0 is 
associated  with  the  negative  energy  cyclotron  mode,  and 
therefore  k2^  > 0 and  k^  will  be  the  corresponding 
perpendicular  wave  number  for  that  mode.  Likewise,  when 


- Y“  - 


I 


A< i)  = - ft,  p2  >0  and  thus  k2^  < 0.  Consequently,  the 


Bessel  functions  in  equation  (57)  of  argument  k a 

j.2 


actually  behave  like  the  corresponding  modified  Bessel 


functions  Iq  and  1^  of  real  arguments.  Furthermore, 


for  |k^al2  = P2a2  ^ these  Bessel  functions  may  be 


expanded  as 


k „a  J (k  a ) 
12  o i2 


J!(kl2a) 


- 2 — ~ k2_a2  » 2 + 


4 12 


1 2 2 

4 P-a 


(62) 


The  functions  and  may  also  be  expanded  for 


6 = - 1,  |p2/p^|  <<  1,  T <<  1 to  yield 


GL  “ 1 


(63) 


°2  * 1 + T 


1 ♦ £ (i * rh) 


so  that  the  dispersion  relation  (57)  becomes 
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kxia  vkna) 

Vkxia) 


(64) 


Since  the  right  hand  side  is  approximately  constant, 
equation  (64)  approximately  determines  the  numerical  value 
of  kila'  and  then  equation  (56)  becomes  the  dispersion 
relation  relating  w and  k^.  The  smallest  solution  of 
equation  (64)  is  given  by 


p2a2 


\ia 


= 2.95  + 


5.66 


+ 2T 


1 + 


) 


» 2.95 


(65) 


Finally,  if  equations  (34),  (42),  (48),  and  (49)  are 

inserted  into  equation  (56),  the  dispersion  relation  for 
linear  waves  becomes 
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Aoj 2 (k2c2 


/ 2u> 2 \ 

(ft2  - — ^ J (k2c2  - oj2) 


GJ 

1 (k2c2  - oj2) 

Y z 
o 


- OJ2) 


0J‘  \ 

Aw2  - ) (k2c2  + oj2  - oj2)2  (66) 

Y o'  P 

where  k2  = k2  + k2  and  Aw  = oj  - k B c.  This  disper- 
il  z z o 

sion  relation  contains  eight  waves:  four  light  waves,  two 

longitudinal  plasma  waves,  and  two  cyclotron  waves.  If 

equation  (66)  is  examined  for  the  case  of  oj2  <<  k2c2  < k2c2 

z 

for  which  the  four  light  waves  are  decoupled,  it  is  quadradic 
in  Aoj 2 , and  the  solution  which  represents  the  two  cyclotron 
waves  may  be  written  as 


t 
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one  obtains 


ft2  - 


(k2c2  + oj‘ 


OJ 

+ -4 


Finally,  for  oj2  <<  k c , y2Qz  one  obtains 
P o 


Au  = u - k B c ^ 
z o 


- n i - 


as  the  dispersion  relation  for  the  negative  energy  cyclotron 

mode.  Here  the  correction  term  - oj2/y2n2  is  due  to  the 

P o 

centrifugal  and  Coriolis  forces  on  the  electrons,  while  the 


correction  + o)2k2  /2k2y2ft2  is  due  to  the  longitudinal  com- 
p il  o 

pression  of  the  electrons  which  couples  longitudinal  plasma 


modes  to  the  cyclotron  modes  through  the  T terms. 


The  dispersion  relation  (70)  is  the  starting  point 
for  all  subsequent  linear  analysis  which  must  be  performed 
to  examine  the  nature  of  the  negative  energy  cyclotron  waves. 

Having  developed  a framework  within  which  the  linear 
mode  analysis  of  cyclotron  and  plasma  waves  may  proceed 
inclusive  of  coupling  effects,  we  will  conclude  by  reiterat- 
ing that  it  is  the  smallness  of  the  and  parameters 

of  equations  (9)  and  (10)  which  validates  the  frequently 
adopted  limits  of  "infinite  magnetic  fields"  (wherein 
6 = 0 (e.  ^ ) <<  1 for  plasma  waves]  or  of  "ultrarelativistic 

beams"  (wherein  T = 0 (e2)  <<  1 for  cyclotron  waves].  In 

addition,  it  is  the  orderings  of  conditions  (13)  and  (14) 

and  the  neglect  of  Ay/y  effects  which  allow  solution  for 

c 

the  beam  eigenfunctions  in  terms  of  Bessel  functions. 
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APPENDIX  I 


LINEAR  THEORY  OF  PLASMA  WAVES  SUPPORTED  BY  A 
RADIALLY  VARYING  ELECTRON  BEAM  IN  A 
MAGNETIZED  WAVEGUIDE 


(ADDendix  I (’nnt-.ai  nq  57  Psapql 
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The  possibility  of  utilizing  the  negative  energy, 
longitudinal  plasma  wave  supported  by  an  electron  beam  for 
collective  ion  acceleration  (i.e.,  the  converging  guide 
accelerator)  has  been  previously  studied  in  some  detail, 
but  mainly  in  the  limit  of  an  infinite  magnetic  guide  field 
and  in  the  limit  of  pencil  beam  geometry,  for  which 

2/3 

Y <<1  + 2 tn  b/a  [where  (y  -l)mc2  is  the  injection 
o o 

energy  of  the  beam  electrons  and  b/a  is  the  ratio  of  the 
radius  of  the  conducting  waveguide  to  that  of  the  beam] . 

One  of  the  difficulties  encountered  in  these  prior  calcula- 
tions [e.g.,  R.  J.  Briggs,  Phys.  Fluids  1_9,  1257  (1976)] 
is  that  the  phase  velocity  of  the  plasma  waves  can  only  be 
made  small  (for  loading  the  ions  into  the  wave)  by  increas- 
ing the  electron  current  extremely  close  to  the  limiting 
current  which  can  be  propagated  without  virtual  cathode 
formation.  It  is  of  interest  to  explore  whether  these 
unfortunate  results  still  hold  for  more  realistic  nonpencil 
beam  geometries  where  the  radial  variation  in  the  equilibrium 

quantities  y,  6 , n is  significant. 

z 

For  this  analysis,  we  shall  adopt  the  infinite 
magnetic  field  limit.  Then  the  beam  equilibrium  calculation 


1 


will  be  just  as  described  in  the  limiting  current  publication 


of  Appendix  C,  to  which  the  equilibrium  analysis  of  Appendix 

G reduces  in  the  limit  e -►  0,  e 0.  To  this  equilibrium 

analysis  we  shall  couple  a linear  theory  calculation  of  the 

plasma  waves,  which  differs  from  that  in  Appendix  H in  that 

we  shall  here  focus  on  the  Ay/y  effects  and  ignore  the 

c 

finite  magnetic  field  effects. 

We  begin  by  reviewing  the  equations  for  the  equi- 
librium quantities.  For  a steady  state  electron  beam  diode, 
two  constants  of  the  motion  for  electrons  are 


(1) 


n 0 = n (r)  0 (2) 

z A A 

where  - e,  m are  the  charge  and  mass  of  the  electrons, 

(y  - l)mc2  is  the  electron  energy  at  the  anode,  0c  is 
A A 

the  axial  electron  velocity  at  the  anode,  n^  is  the 
electron  number  density  at  the  anode,  and  < 0 is  the 

electrostatic  potential  at  the  cathode.  It  is  assumed  that 


the  anode  is  grounded  electrically  to  the  waveguide  walls 

(i.e.,  1>  = 0)  and  that  the  electrons  are  born  at  the 

A 

cathode  with  zero  energy  (i.e.,  y =1).  since  $ is  an 

c c 

equipotential , y must  be  independent  of  r,  and  to  con- 

A 

form  to  prior  notation  we  will  let  y * y , 3 = B * 

A o A o 

(1  - Y~  ) • In  principle,  n^  may  vary  with  r if  a 

way  is  devised  for  the  cathode  to  emit  nonuniformly . There- 
fore we  shall  let  n = n g(r),  where  n is  the  mean 

A o o 

density  of  the  beam,  such  that 


I! 

I 

I 

I 

I 

I 

| 

I 


V = 


TT 


8 


o 


I 

e 

mc3/e 


(3) 


Here  a is  the  beam  radius  (constant  for  a uniform,  infinite 

magnetic  field),  I is  the  electron  current,  and  mc3/e  = 

e 

17  kA.  Consequently,  g(r)  must  have  the  normalization 


a 

1 * ^ | dr  r g (r)  (4) 

o 

Finally,  the  radial  profile  for  the  equilibrium  quantity 
y(r)  is  determined  from  the  radial  Poisson  equation 

3 


I 

i 

I 


T 


r < a (5a) 

a < r i b (5b) 

subject  to  y(r  -*■  0)  being  regular,  and 

Y (r  = b)  = y , (6) 

o 

and  where  S * (1  - y~2) 1/2.  Equations  (5b)  and  (6) 
z 

may  be  combined  to  yield  the  alternative  boundary  condition 


1 d 


dy 


r dr  dr 


4vg  (r) 
a2  0 


(y  + In  - r — ) = y (7) 

\ a dr / o 

' 'r=a 

Techniques  for  the  solution  of  these  equations  are  dis- 
cussed in  Appendix  C.  In  principle,  equation  (5a)  determines 

Y (y  > v,  r/a) , where  y = Y(r  = 0).  Then  the  boundary  con- 
c c 

dition  (7)  determines  y (y  / v,  b/a) , or  equivalently 

o c 

v (Y  , y , b/a).  As  a function  of  y i \>  has  a peak  value 

CO  C 


x 


■ 


v (y  > b/a)  when  y = y , (Y  , b/a) . Several  numerical 
L o C cL  O 


solutions  for  v (y  » Y » b/a)  are  illustrated  in  Figure  1, 

c o 

which  may  be  compared  with  Figure  1 of  Appendix  C. 

The  condition  of  limiting  current  is  3v/3y  = 0, 

c 

which  may  be  seen  to  correspond  to  the  condition 

9y  /3y  = 0,  since  0 = 3y  /3y  + Oy  /3v)(3v/3y  )•  From 
o c o c o c 

equation  (7),  this  limiting  current  condition  becomes 


_1  _d_  dT  4vg  (r ) 

r dr  r dr  Y3S3 

z 


(9) 


subject  to  the  interior  boundary  condi ton 


Appendix  C describes  how  a scaled  version  of  equations  (5a), 

(7),  (8),  (9),  and  (10),  using  the  dimensionless  variable 

y = vr2/a2,  may  be  used  to  achieve  an  efficient  numerical 

solution  for  v , y from  any  initial  choice  of  y , b/a. 

L o cL 

Consider  now  the  equations  for  linearized  waves, 

transformed  with  respect  to  the  basis  functions 

exp  i(£9  + k z - wt) . In  the  limit  of  an  infinite  magnetic 
z 

guide  field,  it  follows  that  8.  = 0 for  the  beam  equi- 

0 

librium  and  8 =0=8  for  the  linearized  waves.  [Here 

r 9 

a tilde  denotes  first  order  (i.e.,  wave)  quantities]. 
Therefore,  the  vector  potential  components  A^  = Aq  = 0, 
and  the  Lorentz  guage  condition  implies  that  A = — — 

Z K C 

z 

As  a result,  the  various  wave  field  components  may  be 
obtained  from  the  perturbed  potential  $ as  follows: 


~ d<J> 

E = - — 
r dr 


(11) 


~ 0)  ~ (0  ~ 

B = - E„  = ; — $ 

r k c 9 k c r 
z z 


~ _ co  ~ _ (o  d<t> 

9 k c r k c dr 
z z 


B * 0 
z 


The  perturbed  axial  equation  of  motion  is 


E 

~ e z 

i A<o  8 = — — r 

z me  y 


where  Aw  = w - k c 8 (r).  The  perturbed  continuity  equation 

z z 


Aw  n - k c n8  = 0 
z z 


Using  equations  (13),  (17),  and  (18)  to  relate  the  perturbed 


charge  density  p = - e n to  <t>,  the  perturbed  Poisson 


8 


equation  may  be  cast  in  the  form 


\k  ± r — - — + Ik2  - — ) \-  1 + 4v^(r)°2  I!  J = 

(rdr  dr  r2  \ z c2  / Y3S  Aw2a2j 


It  is  understood  that  g(r)  =0  by  definition  in  the  vacuum 
gap  region  a < r £ b.  Equation  (19)  for  <t>  is  subject  to 
the  boundary  conditions  that  $(r  -*■  0)  be  regular,  and  that 


$ (r  3 b)  *0 


so  that  E . E , and  B vanish  at  the  conducting  wall. 

0 z r 

In  the  same  fashion  that  the  alternative  boundary  condition 
(7)  was  derived  for  the  equilibrium  differential  equation, 
we  may  combine  an  analytical  solution  to  equation  (19)  for 
<J>  in  the  vacuum  gap  with  equation  (20)  to  obtain  the 
alternative  boundary  condition 


'V  • 


where  I , are  the  usual  modified  Bessel  functions  of 

the  first  and  second  kind,  and  I ' , K'  are  the  derivatives 

X / K> 

of  these  functions  with  respect  to  their  arguments. 

It  may  also  be  shown  that  the  solution  to  equation 
(19)  for  4>  which  is  regular  near  r = 0 must  vary  as 

~ l 

4>  « r . We  will  henceforth  restrict  attention  to  the  case 
of  l = 0,  9-symmetric  plasma  waves,  since  these  eigenmodes 
are  the  ones  that  would  be  utilized  for  the  collective 


10 


acceleration  applications.  For  1 = 0,  without  loss  of 
generality  we  may  assert  the  interior  boundary  condition 


4>(r  = 0)  = 1 


(221 


since  only  the  shape  of  >t>(r)  is  determined  from  equations 
(19)  and  (21). 

The  r-dependence  of  the  coefficients  of  equation 
(19)  is  generally  quite  complex,  since  the  y and  6 

z 

r-dependence  generated  by  equations  (5)  - (7)  appears  not 

only  through  l/y36  , but  also  through  Am.  In  addition, 

z 

g may  have  r-dependence.  However,  we  shall  first  consider 
the  special  case  in  which  all  of  this  r-dependence  is  weak, 
since  analytical  expressions  may  then  be  developed  for  the 
eigenfunction  $ and  for  the  plasma  wave  dispersion  rela- 
tion. Let  us  also  consider  the  case  for  which  b = a and 
Yq  >>  1,  within  the  "fat  beam  limit,"  which  is  opposite  to 
the  "pencil  beam  limit"  within  which  the  previous  analytical 
expressions  were  derived  (e.g.,  by  R.  J.  Briggs,  ibid. ) . Let 
g(r)  = 1 and  v <<  so  that  the  r-dependence  of  the 

11  I 

I 


f 


equilibrium  is  weak.  From  equations  (5a)  and  (7),  one  finds 


Y (r)  = y 


.-fA'-i) 


(23a) 


where 


S (r)  * 8, 
z k 


(23b) 


to  lowest  order,  and  the  value  8^  about  which  we  expand 
may  be  selected  later  for  convenience.  (For  greatest 
accuracy,  8,  should  be  the  value  of  8 at  some  inter- 

K Z 

mediate  value  of  r between  0 and  a.) 

Equation  (19)  may  be  put  in  the  dimensionless  form 


ia's*  kl(R)a!  ;<R)  - 0 


where 


k2(R)a2  = (1  - 62)  - k2a2  + — r 

I A V 3 


Y 8 (8  - 8 ) 

Z Z (p 


R = 


r/a 


<p  = w/k  c 
z 


and  boundary 


conditions  (21)  ^ 


f22)  become 


*(R  = 


0)  = z 


(26a) 


*(R  = l)  = 


Let  k*(R)  = , 2 

i0  *‘R)  and  ;(R)  . .(o) 
rhen  the  zero-o-..-.  ' * «)  + i«> 


(26b) 


order 


ei9enfunction 


obeys 


(R)  . 


lR  dR  R dR  + 


0 (°) 


(R)  , 


(27a) 


*(°\r  = 


0)  = ! 


(27b) 


0 (o) 


(R  * 1)  a 
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(27c) 


i 


I 


whose  solution  is 


~ (o) 

$ (R)  = J (k  a R)  = J (k  r) 

o iO  o 10 


(28) 


where 


k 


10 


X 

on 


2.4, 


(29) 


The  first  order  eigenfunction  must  then  obey 


— — R — + k2  a2 
R dR  dR  Ki0 


5 (1)  (R) 


<Sk2  (R)  a2  i (o)  (R) 


(30a) 


$ (1)  (R  = 0)  = 0 


(30b) 


~ (1) 

* (R  = 1)  = 0 


(30c) 


The  Rayleigh-Ritz  procedure  may  be  used  to  generate  the  dis- 
persion relation  for  plasma  waves,  correct  to  first  order, 
from  equations  (30) : 


14 


Therefore,  the  weak  r-dependence  in  the  coefficient  of 
equation  (24)  should  be  averaged  over  the  zero-order  eigen- 
functions of  the  equation,  to  determine  the  dispersion 

relation.  Let  8,  now  be  chosen  such  that 

k 

<V3>  = Y~3  = (1  - Bj)  3/2  . (34) 


The  dispersion  relation  (32)  may  be  seen  to  agree  with  the 
result  obtained  in  Appendix  H,  equation  (47)  if  one  identi- 


fies y , 6 with  the  constant  values  y , 8 of  Appendix 

J\  K O O 

H.  It  is  a quartic  equation  for  8,  as  a function  of  k . 

<p  z 

If  one  expands  for  8.  = 8.  + 68.,  |68j  <<  8.  , (28,  y 2 ) " 1 , 

<p  k (j)  (J)  k k k 

one  obtains 


CO 

8 = = 6 - 
♦ k z Pk 
c 


2 v 


1/2 


.5/2 


+ Y 


— 2 1 


(35) 

where  y,  - y - 0 (v) , 8,  = 8 - 0 (v) , as  determined  from 

k o k o 
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equations  (23) , (33) , and  (34)  . The  following  features  may 

be  observed  about  the  scaling  of  the  plasma  wave  phase 

velocity  in  equation  (35)  : 

(i)  the  phase  velocity  increases  weakly  as  a 

function  of  k^,  although  this  increase 

will  not  be  noticeable  until  k is  larger 

z 

than  y k ; 
k 10 

(ii)  the  phase  velocity  is  lowest  for  the  lowest 

k modes; 
l 

(iii)  as  v -*•  0,  8 ■*  8,  -*■  8 so  that  the  plasma 

<p  k o 

wave  phase  velocity  approaches  the  electron 
beam  velocity; 

(iv)  as  the  electron  current  is  increased  from 

zero,  the  phase  velocity  decreases  below  the 
beam  velocity,  although  throughout  the  regime 


Thus  8 . (v) 


is  only  a small  correction  to  8 


will  appear  fairly  flat  near  8 


! 


Since  the  lowest  phase  velocities  are  realized  in 
the  long  wavelength  limit,  let's  consider  the  case  of 
Jl  = 0 and  k -*-0.  We  shall  make  the  ansatz  that  6 

Z ( p 

remains  bounded  in  this  case,  or  that  in  -*•  0 as  well. 

Then  for  g(r)  and  b/a  arbitrary,  equation  (19)  reduces 
to 


L 

I 

I 

I 

I 

I 

I 

r 


d 4vg(r|<i - 

r dr  r dr  + y30  (0  - 0 . ) 2 

z z <p 


$ = 0 


while  the  boundary  condition  (21)  reduces  to 


(x  . b d£\ 

^ + _ r 

r=a 


= 0 


and 


4(r  = 0)  = 1 


(36) 


(37) 


(38) 


It  is  now  apparent  that  as  0X  0,  equations  (36)  - (38) 

<P 

for  4>  approach  the  limiting  current  conditions  (8)  - (10) 
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the  phase 


for  the  variable  T.  Consequently,  as  v -*■  v , 

L 

velocity  must  vanish  for  arbitrary  choices  of  the  density 
profile  and  b/a. 

Finally,  numerical  solutions  have  been  determined 

to  equations  (19)  - (22) , and  typical  results  are  illustrated 

in  Figures  1-7  for  cases  with  g(r)  = 1 and  l = 0.  The 

equilibrium  quantities  y(r),  6 (r)  which  appear  in  the 

z 

coefficient  of  equation  (19)  are  determined  from  a numeri- 
cal solution  to  the  equilibrium  equations  (5)  - (7)  . 

In  Figure  2,  the  dispersion  relation  u (k  ) is 

z 

plotted  for  the  lowest  k mode,  for  the  case  y =7  and 

x o 

b/a  = 1,  but  for  three  different  choices  of  y : (a)  y > 

c c 

Y T , (b)  y = y _ , (c)  y < y . These  three  curves  in 

cL  c cL  c cL 

Figure  2 may  be  correlated  with  the  three  points  marked  (a) , 

(b) , (c)  on  the  corresponding  v(yc)  curve  in  Figure  1. 

The  predicted  scaling  of  S vs.  k , v may  be  observed. 

<t>  z 

Equation  (35)  was  also  found  to  be  quantitatively  correct 
when  v <<  v . The  possibility  of  negative  phase  velocities 
for  y < Y is  possibly  associated  with  the  transition  to 

C CL 

virtual  cathode  formation  which  is  experimentally  observed 


19 


to  occur  after  v has  been  raised  to  v and  y reduced 

L c 


“ YcL- 


In  Figures  3,  4,  5,  and  6,  the  linear  eigenfunctions 

<f(r)  are  illustrated  for  the  two  smallest  values  of  k^  and 

for  two  values  of  v.  For  these  figures,  y = 7,  b/a  = 1, 

o 

and  k =0.  It  may  be  observed  that  the  lowest-k  modes 
2 X 

have  no  interior  nodes,  while  successively  higher-k^  modes 

will  have  successively  more  nodes.  The  eigenfunctions  for 

small  values  of  current  (i.e.,  v < v ) are  relatively  broad, 

L 

and  exhibit  the  JQ(k^r)  dependence,  while  the  higher 

current  eigenfunctions  are  more  peaked  near  the  beam  axis. 

Finally,  in  Figure  7,  3,  is  plotted  vs.  v/v  for 

y L 

two  choices  of  y . For  these  data,  k = 0 and  the  lowest 

o z 

k^-mode  is  examined.  The  two  end-points  of  such  a plot  are 

known  to  be  independent  of  b/a,  and  the  rest  of  the  plot 

was  also  found  to  be  remarkably  insensitive  to  b/a  for  the 

modest  values  of  b/a  which  were  computed  (i.e.,  1 £ b/a  £ 2). 

Therefore,  we  are  able  to  conclude  from  this  study  that  the 

difficulty  in  lowering  the  plasma  wave  velocity  seems  to 

persist  even  in  the  fat  beam  regime  where  the  radial  variation 

in  the  equilibrium  quantities  y,  3 , n is  significant.  It 

z 
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remains  necessary  to  increase  the  electron  current 


very  near 


THE  EFFECTS  OF  FINITE  THERMAL  SPREAD  ON  THE 
ELECTRON  CYCLOTRON  BEAM  MODES 


(Appendix  J Contains  16  Pages) 


T 


For  a cold  relativistic  electron  beam  propagating 

with  velocity  u -*•  c along  a magnetic  guide  field  in  the 

z-direction,  the  frequency  uj  of  the  electron  cyclotron 

beam  mode  with  wave  number  k is 

z 


0)  SS  k u - ft 
z 


+ 


a)2  n 

_E 

k2c2 


where 


k1 2 


k2  + k2 
z A 


2 


<< 


1 


£2  and  u)  are  the  relativistic  electron  cyclotron  and 
P 

plasma  frequencies  respectively.  The  perturbed  electron 


perpendicular  velocity  vx  associated  with  the  cyclotron 
beam  mode  is 


1 

V <*  ■ ■ — ■ ■- 

1 u - k v + 
z z 


1 


Thus  a thermal  spread  in  of  magnitude  Av^ : 


v u + Av 
z z 


such  that 


k Av 
z z 


k2C2 


will  modify  the  perturbed  perpendicular  velocities  and 

hence  the  linear  dispersion  relation  of  the  electron 

cyclotron  beam  mode  for  a cold  beam. 

A relativistic  electron  beam  will  acquire  a thermal 

spread  in  v after  passage  through  an  anode  foil.  If  the 
z 

foil  is  thin  enough,  the  scattering  of  the  electron  beam  by 
the  foil  results  in  a negligibly  small  spread  in  energy  but 
rotates  the  electron  velocity  vector  in  velocity  space.  A 
finite  spread  in  propagation  velocities  occurs. 

In  this  appendix,  we  investigate  the  effects  of  a 
finite  thermal  spread  in  v^  on  the  dispersion  relation 
for  the  electron  cyclotron  beam  mode. 


Equilibrium 


We  consider  a slab  equilibrium  of  the  electron 
beam,  infinite  in  extent  in  the  z-  and  y-directions  and 
bounded  by  conducting  walls  at  x = ± a.  The  beam  propa- 
gates along  a magnetic  guide  field  B in  the  z-direction 

oz 

The  electron  beam  equilibrium  is  described  by  the 
equilibrium  distribution  function 


f 

o 


<e,p  ) 

Z 


d 6 ( e - e ) g [p*(e,p  )], 
o o o z 


P > 0 
z 


0 


r P < 0 

Z 


(1) 


where  t and  P^  are  the  energy  and  z-component  of  the 
canonical  momentum  respectively: 


t =>  myc2  - e*  (x) 
o 


P 
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* A (X) 
c oz 
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/•— W*  rri'n  i J--  ' 


at 


7 T- 


4 


■f  SB 


.,/*  ^ 


n> 


P*(e,P  ) = — - m2c2  - P2 
x z c z 


<<  m2Y2c2  is  a measure  of  the  thermal  spread 


The  equilibrium  beam  density  is 


N =*  fd3p  f 
o o 


+00  00 
2ir  f dP 


f V I 
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— i f 

c2  o 


-oo  e (p  ) 

b z 
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Mm  “iip 


. + ? A.J  cJ  + m!o‘l 


1 / 2 


- € 


We  restrict  the  discussion  to  the  case  of  a low 


density  ultrarelativistic  electron  beam  where  e >>  e<l>  , 

o o 

eA  . In  this  limit 
oz 


N 
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2tt  D 
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J z 
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where 


e2 

o 


m2c2  >>  p2 


Equation  (2)  relates  the  normalization  constant  Dq  to 
the  beam  density  N^. 

Similarly,  the  equilibrium  current  is 


J 

oz 


e / d3p  — f 
J my  o 


* - e 


N p 

o o 

my 
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il 

I 

1 
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From  Maxwell's  equations,  the  equilibrium  electric 


potential  4>  and  magnetic  potential  A can  be 
o oz 


determined. 


The  equilibrium  electron  trajectories  may  be 


approximated  by 


where 


p'  - P 
z z 


p"  * p±  COS  {ft  (t'  - t)  + i|/} 


p'  * px  sin  {ft(t'  - t)  + <l>) 


z'  * z + — (t'  - t) 


x'  - x z — — [sin  {ft(t"  - t)  + i|)}  - sin 
tayu 


, r:r..  - 


■ - 

■ i ■■  ■ 

„■  if  - . — 





Dispersion  Relation 

i 

Let  the  perturbed  electromagnetic  fields  have  the 


form 


E(r,t)  -*■  E (x)  exp  [-  iuit  + i k z],  etc. 

z 


The  linearized  perturbed  distribution  function  is 


given  by 


r , „ 

e^bu')] 

dt 

E (x")  + — 

J 

—00 

— myc 

o: 


where  the  time  integral  is  taken  over  the  equilibrium 
trajectories. 

Neglecting  finite  Larmor  radius  terms  of  order 
P*)/(m2Y2a2)  <<  1, 


The  perpendicular  perturbed  currents  are 


(pv 

± i p ) 

I 

= - e / d3p 

f,  “ 
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+00 
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= - e j*  dtp 

1 dPz 

| de  Pl  e±M'  f1  (5) 
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• 09 

e (p  ) 

b z 

The  electron  cyclotron  beam  modes  have  frequency 
oj  ~ (k  P )/(mY)  - ft  so  that 


i 

1 

k P 

> > 

k P 

,,  z z + 0 

z z 

CO  — 

my 

my 

in  equation  (5), 


+ + + 
A = E"  ± i B~ 


The  solutions  of  these  equations  are 


+ + + 

A = a sin  k x + a cos  < x 

J.  <• 


A = ^ sin  k x + S ^ cos  k x 


where 


, cu 2 p2  (k  c - <d) 
oj2  - P o z ,,  Q 

k2  + — 2 + 2 2 r—— {1  - e 

z c 2 c2  k c 

P Z 


{1  - ew  Q E1 (Q)  } 
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,2  _ P 2 

z c2  c2  (k  c - oj  + ft) 
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The  boundary  conditions  on  A at  x = ± a where 


the  beam  is  bounded  by  conducting  walls  are 


0 


'W- 


+ 

A = A 


3 A 3 A 
3x  3x 


= 0 
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For  solutions  antisymmetric  in  A~, 
relation  is 


+ + — — 
k cot  < a + < cot  k a = 0 


If  < a <<  1,  <+a  cot  tc+a  «•  - 1,  and  <+a  * 

The  frequency  to  is  then  given  by  equation 
of  Q,  the  frequency  to  is, 


Q pZ 

o)  = k c - ft  + — k c 
z 2 pz  z 

o 


Q vz 

= k c - fi  + ^ — - k c 
z 2 c z 


In  the  limit  of  Q >>  1, 
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Ei(Q)  - — 


1 _2 
Q + Qs 


1 - r + rr  + 


r ■ ~y, 
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(13) 

the  dispersion 

(14) 
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(11) . In  terms 


(15) 


and  from  equation  (11) , 


2 P‘, 


0) 


(k2C2 


2 


The  frequency  uj  is 


<JL) 


a)  * k c - 8 + 

Z 


ft  - ~2  k c 
c z 


k2c2  *•  - „2 


where  k2  * k2  + k + 2 >> 


— and  we  have  assumed  that 


<*>‘  2 ui 2 

7T§T  <<  We  also  require  that  >>  1 


k2c 


k202  k c 
z 


for 


validity  of  the  asymptotic  representation  of  (Q) . 

A finite  thermal  spread  introduces  a correction 

* 2 
V 

Aw  = - —7  k c to  the  frequency  of  the  electron  cyclotron 
c z 

beam  mode. 

A 5 

, V2 

As  —7  increases,  the  magnitude  of  Q which  solves 

C 

equation  (11)  decreases.  Q remains  real  and  consequently 

k2  v2  o 

w is  real.  In  the  limit  of  — 5—  ■+■  2 , Q -►  0 and 

id2  k c 

P z 

the  frequency  u -*■  k c - (2.  The  neglect  of  finite  Larmor 

z 
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radius  effects  in  this  discussion  imposes  the  condition 

to 2 <<  ft2 . 

P 

In  conclusion,  if  the  thermal  spread  is  low  enough 
k2  v2  2 (2 

so  that  — 2 — K r~  and  finite  Larmor  radius  effects  can 

OJ  K C 

P z 

be  neglected  (k|02)/(fi2)  <<  1,  the  solutions  of  equation 
(11)  for  co  are  real,  and  the  electron  cyclotron  beam  modes 
propagate  undamped  by  thermal  effects. 
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